Abstract 

Jefferson Lab Experiment E04-001 used the Rosenbluth technique to mea- 
sure R = (Tl/cft and F2 on nuclear targets. This experiment was part of 
a multilab effort to investigate quark-hadron duahty and the electromagnetic 
and weak structure of the nuclei in the nucleon resonance region. In addition 
to the studies of quark-hadron duality in electron scattering on nuclear targets, 
these data will be used as input form factors in future analysis of neutrino data 
which investigate quark-hadron duality of the nucleon and nuclear axial struc- 
ture functions. An important goal of this experiment is to provide precise data 
which to allow a reduction in uncertainties in neutrino oscillation parameters 
for neutrino oscillation experiments (K2K, MINOS). This inclusive experiment 
was completed in July 2007 at Jefferson Lab where the Hall C High Momentum 
Spectrometer detected the scattered electron. Measurements were done in the 
nuclear resonance region (1 < W"^ < 4 GeV"^) spanning the four- momentum 
transfer range 0.5 < g2 < 4.5 [GeV'^). Data was collected from four nuclear 
targets: C, Al, Fe and Cu. 
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Introduction 

The study of the strongly interacting particles of nature is the study of hadrons. The 
hadrons make up the atomic nucleus and the proton and neutron account for almost 
all the mass in the known universe. In contrast to their abundance our knowledge of 
the force which determines their interactions - the force responsible for both holding 
the nucleus together, limiting its maximum size as well as determining its dynamics, 
is the least well understood of the four fundamental forces of nature. This should 
not be surprising given the description of the strong interaction in terms of Quantum 
ChromoDynamics (QCD) wherein the nucleons are the lowest mass excitations of the 
complicated blend of quark and gluon condensates which form the QCD vacuum. The 
key features associated with the QCD Lagrangian are exhibited by the nucleon are: 
color confinement, asymptotic freedom and spontaneously broken chiral symmetry. 
A study of the nucleon is a study of the strong interaction and QCD. QCD is a 
renormalizable field theory based upon the principle of local gauge invariance under 
the exchange of color. QCD describes the strong force in terms of fermion fields of 
a given color charge, given the name quarks by Gell-Mann, interacting through the 
exchange of massless gauge bosons known as gluons. The distinctive feature of QCD 
which sets it apart from the electroweak theory is the fact that these gluons possess 
non-zero color charge (unlike the photons in electro-magnetism which carry no electric 
charge) . 

Among the main consequences is the property of color confinement - quarks and 
gluons are not permitted to exist as isolated free particles, but must combine to form 
color-neutral singlets: the three-quark baryons, of which the nucleons are the most 
stable, and the quark/anti-quark mesons. To date, the quantum numbers of all known 
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strongly interacting particles have been accounted for using a model based upon the 
non-Abelian internal SU (3)c symmetry of QCD for elementary quarks of six different 
flavors (these are labeled in order of increasing mass as up, down, strange, charm, 
bottom and top) . Though the static properties of the nucleon can be described in this 
manner, the extreme non-linear nature of the strong force means that the dynamical 
properties of its constituents vary significantly depending on the momentum scale at 
which it is studied. Practitioners are forced to confront a profound chasm: large mo- 
mentum behavior of the nucleon can accurately be described by the quark and gluon 
fields of which it is composed, but this same description fails at low momentum where 
successful models have to rely instead on effective hadronic degrees of freedom. This 
is related to the asymptotic freedom of QCD, responsible for the fact that the force 
between the quarks within the nucleon becomes weaker as they move closer together. 
For high momentum (hard) scattering processes involving the nucleon, perturbative 
QCD (pQCD) can be used to express the physical amplitudes as a perturbative series 
in the strong coupling constant cCs, analogous to Quantum Electrodynamics (QED). 
Renormalization of the field theory is required for such a series to converge and is 
achieved by the introduction of an arbitrary mass scale. This, in turn, leads to the 
running of the coupling strength of the strong interaction. Said another way, the 
effective strength is not constant as in the electromagnetic case but depends on the 
momentum scale, decreasing logarithmically at large momentum transfers. Hence an 
increase in momentum not only implies a shorter distance scale but a decrease in the 
effective coupling, allowing the perturbative expansion to converge and enabling rig- 
orous predictions based upon the quark and gluon fields alone. Asymptotic freedom 
also means that the relatively large running coupling strength associated with the 
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low momentum (soft) regime makes a perturbative expansion in of little value at 
large distance scales, such as typical hadronic sizes or larger. Rigorous solutions are 
not available for non-perturbative QCD and nucleon behavior at comparatively low 
energies leaving effective theories the only viable alternative. These are based upon 
identifying appropriate hadronic degrees of freedom for a particular momentum range 
and determining their characteristics. 

A variety of methods, each limited in their applicability, have been in use for 
some time: QCD sum rules [U El |3l HI Ej rests on the exploitation of the underlying 
symmetries of the strong interaction; effective Lagrangians are the basis for many 
model which approximate QCD at low energy; computationally intensive calculations 
of physical quantities on a discretized space-time lattice are becoming commonplace. 

Common to all these approaches is the special role played by the non-perturbative 
vacuum - it links low energy QCD to the phenomenon of spontaneously broken chi- 
ral symmetry. Light quarks are approximated as massless (when compared to the 
nucleon) and the conservation of their left and right-handedness introduces another 
layer of symmetry into the effective Lagrangian. This chiral symmetry is not exact 
and can be spontaneously (or dynamically) broken by the non-zero mass of confined 
quarks, leading these quarks to develop large effective masses through interactions 
with the vacuum. The result is that a nucleon consists not only of three valence- 
quarks but includes a complex sea of quark/anti-quarks pairs and gluons. This sea's 
characteristics are determined by the momentum scale involved. The electromagnetic 
response of the nucleon has been studied for more than fifty years now with measure- 
ments largely focused towards to the extreme upper and lower ends of the spectrum 
of momentum transferred to the nucleon. The intermediate region involving moder- 
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ately large momentum transfer, where there is significant give and take between soft 
and hard physics and therefore between quark-gkion and meson-baryon degrees of 
freedom, has long since been identified as a prolific testing ground for models of the 
strong interaction. 

Deep inelastic scattering (DIS) has proved to be a powerful tool to study the 
structure of matter. The experiments at SLAG O |71 E] showed for the first time the 
absence of scale dependence in inelastic electron-proton scattering. Later experiments 
at larger four-momentum transfer revealed logarithmic scaling violations, which was 
explained in terms of QCD and proved QCD as being the correct theory of strong 
interactions. In a few GeV energy range where hadronic degree of freedom is dominant 
the strong coupling constant becomes large and pQCD becomes inapplicable. This 
energy range is where the effect of confinement (no free quarks) make strongly-coupled 
QCD highly non-perturbative and it is more easy to work with mesons and baryons. 
Despite the significant difference between low and high energy description of the 
nuclear dynamics a phenomenon of so called quark-hadron duality is observed which 
can be though as a link between this two regimes. The quark-hadron duality was 
first observed in data from the early SLAC experiments by Bloom and Oilman [9]. 
The observations showed a striking similarity between the F2 structure functions 
measured in resonance and DIS. Based on this data Bloom and Oilman proposed 
that the resonances are not a separate entities but are an intrinsic part of the scaling 
behavior of F2. In order to analyze the degree of similarity of the F2 scaling function. 
Bloom and Oilman determined a Finite Energy Sum Rule [TT] which allowed one 
to investigate the phenomenon of quark-hadron quantitatively. 

With the development of QCD in the early 1970s, Bloom-Oilman duality was 
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reformulated [12] in terms of an operator product expansion (OPE) of moments of 
structure functions. This has allowed an expansion of the structure function integrals, 
called moments, in terms of the hard scale clarified many aspects of quark- 

hadron duality and it's violation. However, the OPE is not able to adequately explain 
why at low energies the correlations between partons were suppressed, and how the 
scaling worked for the resonances. 

The advent of high energy, high duty-factor electron accelerators, like the 6 GeV 
electron accelerator at Jefferson Lab, has provided quality data in recent years to 
test and constrain the theoretical calculations in this field. One of the significant 
observation in the JLab data is that Bloom-Gilman duality appears to work at 
values as low as 1 GeV^ or even lower. At low the strong coupling as is relatively 
large compared to it's value in DIS, but on average the inclusive scattering process 
appears to mimic the scattering of electrons from nearly free quarks. These observa- 
tions led to renewed interest in quark-hadron duality. In particular, measuring the 
F2 structure function in the resonance region at high Bjorken x {x > 0.7, where x is 
the longitudinal momentum fraction of the hadron carried by a parton in the infinite 
momentum frame) were of great interest since it allowed precision study of the F2 
structure function in x > 0.7. 

The observation of duality in electron scattering is thought to be a fundamental 
property of quarks and hadrons, which raises the question of the existence of du- 
ality in neutrino scattering. This is expected because scattering takes place from 
the same constituents as in electron case. Weak scattering can provide complemen- 
tary information on the quark structure of hadrons, not accessible to electromagnetic 
probes. Neutrino-induced reactions also can provide important consistency checks on 
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the validity of duality. Since the weak interaction mechanism is different from that of 
electromagnetic one the observation of duality may provide additional information. 
The MINERvA [13] experiment intends to study duality in neutrino scattering and 
can provide data to answer at what kinematic range duality works, in what structure 
functions and in what reactions. 

Duality in electron scattering have focused mostly on proton. There have been 
some experiments that performed measurements on deuterium and heavy nuclei in the 
high-x and low to moderate region [T5l [T6] . The results of these experiments 
revealed additional information about duality. Scaling in nuclei was observed and 
was interpreted in terms of local duality where the averaging over the resonances was 
accomplished by the Fermi momentum of nuclei inside the nucleus. 

The primary motivation of this experiment is to study quark-hadron duality in 
electron scattering from nuclear targets in the resonance region as well as to study the 
nuclear dependence oi R = olI ot- A significant nuclear dependence of R is predicted 
by Miller [17] due to nuclear pions at low values of ~ 0.3 GeV^. The results of 
this experiment, combined with the results of another experiment, E02-109 [18], will 
allow the extraction of through a precise Rosenbluth separation [19j and test if 
the pions are really the carriers of the nuclear force as predicted by Miller. Also, 
the data will allow us to model the electron— nucleon scattering cross section in the 
resonance region which are of great importance for neutrino scattering experiments 
since the neutrino scattering structure functions can be extracted if the vector part 
is provided from electron scattering. 



Chapter 1 
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Nucleon Structure 

In this section we provide a simple review of the progress made in understanding the 
nucleon, introducing the necessary formalism and key concepts along the way. An 
excellent overview of the field of nucleon structure can be found in Ref. [20] . 

1.1 Inclusive Lepton-Nucleon Reactions 

Inclusive lepton-nucleon scattering is an important tool with which to study the struc- 
ture of the nucleon. In the energy range accessible to modern accelerators leptons are 
considered structure-less. Therefore leptons are ideal to probe the the nucleon since 
the only unknown is the structure of the nucleon. In inclusive scattering, eA^ — )■ e'X, 
only the final state electron, e', is registered while the hadronic final state X remains 
undetected. In the target rest frame (laboratory frame) the incident electron with 
energy E scatters from the target, the scattered electron have angle 6 relative to the 
direction of the incident beam and energy E'. In the one photon exchange (Born) 
approximation, shown in Fig. [T| the scattering takes place by exchanging a virtual 
photon (or or Z boson in neutrino scattering). The energy of virtual photon is 
given hj u = E — E' and the 3-momentum transfer in q. Since the photon is spacelike 
the virtuality of photon = — g* ^ is negative. It is convenient to use positive 
variable = — g^, which is related to the initial and final energy and angle of the 
electron. 

Q^ = AEE'sm^^ (1.1) 
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Figure 1: Leading order Feynman diagrams for (a) elastic ep scattering and (b) 
deep inelastic scattering. In each respective case, the unknown aspects of the proton 
structure are parametrized in terms of the form factors, Fi[Q'^) and F2{Q'^), and 
parton distributions, q{x) and Aq{x). 

The invariant mass squared of final hadronic state is given by 



In general the cross section for inclusive scattering depends on two variables u and 
Q^. In the Deep Inelastic Scattering (DIS) regime, which will be discussed later, yet 
another important variable Bjorken x 



is introduced. In the DIS region this variable is the fraction of nucleon momentum 
carried by the struck parton in the "infinite momentum frame" — the frame where the 
electron is in rest while the nucleon is speeding toward it. For the nucleon < x < 1. 
Since the electron-nucleon scattering is electromagnetic in nature the one photon 



(1.2) 



9 

exchange approximation is a very good approximation. Therefore the differential 
cross section for electron nucleon scattering from unpolarized target can be written 
in as 

-^^^^L^'^W (14) 

where a is the fine structure constant. The leptonic tensor L^j, is calculable from 
Quantum Electrodynamics and can be written as 

— 2 [h^hy + h^J^n — Qjxuk ■ k ) (1-5) 

where k and k' are initial and final electron momentum, respectively and g^i, is the 
metric tensor. For an unpolarized initial nucleon the general form of hadronic tensor 
can be written as 



W^'' = W.g'^^ + ^p^p^ + ^q^q^ + ^ (p'^g^ + gV) (1-6) 



where Wj are unknown functions of u and Q"^. The structure function 14^3 and Wq 
do not contribute in the electron scattering, since Ws violates parity conservation 
and We's antisymmetric component is multiplied by a symmetric lepton tensor L/^^. 
Based on Lorentz and gauge invariance, together with parity conservation in electron 
nucleon scattering the hadronic tensor can be written as 

H^- = W.W, Q^) - r) + {p^ + '^,") + '^f) (1.7) 

where Wi and W2 are scalar functions of i> and Q^. After contracting L^j^yW^^ the 
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cross section can be written as 



dQdE' 



O'Mott 



W2{u, Q'') + 2Wi{u, g2)tan2 



:i.8) 



where aMott is the Mott cross section for the scattering from a point particle of charge 
e and is given by the following formula 



(7 Mott 



[1.9) 



The structure functions Wi and W2 contain all the information about the struc- 
ture of the nucleon and are measured by experiment. These can also be written as 
dimensionless functions as 



Fi {x, = MWi {u, Q^), and 



;i.io) 



F2{x,Q^)^iyW2{u,Q^ 



(1.11) 



The two structure functions Wi and W2 are related to the photoabsorption cross 
section. This is due to the fact that the virtual photon has two states of polarization, 
transverse (helicity ± 1) and longitudinal (helicty 0). Based on this the inclusive 
cross section can be written in terms of ar and ctl, 



dQdE' 



:i.i2) 
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where F is the flux of virtual photons and is given by 



aK E' 1 



;i.l3) 



The factor K is given by 



2Mu - 
2M 



(1.14) 



and is defined as "equivalent photon energy". In other words, K is the energy that 
would be required to form the final hadron state of the same mass that would be 
created by the real photon with energy v. The ratio of longitudinal to transverse 
virtual photon polarization given by 



1 -1 



(1.15) 



and ranges between e = and 1. 

The structure functions Fi and F2 can be written in terms of ar and ctl as 



KM 

Fi = — ^(Tt, and 



;i.i6) 



ATT^a (1 + v^lQ^) 



The ratio of longitudinal to transverse cross section is given as 



The structure function Fi is purely transverse as it is related to gt only, while the 
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F2 is a combination of both and a^- A purely longitudinal structure function 
is defined as 

Fl= {l+'^^F2-2xFi (1.19) 

and R can be written as 

i?=^%. (1.20) 

Using the ratio i?, the F2 structure function can be extracted from the measured cross 
sections a according to 

a l+R , , 

-ye-—— 1.21 



c" Mott 1 + e-R 

provided that R is already known. 

Inclusive neutrino-nucleon scattering takes place by the exchange of W bosons. 
These are the charged current reactions involving neutrinos: vN — )■ e~X or 
vN — > e^X. The cross section can be written as |2T] 



dndE 

where mw is the mass of the W and the tensors correspond to the same definition as 
in Eq. 1.4 After contracting the L^^}W^^ the differential cross section is written as 



,^ ,^ - o 1 9 I 1 21^1 sin^ - + W2 cos^ - =F PV's — 7-^ sin 

(1.23) 

The difference between electromagnetic and weak scattering is the existence of the 
third structure function. This is a result of parity violation in weak interactions. 
The structure functions of neutrino scattering can also be written in terms of the 
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photo- absorption cross sections. 



K 



(1.24) 



W2^ 



K Q2 



^ {(Tr + (7l + 2(7s) , and 



(1.25) 



K 2M 



(1.26) 



W^3 = 



In the case of electron scattering parity invariance forces ctr = ctl (L and R stand 
for left and right, S stands for longitudinal polarization), so there was only one 
term gt- One of differences between electron and neutrino scattering is that the 
electron— nuclear structure functions have contributions only from vector- vector (VV) 
terms, while the neutrino structure functions are a combination of W and axial-axial 
(AA) terms. The vector-axial vector interference term is given by {W^ = in 
electron scattering since (Tr — ul)- The significance of comparing neutrino scattering 
to electron scattering will be discussed later in next section. 

1.1.1 Deep Ineleistic Scattering and Queirk Parton Model 

Elastic electron-nucleon scattering happens when the struck nucleon recoils without 
disintegration and the final hadronic mass W is equal to the nucleon mass, W — Mj^. 
Nuclear resonances are excited states of the nucleon and can be seen in the spectrum 
of inelastic electrons when W <2 GeV and < 4.5 GeV^. At fixed W nuclear the 



structure functions go to zero when — )■ oo, 
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^0. 



;i.27) 



If instead of keeping W fixed, the variable x (defined in Eq. 1.3) is kept constant for 
> \ GeV^ and W > 2 GeV the remarkable phenomenon of Bjorken scaling [22J 
is observed, see Fig. [2j The structure functions Wi and W2 become independent of 
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Figure 2: World data on the structure function F2{x,Q'^) of the proton from the 
Particle Data Group 
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the mass scale. It was the measurement of Wi and W2 at SLAC [S] that uncovered 
the first evidence that the nucleon consists of structure-less particles, and the cross 
section was an incoherent sum of individual elastic scattering cross sections from these 
constituents. The data showed that at values above a few GeV, the structure 
functions depend only on this new variable x: 

m^W,ix,Q')^Fnx), 

uW2ix,Q')^Fr''{x), (1.28) 

Recalling that the structure function for any object with internal structure must 
be Q^-dependent, this observation implies that the virtual photons in DIS must be 
scattering from point-like, structure-less objects inside the proton. These internal 
particles were given the name of partons, it was only later that they were associated 
with the quarks and gluons of QCD. The Q^-independent behavior is the property 
of Bjorken scaling and arises as a consequence of the asymptotic freedom associated 
with QCD at short distances and imphes independence of the absolute resolution 
scale and hence point-like substructure. The physical significance of x becomes clear 
in the infinite momentum frame [21] where nucleon is moving with a momentum 
approaching to 00 in the z-direction. In such a frame, relativistic time dilation implies 
that during the interaction with the virtual photon, the nucleon can be considered as 
a collection of non-interacting partons (the impulse approximation [2.1]), each with 
different fractions of the total nucleon longitudinal momentum. For one of these 
partons to absorb a virtual photon of energy u and mass Q^, it must carry exactly the 
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momentum fraction x. It is then clear that DIS, in addition to providing evidence for 
internal structure of the proton, can be used to measure the momentum distribution 
of the partons inside the nucleon (for a review see Ref. pl|). The unpolarized parton 
distribution q{x) for both spin-1/2 and spin-0 partons can be related to the nucleon 



structure functions of Eq. |1.28| by the relation: 

F^{x) = xJ2eV^{x) (1.29) 

i 

where the sum is over the various species of partons of charge ef and fi{x) is the 
probability that the parton has momentum in the interval x — )■ x + dx. The Fi{x) 
structure function is related to parton distribution functions by the following formula 

Frix) = F,ix)/2x = lY.^U^{x). (1.30) 

i 

This relation is known as Callan-Gross relation. From Eq. |1.19 and Eq. 1.20 the R 



defined in Eq. |1.20 can be written as 



The above equation shows that when — )■ oo, i? — ?■ 0. Since the virtual photons 
have helicity ± 1, they can only be absorbed by partons with non-zero spin (this 
follows from conservation of angular momentum) otherwise after absorbing a photon 
the parton will still have zero angular momentum, violating conservation of angular 
momentum, since the total angular momentum of initial state was not zero). This 
prediction, that R is small at large is in agreement with data |25l [26], |27l |28l [29] , 
but the predicted quadratic reduction of R at small x is not observed. R being small 
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supports the idea of partons having spin ^. 

All of the above discussion assumes that partons don't have mass or transverse 
momentum. Including these effects [30] R can be written as 



where is transverse momentum of partons, m is the mass of partons, A is binding 
energy of partons. This form is significant since it predicts large values for R in 

=1-5 GeV range. A very nice review about QPM can be found in Ref. [30] . 

The quark and gluon parton distributions have been accurately measured over 
a wide range of kinematics, and have led to important revelations concerning the 
internal nature of the nucleon. The most significant of these discoveries is the fact 
that the sum of quark momenta only amounts to around 45 % of the nucleon's total 
momentum. If instead of considering the spin-averaged structure functions one intro- 
duces the spin-dependent structure function gi, it is possible to define the polarized 
parton distributions Ag(a;) and Ag(x) in a similar way: 



The polarized parton distributions are defined as the difference between quark dis- 
tributions with spin parallel/anti-parallel to the proton spin Ag(x) = q{xY — q{xY. 
Measurement of these entities led to another discovery: the spin of the quarks only 
contribute about 25 % to the total spin of the nucleon. This spin crisis was discovered 
at CERN in the EMC experiment [31] and has led to the question: if the quarks only 
carry a quarter of the nucleon's intrinsic spin, where does the remainder come from? 



R 



4/c|, + rri 



2± A 



(1.32) 




(1.33) 
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To answer this, the nucleon spin Sn can be decomposed into three components: 

5^ = ^AS + AG + L, = ^ (1.34) 

where AS is the total spin of the quarks, AG the total spin of the gluons and is 
the quark's orbital angular momentum. Recent results of DIS data [32] showed that 
contribution of AG is compatible with zero. 

1.1.2 Scaling Violation in DIS 

In the simple quark-parton model partons are assumed to be asymptotically free at 
high Q^. In addition if one assumes that number of partons do not change with the 
scaling should not be violated at all. However, it can be seen from Fig. [2] that scaling 
is violated since function F2 has some dependence. The observed scaling violation 
is due to interactions between partons and a change in their numbers. At low values 
of X the structure function F2 rises as increase while at large x it decreases. The 
scaling violation at low x is caused by gluon splitting to quark-antiquark pairs. As a 
consequence, the number of quarks with low values of x increase, which is equivalent of 
F2 increasing as increase. At high values of x where valance quarks are dominant, 
gluon radiation shifts quarks to lower x causing F2 to fall with increasing Q^. 

This scaling violation can also be explained in terms of the resolution effects 
of the virtual photon. At some value of virtual photon can probe structure of 
A ~ Xj \fQ\. At low a;, where sea quarks and gluons carry most of the momentum 
of the nucleon, only part of them can be seen by the virtual photon. At Q\ > Q\ 
resolution of the virtual photon is higher and is able to resolve more quarks and 
gluons increasing the value of F2. At the same time, the higher resolution of virtual 
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photon at high x (valance quarks are dominant) means that instead of seeing quarks 
with high X the virtual photon can see quarks which radiate a gluon and have their 
X lowered by some amount, effectively decreasing F2. 

The structure functions are not calculable in the framework of QCD. However, 
it is possible to calculate the evolution of structure functions in perturbative 
QCD. This evolution is done according to Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
equations (DGLAP) [531 [MIES]. 

Besides scaling violations caused by the self-coupling of gluons, which has a log 
dependence, there exist scaling corrections of the form (l/Q^)". These are called 
power corrections. These power corrections to perturbative QCD are due to the non- 
vanishing target nucleon mass and higher twist effects. The higher twist effects are 
caused by the interactions between the struck quark in the electron-quark scattering 
process and the other quarks in the nucleon (Final State Interactions). At large 
Q^, the target mass and higher twist effects are negligible. But, these effects may 
be significant at low Q"^. In order to study QCD scaling violations in F2 structure 
function Nachtmann introduced a scaling variable 



As it can be seen from the expression of Nachtmann scaling variable, as increase 
^ — )• X. Using the Nachtmann scaling variable ^ instead of x the scaling of the F2 
structure function can be extended to lower Q^. The lower limit of at which the 
QCD is still applicable is not exactly known. However, a comparison of the data 
[37] to a QCD analysis suggests that evolution work well for values as low as 




(1.35) 
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1.2 F2 and R in the Nuclear Resonance Region 



In the nucleoli resonance region < 4 and = 0.5-4.5 GeV^, the non-perturbative 
higher twist and target mass effects may be quite significant. This is expected since 
the resonances are bound states of quarks and gluons and necessarily involve higher 
twists. Therefore, the behavior of F2 and R in the resonance region may be different 
from that in the DIS region. Even at high values of Q^, R may not be small due to 
strong gluon binding effects (see Ref. [30]). 

The nuclear dependence of R in the resonance region has never been determined 
experimentally. The only data that exists is in the DIS region. In Fig. [3] the world 
data on R is shown for all nuclear targets. The left plot in Fig. |3] shows the world's 
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Figure 3: Left: A compilation of the world's published data on R at high in 
the DIS for both nucleons and nuclei. Right: The SLAG E140 data on the nuclear 
dependence of R in the DIS region presented in the form Ra — Rd (The magenta lines 
in the figure on the right represent the projected uncertainties of this experiment.) 



data on R for all nuclear targets in the DIS region. The right plot in Fig. |3] shows 
all available data on the nuclear dependence of R in the DIS region. The errors are 
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very large and none of the data are in the resonance region. Note that the data in 
Fig. [3] (right) were taken at high values of W"^ in the DIS region. At these high W and 

values, the overall value of R is generally quite small. Therefore, it is difficult to 
discern any nuclear effects. The measurements of the present experiment are at low 
Q^, where R is larger and a nuclear dependence should be easier to discern. The lines 
in Fig. [3] (right) represent the small projected uncertainties of the current experiment. 

Precision measurements of R in the resonance region at moderate to low will 
greatly aid efforts to develop a reliable global description of R. The new data on R, 
F2, Fi, and in the resonance region will also help efforts to develop global models 
to describe all unpolarized cross sections, covering the range from large x, low to 
smaller x, low Q^. (These models are useful for electron— nuclear scattering model 
development, and the extraction of parton distribution functions.) A detailed study 
of F2 and R = ctl/<^t on nuclear targets in the resonance region is an important 
ingredient in forming an integrated description of charged lepton and neutrino scat- 
tering cross-sections. High rate neutrino beams now under construction or planned at 
Fermilab and J-PARC will allow the first precision experimental comparisons of elec- 
tron and neutrino cross sections and present and planned future neutrino oscillation 
experiments will use these results to predict event rates. 

The full program of studies requires first additional precise electron scattering 
data, in particular ctl and ar (or equivalently F2 and R) on nuclear targets (ma- 
terials suited for future neutrino oscillation detectors, including water [38J, hydro- 
carbons [39j, liquid argon and steel), where the most precise high energy neutrino 
cross sections have been measured ^]) in the relevant kinematic regime. Later, as 
the new generation of high rate neutrino beams at Fermilab and J— PARC become 
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available, the approach can be directly validated with comparisons to data from high 
rate neutrino cross— section experiments on the same targets |13j . 

Additionally, R will be used to investigate the phenomenon of quark— hadron dual- 
ity (discussed in the next section), which suggests that, on average, electroproduction 
structure functions in resonance region scale the same way as in the DIS region. 
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1.3 Quark Hadron Duality 

At high energies perturbative QCD methods fully describe experimental results and 
at low energies chiral perturbation theory allows the study of low-energy dynamics of 
QCD. A phenomenon of quark-hadron duality, first discovered by Bloom and Oilman, 
is observed connecting these two regimes. In intermediate kinematics, a wide variety 
of reactions can be described simultaneously by single particle (quark) scattering, and 
by exclusive resonance (hadron) scattering. 

1.3.1 Bloom-Gilman Duality 

In the early 1970s, when analyzing inclusive electron-proton scattering data from 
SLAC, Bloom and Oilman observed [HI SI] a remarkable connection between structure 
function vW2{i'i Q^) in the nucleon resonance region and in the deep inelastic region. 
It was found that the resonance structure function averaged over the scaling variable 
u' was approximately equal to deep inelastic structure function. The scaling variable 
is given by the following formula 



, 2Mu + M^ , , 

^= =^^Q^=^^Q^ (1-36) 



where M is the nucleon mass. The choice of u' over u = 1/x was not clear at the 



time it was used. If one looks at the Nachtmann variable defined in Eq. 1.35 it will 
be clear that using u' for scaling parameter is somewhat equivalent of using 1/,^. 
This new variable u' was found without theoretical guidance but an extraction of the 
Nachtmann variable revealed a similarity. The Bloom and Oilman scaling variable 
u' included additional dependence that incorporates the finite target mass effects. 
Using the variable u', Bloom and Oilman were able to make the first quantitative 
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observations of quark-hadron duality in inclusive electron scattering. 

The original data on the proton vW'^{v^Q'^) structure function in the resonance 
region are illustrated in Fig. |4] for several values of from 1.5 to 3.0 GeV^. In 
the resonance region there are about twenty nucleon resonances, but only three of 
them are visible in measured inclusive electron-proton cross section. The first peak 
corresponds to a single resonance P33(1232), while other resonances are composed of 
overlapping states. The second resonance region comprises from the Sii(1535) and 
Di3(1520) resonances. The data presented in Fig. |4] are from inclusive measurements, 
and may contain tails from heavier resonances as well as some nonresonant compo- 
nents. The vW'^{v,Q'^) structure function data were extracted from the measured 
cross sections using a fixed value of the longitudinal to transverse cross section ratio, 
R = olIot = 0.18. 

The scaling curve shown in Fig. |4] is a parametrization of the high-W (high-Q^) 
data available in the early 1970s [12]. It can be seen that the resonance data are clearly 
seen to oscillate about, and average to, the scaling curve. Based on the similarity of 
the structure functions in the resonance and DIS region Bloom and Gilman concluded 
that the resonance data are, on average, equivalent to the scaling curve. Also, the 
agreement of resonance region data gets better with increasing Q^. These observations 
led Bloom and Gilman to conclude that resonances are not a separate entity but are 
an intrinsic part of the scaling behavior of vW2{i'-, Q^) [H]. 

In order to quantify the similarity of scaling functions. Bloom and Gilman wrote 
a finite energy sum rule jH] at fixed for i'W2{i^,Q'^) ■ They equated the integral 



26 



0.4 

as 

^ 0.2 

0.1 


0.4 

Q3 

0.1 





q2 = l.50 



l/W; 



0.4 

0.3 

'2 0.2 



0.1 - 



i l / 1 I I I \ I I 



= 1.75 





0.4 

0.3 



q2 = 2.75 




2 0.2 



0.1 



1.0 2.0 3.0 4.0 5.0 60 7.0 8.0 9.0 lOO 



q2=3.00 




1.0 2.0 30 4.0 5.0 6.0 7.0 8.0 9.0 10.0 



Figure 4: Early proton vW"^ structure function data in the resonance region, as a 
function of w', compared to a smooth fit to the data in the scahng region at larger 
Q^. The resonance data were obtained at = 1.5, 1.75, 2.75, 3.0 GeV^, for the 
longitudinal to transverse ratio R = 0.18. (Adapted from Ref. |i41j.) 



over u of in the resonance region, to the integral over u' of the scaling function: 



2M 



uW2{uj')duj'. 



;i.37) 



In this equation the upper limit on the z/ integration, Prn = (W^m ~ + '5^)/2^, 
corresponds to the maximum value of = 1 + W^/Q"^, where Wm ~ 2 GeV, so that 
the integral of the scaling function covers the same range in u' as the resonance region 
data. The finite energy sum rule allows the area under the resonances in Fig.|4]to be 
compared to the area under the smooth curve in the same u' region to determine the 
degree to which the resonance and scaling data are equivalent. 



A comparison of both sides in Eq. 1.37 for Wm = 2 GeV showed that the relative 
differences ranged from ~ 10% at = i GeV^, to ~ 2% beyond = 2 GeV^ [31]. 
This demonstrates the equivalence on average of the resonance and deep inelastic 
regimes. Using this approach. Bloom and Gilman's quark-hadron duality was able to 
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qualitatively describe the data in the range 1 < < 10 GeV^. The resonances in 
inclusive inelastic electron-proton scattering do not vanish with increasing relative 
to the "background" beneath them, but instead fall at about the same rate with 
increasing . The nucleon resonances are therefore strongly correlated with the 
scaling behavior of vW2- 

1.3.2 Duality in Nuclei 

Duality studies have been studied mostly on protons. There are handful of experi- 
ments in the high— x and low to moderate region that tried to study duality on 
deuterium and heavy nuclei [T^ |T3]. Inclusive electron-nucleus experiments at 
SLAG designed to probe the x > 1 region in the structure function concluded 
that the data began to display scaling indicative of local duality [H]. An explanation 
for the origin of ^ scaling was proposed by the Benhar and Liuti |13]. The authors 
suggested that the observed scaling might instead come from an accidental cancella- 
tion of dependent terms, and would occur only for a limited range of momentum 
transfers (up to ~ 7 GeV^). 

In order to study this effect and confirm the observation of the experiments done 
at SLAG an experiment was done at Jefferson Lab [44j. A plot shown in Fig. |5] 
is shown to illustrate the results of that experiment. Here, F^^ j A is plotted as a 
function of ^ at several values of [x = 1). An important thing to realize from the 
first glance is the absence of resonance structure, which is clearly observed for the 
free nucleon. Also, the quasi-elastic peak is not visible. The absence of a distinct 
quasi— elastic and resonance structure is a result of Fermi smearing. It was discussed 
in the previous section that structure functions averaged over the nuclear resonance 
region resembles the structure function of DIS. In the case of nuclei, it appears that 
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Figure 5: The vW2'^ = F2 structure function for iron (per nucleon) as a function 
of ^. The data were obtained at fixed electron scattering angle, and the quoted 
(in units of GeV^) are the values for x = 1. The arrows indicate the values of ^ 
corresponding to the quasi-elastic peak for each setting. (Figure from Ref. (TB]) 



averaging is done by Fermi motion of nucleons inside the nuclei. As the resonance 
structure is smeared by Fermi motion the scaling can be observed at all ^ and it 
is impossible to differentiate DIS and resonance regimes other than by calculating 
kinematics. It can be seen from Fig. [s] that all data except the lowest fall on a 
smooth scaling curve. 

Qualitatively, the nuclear effects in the resonance region appear to be similar to 
those in the deep inelastic region. The nuclear dependence of the scaling structure 
functions are not expected to be the same as the nuclear dependence of resonance 
production. This is somewhat surprising, but perhaps it is another hint that quark- 
hadron duality can be applied in inclusive nuclear scattering too. Based on the data 
of ^-scaling in nuclei, it can be implied that Fermi motion is averaging the nucleon 
electromagnetic response over a finite energy range (for the proton this was done 
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using Eq. 1.37). 

There is limited data available in the nuclear resonance region and at the same 
time the heavy model dependence at large x does not allow duality studies for nuclei 
at the same level as it is done for proton. The results of the current experiment will 
complement the existing data and allow precise studies of the nuclear dependence of 
quark-hadron duality. 

1.4 Theoretical Basis of Duality 

At the time when Bloom-Gilman duality was observed QCD was not yet fully de- 
veloped. In order to give theoretical explanation to it, phenomenological models or 
models based on hadronic degree of freedom were developed. After the development 
of QCD and the recognition that it is a real theory of strong interactions, the phe- 
nomenon of duality was reanalyzed using the Operator Product Expansion (OPE). 
This is discussed in the next section. 

1.4.1 Moments of Structure Functions 

Before starting a discussion of the Operator Product Expansion let's introduce 
Cornell-Norton [13] moments first. The n-th moments of the spin-averaged Fi, F2 
and Fl structure functions are defined as: 



M^{Q^)= [ dx x''-^Fi{x,Q^) (1.38) 
Jo 



MIl{Q'')= I dx x^-''F2,l{x,Q^). (1.39) 
Jo 
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In this definition the n = 1 moment of the Fi structure function in the parton 
model counts quark charges, while the n = 2 moment of the F2 structure function 
corresponds to the momentum sum rule. In the Bjorken limit, the moments of the 



Fi and F2 structure functions are related via the Callan-Gross relation, see Eq. 1.30 



as Mg""* = 2M["'\ The Cornwall-Norton moments defined in terms of the Bjorken 
X scaling variable are appropriate in the region of kinematics where is much 
larger than typical hadronic mass scales, where corrections of the type M^/Q"^ can 
be neglected. In the next section it will be clear how the structure function moments 
can be utilized to test the validity of quark-parton duality at small where the 
contributions from higher twist effects are not negligible. 

1.4.2 Operator Product Expansion 

The theoretical basis for describing Bloom-Gilman duality in QCD is the operator 
product expansion (OPE) of Wilson [46j. Based on OPE, sum rules are derived which 
can be used to analyze the structure functions and their moments. The important 
feature of these sum rules is model independence. They are derived using only some 
very general results from quantum field theory. 

The OPE allows one to evaluate a products of operators in the asymptotic limit. 
The product of two operator- valued fields, A{x) and B{y), can be expressed as an 
infinite series 

A{x)B{y) = Cn{x - y)On{x) (1.40) 

n 

where C„(x — y) are analytic functions of x — y and On{x) are local fields. This 
expansion is valid as long as a; — y is small enough or equivalently, is large enough 
compared to relevant mass scale. At any fixed values of x — y only finite number of 
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terms contribute in the above expansion. 

According to the OPE, at large » Ag^^, the moments of the structure 
functions can be expanded in powers of 1/Q^. The coefficients in the expansion 
are matrix elements of quark and gluon operators corresponding to a certain twist, 
r, defined as the mass dimension minus the spin, n, of the operator. For the n-th 
moment of the F2 structure function, M2"'\ for example (see Eq. (43)), one has the 
expansion 

Mt\Q')= E ^'"'^?if^\ n = 2,4,6... (1.41) 

r=2,4,... ^ 

where are the matrix elements with twist less than r. The dependence of the 
matrix elements can be calculated perturbatively, with expressed as a power series 
in Q!,(Q2). 



The moments of structure functions defined in Eq. 1.39 can be used in a QCD 
analysis. It must be emphasized that OPE analysis is perturbative since the structure 
functions are expanded in terms of a hard scale, l/Q"^- The usefulness of OPE is that it 
allows a separation of "soft" , nonperturbative physics contained in parton correlation 
functions, from the "hard" scattering of the probe from the partons |47j . 

The relation between the higher-twist matrix elements and duality in electron 
scattering was explained in the paper of De Rujula, Georgi and Politzer fl2[ H8] . 
They explained the empirical observation of Bloom-Gilman duality in terms of the 



twist expansion of the structure function moments in QCD. Examining Eq. 1.39 they 
saw that the lowest moment corresponding to n = 2 is the Bloom-Gilman integral 



defined in Eq. 1.37 Scaling (Q^ independence of the F2 structure function) is only 
observed at high enough values of Q^, which means that at low duality should be 
violated. At high Q^, where the scaling of the structure function is observed, moments 
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become independent of . In order for the integrals of the structure functions at 
higher momentum to be equal to the integral of the structure functions at low Q^, the 
higher-twist contributions must be small or cancel. This is equivalent to saying that 
duality can exist if the interactions between the scattered quark and the spectator 
system is suppressed. 
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1.5 The Nuclear Dependence of R 

Pions are thought to be the carrier of the nuclear force at nucleon distance scale 
and hence are the dominant agent responsible for the binding of nuclei. The Yukawa 
interaction can be used to describe the strong nuclear force between nucleons carried 
by pions. This was strengthened via the observation of very significant effects of pion 
exchange currents in a variety reactions involving electromagnetic probes of nuclei. 
In the EMC experiment [l9j it was discovered that the structure function F2 depends 
on the mass number A of the target. This effect has been studied by comparing 
F2 measured on bound nucleons in nucleus A and deuterium. This implies that 
there is a significant difference between the parton distributions of free nucleons and 
nucleons in a nucleus. There were several attempts to understand the observed effect 
in terms of nuclear pions but all of them failed to explain the experimental data at 
all kinematics. These failures to observe the influence of nuclear pions caused a crisis 
for nuclear theory |50j . 

One possibility for verification arises from measurements of the ratio R, of scatter- 
ing of virtual photons in a longitudinal or transverse polarization state , R = ai/ax 
from nuclei. A large value of ai, and the corresponding violation of the Callan-Gross 



relation Eq. 1.30, indicates the presence of nuclear bosons as fundamental constituents 
of nuclei [TT] . A measurement of R in the region covered by this experiment combined 
with the data of experiment E02-109 [18] can provide information about the dynamics 
of the internucleon force inside nuclei as predicted by Miller [21]. The calculations 
performed by Miller [51] predict a significant pion excess enhancement in the ctl/'^l 
ratio at low Q2 and moderate shown in Fig. |6j The phenomenology is de- 

scribed in the reviews [521 ESI SS] • The main point is that the quark and anti-quark 
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Figure 6: Plot from G. A. Miller [51] showing the predicted sensitivity of the inclusive 
longitudinal cross section ratio of iron to deuterium due to pion excess. 

distributions can be given as convolution of the g, q distributions in a given nuclear 
hadronic constituent with the light cone distribution functions: fn/A{y), fw/Niy), the 
probability to find an excess pion in the nucleus (A) or nucleon (N), with a plus- 
momentum given by yrriN- The function fT,/N{y) is constrained by the data on the 
nucleon sea which restrict the u and d distributions to be similar [311 ES] . 

The longitudinal cross sections for nuclear and deuterium targets are related by 
the formula [51] 

g-l+4/.K) (/_^^,)^g^ (l + fi.) (1.42) 

where is the pion form-factor. This equation is plotted shown in Fig. [6] for iron 
an target. The largest sensitivity to pions is predicted at = 0.3 GeV^ and quickly 
vanishes with the increase of the due to the falloff of the pion form-factor. The 
results of this experiment, combined with the results of the E02-109 [T8] experiment, 
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will allow the extraction of (xf and af (Q^ ~ 0.3 GeV^, < 4.0 GeV^). The 
expected enhancement in the af^ / erf ratio at low will be regarded as an indication 
of pion excess in heavy nuclei as suggested by Miller [5T] . 
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1.6 Neutrino Scattering 

Neutrino interactions with nuclei are important for the determination of the internal 
structure of the nuclear matter and complement deep inelastic electron scatterings ex- 
periments. Electron scattering experiments are sensitive to the average of the square 
of the electric charge of the nucleon constituents, while neutrino experiments measure 
their weak structure. In particular, neutrino(antineutrino)— nucleus scattering can be 
used to understand the flavor composition of nuclei as neutrino scattering is sensitive 
to specific flavors of quarks. 

It is also interesting to test if Bloom-Gilman duality is observed in neutrino scat- 
tering. It was argued in Ref. (56] that duality should also exist in weak structure 
functions. This can only be confirmed by experiment and there is not enough data 
in the resonance region to draw a final conclusion. The interaction of neutrinos with 
the nuclear matter is weak and cross sections are small which requires the use of 
nuclear targets. Also nuclear effects in DIS have been studied using only muon and 
electron beams. Neutrino scattering experiments were only done using heavy nuclear 
targets such as iron target-calorimeters. Results of these experiments indicate that 
the nuclear corrections for e — A and v — A are different. Among these differences is 
evidence for quark- flavor dependent nuclear effects [FT] . 

The MINERvA [57j experiment seeks to measure low energy neutrino interactions 
both in support of neutrino oscillation experiments and also to study the strong 
dynamics of the nucleon and how the nuclear medium affects these interactions. 
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1.6.1 Structure Functions in DIS 

Deep inelastic neutrino structure functions are directly related to the parton distri- 
bution functions and the same parton distribution functions determine the charged 
lepton scattering. For electron scattering the unpolarized parton distribution q{x) for 
both spin-1/2 and spin-0 partons is related to the nucleon structure functions accord- 



ing to Eq. 1.29 where the sum is over the various species of partons of charge e, and 



fi{x) is the probability that the parton has momentum in the interval x — )■ x + dx. 



Assuming that partons have quantum numbers of quarks and using the Eq. 1.29 F2 



electron scattering structure functions of proton and neutron can be written as: 



l^f = ^[uP{x) + uP{x)] + hdP{x) + dP{x)] + hs^ix) + sP{x)] (1.43) 
X t) 9 9 



-Fr = ^K(^) + «"(^)] + ^K(^) + d'^ix)] + hs^ix) + rs"ix)] (1.44) 
X i) y y 

where proton is indicated as p and neutron as n. Since u, d quarks and proton, 
neutron both form isospin doublets one can see that = c?" (call this u), d^ = c?" 
(call this d), = (call this s). Here the contribution of other quarks is small 
and is not considered. Taking into account the isospin symmetry described above the 
structure functions can be written as: 

l^f = ^(u + u) + hd + d + s + s) (1.45) 
l^f = l(ci + J) + l(u + ^ + s + s). (1.46) 
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For neutrino nucleon scattering F2 can also be written in terms of quark distribu- 
tions. Charged current neutrino scattering is different from the electron scattering 
since it only couples to specific flavors of quarks, in particular for neutrino proton 
scattering only (i, s, m, c quarks contribute, and in neutrino-neutron scattering only 
c, (i, s quarks contribute. Neutrino scattering structure function can be written as 

F2(x) = x^/,(x), (1.47) 

i 

where only specific quarks contribute as discussed above. At energies where the charm 
quark production is absent or negligible and the Cabibo angle [58j (The Cabibo angle 
is the probability that one flavor of quark will change into other flavors under the 
action of the weak force) is zero, the proton structure function can be written in terms 
of quark distributions as 

-Ff{x) = 2[d{x) + u{x)] (1.48) 



X 



-K^(x) = 2[u{x) + d{x)] 



;i.49) 



where the factor 2 comes from the fact that weak current contains both vector and 



axial parts. Using equations 1.45 1.46 , 1.48 , 1.49 the structure functions of electron 
and neutrino scattering can be related to each other as 



+ Ff 

puN ^ pyP 



2 



2 J 



[X 



'u + u + d + d + 2/^{s + s) 
u + u + d + d 



;i.5o) 
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If the contribution of strange quarks is negligible [x > 0.2) the structure functions are 
related by a factor of 5/18, known as the 5/18ths rule. The observation that electron 
scattering and neutrino scattering structure functions are approximately related by a 
factor of 5/18, was a significant triumph for the QPM. An interesting and very useful 
review about this QPM can be found in [2T] . 

1.6.2 Neutrino Scattering in Resonance Region 

In the previous section relating the structure functions of electron and neutrino scat- 
tering from the nucleon in DIS region was rather straightforward. The same can not 
be said in the nuclear resonance region. In the resonance region higher order QCD 
effects are not negligible as in DIS and it is not possible to consider nucleons as non- 
interacting group of partons. Since the electromagnetic and weak interactions are 
sensitive to different types of partons (quarks) the relationship between the structure 
functions in these reactions can only be established if specific conditions are satisfied. 
For example a neutrino beam can convert a neutron into a proton, but it cannot 
convert a proton into a neutron (and vice versa for an antineutrino beam). From the 
conservation of the vector current, the vector structure functions measured in electron 
scattering can be related to their counterparts in neutrino scattering for only specific 
isospin final states. The axial structure functions in neutrino scattering can only be 
determined if the vector part is provided from electron scattering. 

In order to give a general idea how the structure functions in electron scattering 
can be related to their counterparts in neutrino scattering, electron scattering cross 
section is written in terms of helicity amplitudes. 
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where T-r is the virtual photon flux, e^, describes the degree of longitudinal polarization 

1 /2 

of the photon, (Tj! is the total transverse absorption cross section with helicity 1/2 
for the photon-nucleon system, and a^^"^ is the helicity 3/2 cross section. Here ctl is 
the total cross section for the absorption of a longitudinal (scalar) photon. Analyses 
of electroproduction data give numerical values for cross sections at the peak of each 
resonance [591 EO] , 

ar{W = Mn) = {Aj/, + Al^,) (1.52) 



a,{W = Mn) = ^%Sl/, (1.53) 
TrMr qi ' 

where Tr is the width of the resonance, and Mr is the mass of the resonance. Using 
these formulas one can relate electromagnetic to weak form factors using isotopic 
symmetry. Photons can have two isospin states |/, J3 >= |1, > and |J, Is >= |0, >. 
The isovector component belongs to the same isomultiplet with the vector part of the 
weak current. Each of the amplitudes A3/2, ^1/2, 'S'1/2 can be further decomposed 
into three isospin amplitudes [61] . A general helicity amplitude on a proton {Ap) and 
neutron (A„) target has the decomposition 



Ap = Api^p ^R+) = b-^^a' + y^a^ (1.54) 



An = Aniin ^R^)=b+ ^J^a' + ^J'^a^ (1.55) 
For the weak current we have only an isovector component of the vector current. 
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therefore the b amphtude never occurs in weak interactions. The and can be 
found from weak interactions according to the following formulas 



A{W+n /?(^)+) = -^a\ (1.56) 



AiW+p = V2a\ (1.57) 



A{W+n i?(=^)+) = ^a^. (1.58) 



Comparing Eq. 1.56 with Eq 1.54 it can be seen that for the isospin 1/2 resonance 
the weak amplitude satisfies the equality A{W~^n R^^^+) = A\ — A^. Since the 
amplitudes are linear functions, the weak form-factors are related to electromagnetic 
form-factors with a simple relation 

1 = 1/2: = - Cf. (1.59) 

For the isospin 3/2 resonances the Al{W^n -> i?(^)+) = A'i{W~p -> i?(^)°) = 



y^2/3a'^. The weak form-factors are 



J = 3/2 : = = Cf. (1.60) 

The above description of relating the weak and electromagnetic structure functions 
in nuclear resonance region is thoroughly described in Ref. [61]. The article is self- 
contained and facilitates writing simple programs to reproduce the cross sections. 

The MINERvA experiment will perform a high statistics neutrino-nucleus 
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scattering experiment using a fine-grained detector specifically designed to measure 
low-energy neutrino-nucleus interactions accurately. The high-luminosity NuMl beam 
line at Fermilab will provide energies spanning the range ~ 1 — 15 GeV, over both the 
resonance and deep inelastic regimes, making MINERz/A a potentially very important 
facility to study quark-hadron duality in neutrino scattering. The results of this 
experiment will be used in future analysis of neutrino data (including MINERvA) in 
order to investigate quark hadron-duality in the axial structure functions of nucleons 
and nuclei. 



Chapter 2 
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Experimental Apparatus 



2.1 Overview 

The Jefferson Laboratory Hall C experiment E04-001 [22] was performed in May- June 
2007. The purpose of this experiment is to measure the longitudinal-transverse (L-T) 
separated structure functions F2 and R = a^/ ax from nuclear targets in the resonance 
region. The targets used include carbon-12, aluminum-27, iron-56 and copper-64. The 
data were taken at 2.1, 3.2, 3.3, 4.1, 5.1 GeV beam energies and electron scattering 
angle ranging from 12° to 76°. The scattered electron was detected in the High 
Momentum Spectrometer (HMS). The kinematic coverage of the experiment is shown 
in Fig. [7] and is listed in Appendix. 




2.097 GeV 
3.117 GeV 
E^= 3.267 GeV 
E^= 4.065 GeV 
Eb= 4.137 GeV 
Eb= 5.157 GeV 











2 



4 



(GeV^) 



Figure 7: Kinematic coverage of E04-001 experiment. 



2.2 The Thomas Jefferson Accelerator Facihty 
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Figure 8: Lay-out of the CEBAF facility. The electron beam is produced at the 
injector by illuminating a photocathode and the produced photoelectrons are then 
accelerated to 67 MeV. The beam is then further accelerated in each of two super- 
conducting linacs, through which it can be recirculated up to five times. The beam 
can be extracted simultaneously to each of the three experimental halls. 

The Continuous Electron Beam Accelerator (CEBAF) of JLab [63] was designed 
to deliver a continuous beam of electrons simultaneously to three experimental end- 
stations. A diagram of the racetrack shaped accelerator is given in Fig. |8} The source 
of the injector is a 100 kV photocathode guns [M] with a maximum beam current 
of a few mA. Next, the beam is incident on a chopping aperture which contains slits 
of different sizes, one for each of the experimental Halls A, B and C. The width of 
these slits determines the beam current that is delivered to each Hall. The chopper 
sweeps the beam over the slits with a rotating electric field with a frequency of 1497 
MHz. The beam then enters the first superconducting accelerator section, where it 
is accelerated to 67 MeV and then injected into the North Linac. The North Linac is 
a string of 20 cryomodules, with each cryomodule containing eight super-conducting 
niobium cavities. These cavities are kept super-conducting by 2 K Helium coolant 
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from the Central Helium Liquefier (CHL). Klystrons excite the cavities and establish 
electric fields which accelerate the electrons as they travel through them. At the end 
of the North Linac, the electron beam has a nominal energy of 600 MeV, although by 
careful tuning of the accelerating electric field of the cavities, this energy can be raised 
or lowered. Next, the beam enters the east magnetic recirculation arcs, where it is 
bent in a semi-circle to the South Linac. Here again the beam is accelerated through 
a string of 20 cryomodules. At the end of the South Linac, the beam can be extracted 
for use in any of the experimental halls or it can proceed around the west recirculation 
arcs for another pass around the accelerator. The beam can traverse the accelerator a 
maximum of five times, gaining a nominal 1200 MeV of energy with each pass around 
the machine. In the North and South Linacs the different energy beams, resulting 
from each pass around the accelerator, travel in the same beam-line. However, the 
different energy beams require different bending fields in the recirculation arcs. When 
the beams reach the arcs they are separated by momentum and each one goes through 
a different arc. At the end of the arcs, the beams are recombined into the same beam- 
line again. When the beam is of the energy requested by the experimental halls, it is 
extracted from the accelerator to the Beam Switch Yard (BSY). There the three sets 
of beam bunches are separated into the appropriate experimental hall beam-line by 
defiecting cavities operating at 499 MHz. Each hall receives a short (1.67 ps) train of 
pulses with a frequency of 499 MHz. The beam has a very small transverse size ( > 
200 //m (FWHM) at 845 MeV). The fractional energy spread (AE/E) is at the 10"^ 
level. 
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2.3 Beam-line Apparatus 

Electron beam parameters are measured by several devices (see Fig. [9]) situated along 
the beam-line, upstream of the target, near the entrance to the hall. Measurements 
of position, current, energy and longitudinal polarization are possible, with some 
intentional redundancy to allow cross checking. Some typical beam parameter values 
for the present experiment and the uncertainties associated with their measurement 
are given in Table [l] 



Beam Parameter 


Beam-line 


Measured 


Accuracy (absolute) 




Device (s) 


Value 




Position X (at target) 


BPM/ Superharp 




<200/im 


Position y (at target) 


BPM/Superharp 




<200/im 


Current 


BCM 


30-80 /xA 


<3xlO-VA 


Energy 


ARC 




2x10-^ GeV 



Table 1: Typical values measured for the various electron beam properties described 
in the text, together with the associated accuracies. 

To measure the absolute value of the energy of the electrons in the beam, eight 
dipole magnets between the beam switch-yard and the hall entrance deflect the elec- 
trons through a nominal angle of 34.3°. Wire scanners before and after the magnets 
accurately determine this angle, which along with the measured arc field integral 
J Bdl (magnetic field integral over the path of the electron beam), is used to calcu- 
late the electron momentum. 

2.3.1 Beam Position Monitors 

So-called superharps make a destructive but extremely precise measurement of the 
beam position and profile. A superharp consists of three tungsten wires, two vertical 
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HallC 



2 Siiperharps, 
1 BPM 



2 Superharps, 
1 BPM 



2 Superharps, 

3C Arc 1 BPM 




Figure 9: Schematic lay-out of Hall C, indicating the location of the the raster, the 
beam energy measurement system, the beam current monitors (BCM) and the beam 
position monitors (BPM) upstream of the target. 

that measure the horizontal beam profile and one horizontal wire that measures the 
vertical beam profile, mounted in a frame which is connected to an arm that can be 
moved in and out of the beam. Analog-to-Digital converters (ADCs) connected to 
each wire read the signals on the wires as the frame is moved in and out of the beam, 
while a position encoder determines where the wire intercepts the beam. With the 
position information and the ADC measurements the position and profile of the beam 
can be measured. More detailed information about the superharps can be found in 
Ref. [eg. 

The beam position is also measured with three Beam Position Monitors (BPMs) 
which provide non-destructive beam position information along the beam line which 



are used during data taking; see Fig. 10 A BPM is a cavity containing two pairs 
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of antennae perpendicular to each other and inchned by ±45° with respect to the 





Y 

,'' beam 


V'''' 


X 











Figure 10: Diagram of the orientation of the BPM antennae. 



horizontal plane. When the beam passes through the cavity parallel to the symmetry 
axis both pairs of antennae pick up the frequency of the beam. Each antennae gives 
a signal proportional to the distance that the antennae is positioned from the beam. 
The beam position is derived from the difference over the sum of the signals from an- 
tennae on opposite sides of the beam. This method provides a position measurement 
independent of beam current. Since this is a position relative to the central axis of the 
beam line these measurements must be compared to the Superharp measurements to 
determine the absolute position of the beam during data taking. The final accuracy 
of the beam position measurements is ±1.0 mm, with a relative position uncertainty 
of 0.1-0.2 mm. More detailed information about the BPMs can be found in Ref. 



2.3.2 Beam Current Monitors 



The Hall C beam line has four devices installed to measure the beam current: three 
microwave cavity Beam Current Monitors (BCMs), and a parametric DC current 
transformer (Unser monitor). 

The BCMs measure the integrated beam current in two second intervals, and their 
output voltage, as well as the output voltage of the Unser, is converted into frequency 
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with a voltage-to-frequency converter, and read out by a scaler. A BCM consists of 
a cylindrical wave guide, mounted in the beam line so that the beam travels along 
the axis of the cylinder. The dimensions of the cylinder were selected so that the 499 
MHz structure of the beam excites the 1497 MHz TMqio mode in the wave guide. 
The resonance frequency is picked up by wire loop antennae, and converted to a DC 
voltage through a RMS-to-DC converter. The cavities have a stable gain and offset, 
and a high signal/noise ratio, but cannot measure the absolute current, because the 
output power as a function of the measured beam current depends on the cavity 
impedance, quality factor, and the signal cable attenuation. The absolute calibration 
is done with the Unser monitor. 

The Unser monitor is installed between BCMl and BCM2. It cannot be used for 
charge measurements because it is sensitive to thermal fluctuations, resulting in large 
drifts in its zero-offset. However, because the gain is stable and well measured, the 
Unser monitor is used to calibrate the gain of the BCMs. More information about 
the Unser monitor can be found in Refs. |67t [68] . 

2.3.3 Beam Raster 

The electron beam generated at CEBAF is a high current beam with a small trans- 
verse size (< 200 fim FWHM). It can deposit a large amount of power in the target. 
In order to prevent local boiling in the cryotargets or melting of the solid targets the 
beam is rastered before striking the target. For that purpose the beamline is equipped 
with a pair of fast raster magnets, located 25 meters upstream of the target. The first 
set rasters the beam vertically and the second horizontally. The current driving the 
magnets is varied sinusoidally, at 17 kHz in the vertical direction, and 24.2 kHz in 
the horizontal direction. The frequencies are chosen to be different so that the beam 
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motion does not form a stable (Lissajous) figure at the target, but it moves over a 
square area. The amphtude of the raster pattern on the target was ±1.0 mm in both 
directions during the present experiment. More details about fast raster can be found 
in Ref. [05]. 



2.4 Experimental Hall C 

Hall C is the second largest of the three experimental halls at Jefferson Lab. It 
contains two spectrometers the High Momentum Spectrometers (HMS) and the Short 
Orbit Spectrometer (SOS) located on either side of the electron beam-line. The 
spectrometers can be moved clockwise or counter-clockwise about the target position 
over a wide range of angles. The HMS has a minimum angle of 10.5° with respect 
to the beam-line and a maximum central angle of 165°. The SOS has now been 
dismantled. Between the beam switch-yard and the target (which sits at the center 
of the hall) there are various devices to monitor and measure beam current, beam 
position, beam energy and beam polarization. A schematic of Hall C is shown in 



Fig. 11, At the pivot point of the two spectrometers sits the cryotarget, encased 
within a cylindrical aluminum scattering chamber. The electron beam is incident 
on the target through the beam— line. The targets used for this experiment were 
deuterium, carbon, aluminum, iron and copper. Scattered particles were detected by 
HMS detector. The primary beam continues along the path of the beam-line to the 
beam dump. 
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SOS 



HMS 



Beam line 




Figure 11: Schematic layout of Hall C. 



2.5 Target 



The targets are located in a cylindrical aluminum scattering chamber, which is in- 
stalled at the spectrometer pivot. The scattering chamber has an inner radius of 
61.6 cm and a height of 150 cm. The beam exit windows are made of 0.4 mm and 
0.2 mm thick aluminum foils on the sides facing the HMS and SOS spectrometers, 



respectively. The target ladder (see Fig. 12 ) is located inside the scattering chamber. 
The solid targets are BeO, carbon, copper, iron targets (See Table |2]). The targets 
are mounted on the target ladder and can be exchanged in few minutes by the lifter 
mechanism. The mechanism permits accurate, reproducible positioning of any of the 
targets at beam height. 



The cryotarget consists of three loops for circulating cryogenic liquids, see Fig. 13 
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Solid targets 



Dummy target 



Optic targets 




Liquid targets 



Beam direction 



Figure 12: Hall C target ladder. 



This experiment and E06-009 |69] experiment ran simultaneously and used the 



same target ladder shown in Fig. 13 The E06-009 used loop 1 filled with deuterium, 
while the E05-017 experiment, which also ran during the same time period, used loop 1 
filled with liquid hydrogen. The liquid hydrogen (LH2) was cooled down to 19.0 K and 
held at a density of 0.0723 ± 0.0004 g/cm^. The liquid deuterium (LD2) was cooled 



Target 


Purity 


Thickness (g/cm^) 


Radiation Length (%) 


BeO 


99.00% 


0.2918 ± 0.00030 


0.45 


Carbon 


99.95% 


0.35525 ± 0.00030 


1.06 


Copper 


99.95% 


0.17775 ± 0.00015 


1.81 


Iron 


99.95% 


0.11870 ± 0.00014 


1.00 



Table 2: Thicknesses and radiation lengths of the solid targets. 
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Figure 13: Hall C target ladder. 



to a temperature of 22.0 K and a density of 0.1670 ± 0.001 g/cm^. The temperature 
and pressure of the liquid targets were monitored by target control software which 
was programmed to sound alarms when the temperature or pressure in the targets 
changed beyond the predetermined limits. 

The liquid target cells are made of 0.0127 cm aluminum. The 4 cm dummy target 
(A16061-T6) consists of two dummy endcaps to simulate an empty target for target 
wall background measurements. The thicknesses and radiation lengths of the targets 
in the cryotarget ladder are listed in Table |4} 



54 



Foil Number 


Position 


Thickness (g/cm^) 


Radiation Length (%) 


Foil 1 


Upstream foil 


0.2658 ± 0.0035 


1.10 


Foil 2 


Downstream foil 


0.2549 ± 0.0034 


1.06 



Table 3: Dummy target thickness and positions. The alloy is A16061-T6. 



Target 


Thickness (g/cm^) 


Radiation Length (%) 


4 cm LD2 


0.6570 ± 0.0039 


0.56 


Cell Walls (27A1) 


0.0340 ± 0.0035 


0.14 


4 cm Dummy (^^Al) 


0.5210 ± 0.0017 


2.11 



Table 4: Thicknesses and radiation lengths of the targets in the cryotarget ladder. 
Cell walls represent entrance and exit foils. 

2.6 High Momentum Spectrometer 

The High Momentum Spectrometer (HMS) was one of the two standard spectrometers 
in Hall C. It was designed to have moderately large acceptance, good position and 
angular resolution in the scattering plane, an extended target acceptance, and a large 
angular range. During the current experiment the spectrometer was situated on the 
right-hand side of the beam-line (with respect to incoming beam direction) and was 
used to detect the scattered electron. 



2.6.1 HMS Optics Design 

The magnetic elements of the HMS consist of three superconducting quadrupoles and 



a superconducting dipole, in a QQQD arrangement (Fig. 14). The dipole magnet is 
a superconducting, cryostable magnet. Its basic parameters are an effective length of 
5.26 m, a bend radius of 12.06 m and a gap width of 42 cm. It was designed to achieve 
a 25° bending angle for 7.4 GeV/c momentum particles. It provides defocusing in 
the vertical plane to achieve good momentum resolution. Quadrupole Q2 focuses in 
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Detector Hut 




Figure 14: Schematic lay-out of a HMS magnetic elements, showing the geometrical 
configuration of the three quadrupole and the dipole magnets. 

the horizontal plane, whereas Qi and Q3 both provide vertical focusing. Horizontal 
focusing of Q2 provides a large momentum bite, solid angle and extended target 
acceptance. The HMS characteristics are given in Table [5) 



2.7 HMS Detector Package 

The detector package of the HMS consists of two drift chambers, DCl and DC2, two 
pairs of scintillator hodoscopes, SI and S2, a gas Cerenkov detector and a lead-glass 
calorimeter. The detectors are mounted on frames that connect to the carriage that 
supports the magnets. This design insures that the detector package and magnets stay 
unmoved relative to each other. The detector package has the following functions: 
triggering, tracking and particle identification. A schematic view of the HMS detector 



package is shown in Fig. 15 
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Maximum central momentum 


7.4 GeV/c 


Momentum acceptance 


±10% 


Momentum resolution 


<0.1% 


Solid angle 


6.7 msr 


Scattering angle acceptance 


±40.0 mr 


Out-of-plane angle acceptance 


±80.0 mr 


Extended target acceptance 


10.0 cm 


In plane angle resolution 


0.4 mr 


Out of plane angle resolution 


0.9 mr 


Useful target length 


10.0 cm 


Vertex Reconstruction Accuracy 


2.0 mm 



HMS Vacuum pipe exit 
DC1 



Table 5: HMS performance characteristics. 



Cerenl(Ov 



Calorimeter 



DC2 



S1X S1Y 



S2X S2Y 




Figure 15: Schematic side view of the HMS detector package. 



2.7.1 Drift Chambers 

A drift chamber is a particle tracking detector that measures the drift time of ioniza- 
tion particles in a gas to calculate the hit coordinates of ionizing particle. Combined 
with knowledge of the optical transfer properties of the spectrometer, the drift cham- 
ber hit coordinates are used to reconstruct the particle trajectory, reaction vertex and 
momentum at the target. 

The drift chambers are spaced 81.2 cm apart, and each has an active area of about 



57 



Amplifier/Discriminator 
Y Y' Cards 




Incident particle direction ("Z") 
(into the page) 

Front View 



] 

] ^ 

] 

] 

] 

] 



X, X' 



X Y U V Y' X' 



Incident particle direction 



Side View 



Figure 16: Drift chamber 

113 cm (x) by 52 cm (y). They consist of six separate planes of sense wires (anodes) 
of 25 fim diameter gold-plated tungsten spaced 1.0 cm apart in a gas mixture of argon 
and ethane, and field wires (cathodes) of 150 /xm gold-plated copper beryllium wires. 
The planes are spaced 1.8 cm apart. Between each sense wire is a field wire which is 
held at a negative potential (-1800 V to -2500 V). On both sides of the sense planes is 
an additional plane of field wires, held at the same negative potential. The planes are 
ordered X, Y, U, V, Y , X as seen by incoming particles. The X and X planes provide 
two measurements of position of the particles in the dispersive direction. The Y and 
Y' planes offer two measurements of position in the transverse direction and the U 
and V are rotated ±15° from the X (X ) planes to prevent the right-left ambiguity. 



see Fig. 16 



When a charged particle passes through one of the chambers, it ionizes the gas 
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atoms and produces a trail of electrons and ions. Far from a sense wire, where 
the electric field between wire and cathode is uniform and parallel, electrons drift 
toward the wire with a constant velocity. As they get closer, the field has a stronger 
radial nature, causing the electrons to accelerate and produce a secondary electron 
avalanche. This avalanche generates a negative pulse on the wire. The signals from 
each wire are amplified and discriminated on the cards attached directly to the drift 
chambers and then sent to the TDCs (Time-to-Digital Converters) located in the 
back of the detector hut. The measured drift time and known electron drift velocity 
can then be used to calculate the perpendicular distance between wire and particle 
track. A typical track generates signals in about five wires per plane. The position 
resolution at the focal plane is approximately 280 fim per plane. More information 
about the HMS drift chambers can be found in Ref. [70] . 

2.7.2 Hodoscopes 

The HMS has 4 hodoscopes, which provide the trigger for detector read-out and 
allows the identification of heavy particles through time-of-fiight (TOF), though this 
was not used in the present experiment. The hodoscopes are paired in two horizontal- 
vertical X-Y sets (SIX, SlY and S2X, S2Y). The sets are separated by 220 cm. Each 
X hodoscope consists of 16 horizontally oriented scintillators (paddles), while the Y 
hodoscopes consist of 10 vertically oriented scintillators. The paddles of the X and Y 
hodoscopes are all 1 cm thick and 8 cm wide, but they have different lengths. The 
X paddles are 75.5 cm long, and the Y paddles are 120.5 cm long. Each scintillator 
paddle is wrapped in one layer of light-tight aluminum foil and two layers of Tedlar. 
At both ends of each paddle are Photomultiplier Tubes (PMTs) attached to the 
paddle through Incite light guides. The paddles are staggered in the beam direction 
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with 0.5 cm overlap between the paddles to avoid gaps, see Fig. [T7| When charged 
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Figure 17: Hodoscope geometry. 



particles pass through the paddles they excite the atoms of the scintillators. These 
atoms emit light as they return to their ground state. The light is detected by PMTs 
at the ends of the paddles. The light that is not emitted along the length of the 
paddle is reflected internally through the scintillator and ultimately also detected by 
the PMTs. 

The signal pulses from the PMTs are sent to the counting house where they 
run through the splitter, giving two signals with 1/3 and 2/3 of the amplitude of the 
original input signal. The smaller signal goes to Analog-to-Digital Converters (ADCs) 
that measure the integral of the signal. The other part of the signal is discriminated 
and one set of outputs is sent to TDCs (for timing information). The other set of 
outputs is sent into a logic module. The logic module first generates the logical OR of 
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all the discriminated signals from the tubes on one side of a given plane, for example: 
S1X+ = (S1X1+ OR S1X2+ OR ... S1X16+). There are equivalent sets of signals for 
the — side of each plane. Then, these sets of signals are combined into six outputs: 

a) SIX = (S1X+ AND SIX-) and analogously for SlY, S2X, S2Y. These four 
output logic signals indicate which of the hodoscope planes are active and make a 
new logic signal in case at least three of them have fired which is called SCIN. 

b) The X— Y pairs are further combined to form SI = (SIX OR SlY), and S2 = 
(S2X OR S2Y). These signals indicate whether the pairs are active and make a logic 
signal called STOF. 

These logical outputs are then sent to the main trigger logic and to the scalers to 
be recorded. A more detailed description of the hodoscopes can be found in Ref . [71] . 

2.7.3 Gas Cerenkov Detector 

The HMS gas Cerenkov detector provides particle identification by operating as a 
threshold detector. It consists of a large cylindrical tank (diameter of 150 cm, length 
of 165 cm) situated in the middle of the detector stack between the hodoscope pairs 



SI and S2, see Fig. 15 A pair of front reflecting spherical mirrors mounted vertically 



with 1 cm overlap at the rear of the detector is rotated over 15 degrees to focus 



the light on a pair of PMTs, see Fig. 18 The detector is filled with C4F10 gas. The 
Cerenkov detector measures the light emitted when a charged particle travels through 
the gas with a velocity above the speed of light in the gas. This is known as Cerenkov 
radiation. The light will be emitted with an angle cos{6) = l/n/3, where /3 is the 
velocity of the particle relative to the speed of light and n is the index of refraction of 
the material. If n/3 < 1 no light will be emitted. The light is reflected from focusing 
mirrors to PMTs, which generate a signal proportional to the number of Cerenkov 
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photons. The gas (and here the index of refraction) is chosen such that electrons at 
the spectrometer momentum will emit Cerenkov radiation and pions will not. The 
refractive index of the C4F1Q gas at 1 atm is n = 1.0006 and this gives a pion threshold 
above 4 GeV/c and an electron threshold of about 15 MeV/c. The average measured 



signal from an electron is about 10 photoelectrons, see Fig. [T9j However, it is still 
possible for a pion to be misidentified as an electron when it produces a knock-on 
(5-electron that fires the Cerenkov detector. 

The signal from each PMT is sent to the counting house where each signal is split 
in a fashion similar to the hodoscope signals. One pair of signals is sent to the ADC 
and the other pair is summed and put through the discriminator to give signals for 
the TDC and trigger logic. More information about the Cerenkov detector can be 
found in Ref. [71]. 
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Figure 19: HMS Cerenkov spectrum. Most of the pions appear at zero photoelectrons. 



2.7.4 Electromagnetic Calorimeter 

The HMS has a lead-glass calorimeter which is used to discriminate between electrons 



and pions. Fig. |20] depicts the calorimeter which consists of 52 TFl type lead glass 
blocks 10 cm x 10 cm x 70 cm, with a PMT on one end. The blocks are arranged 
in four layers with 13 blocks per layer, giving a total thickness of the calorimeter 



along the direction of particle motion of 16 radiation lengths. As shown in Fig. 15 the 
calorimeter is rotated 5° with respect to the dispersive plane to prevent particles from 
passing between the blocks. Electrons interacting with the lead glass radiate photons 
in the calorimeter, which in turn produce electron-positron pairs (when the photons 
are energetic enough). These pairs in turn also radiate photons and so a shower of 
particles (e"'",e~,7) is produced in the calorimeter. The charged particles produce 
Cerenkov radiation which is detected by photomultiplier tubes. The produced signal 
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Figure 20: Calorimeter geometry. 



is proportional to the total track length of the particles in the calorimeter which is in 
turn proportional to the energy of the initial electron. Electrons (positrons) entering 
the calorimeter deposit their entire energy, and thus the ratio of deposited energy 
of electrons (positrons) in the calorimeter to the detected energy derived from the 



particle bending in the spectrometer is unity, see Fig. 21 



Pions normally deposit about 300 MeV through ionization in the calorimeter so a 
peak in the Ecaio/E' spectrum can be observed at 0.3 GeV/E'. However, pions can 
have a charge-exchange reaction and produce a neutral pion, which in turn decays 
into two photons, the full energy of which will be deposited in the calorimeter. This 
leads to a high-energy tail for pions, which can result in pion misidentification. 

The signals from the calorimeter PMTs are sent to the counting house where they 
are split 50/50. One half is sent to the ADC and the other half to the linear modules 
to be summed. The sum in the first layer (PRSUM) and the sum of the entire 
calorimeter (SHSUM) are discriminated to give three logic signals for the trigger. 
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Ecalo/E' 



Figure 21: Calorimeter energy distribution divided HMS momentum 
(Phms=1.0 GeV/c). 

PRSUM is the sum of all signals from the first layer of the calorimeter. The SHSUM 
signal, obtained by summing the signals from all lead-glass blocks, represents the 
total energy deposited in the calorimeter. PRSUM are used to form the high (PRHI) 
and low (PRLO) thresholds on the energy in the the first layer of calorimeter. SHLO 
is a cut on the total energy in the calorimeter. 

2.8 Data Acquisition System 

The standard Hall C DAQ system consists of a variety of CAMAC, NIM, Fastbus 
and VME electronics, as well as modules custom built by the JLab electronics group. 

The data acquisition is handled by the CODA (CEBAF Online Data Acquisi- 
tion) [72j software package running on a PC. Data for each run are written directly 
to rotating memory and consists of three types of events: 1) detector information 
handled by the ADCs and TDCs, 2) scaler information, and 3) information from the 
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EPICS [73] database. The ADCs and TDCs are read-out for each event, while the 
scalers are read every 2 seconds. The EPICS database contains information such as 
magnet settings, beam position and target temperature and pressure. These quanti- 
ties are read out every 30 seconds. More information about data acquisition can be 
found in Ref. [72]. 



2.8.1 Triggers 



Triggers generated by the PMT signals from the scintillator planes and from calorime- 
ter blocks are sent to DAQ. 

The physical information of interest is ultimately obtained also from raw data 
(digital numbers) read out from hundreds of ADC and TDC channels. 
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Figure 22: Schematic diagram of the HMS trigger logic. 
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Charged particles passing through the spectrometer produce triggers in one or 
more of the detectors described earher. Certain combinations of these triggers are 
used to form the pretrigger. An electron trigger (ELREAL) can be produced in two 
ways: 

1) The low — level electron trigger (ELLO) requires a Cerenkov signal (C), plus at 
least two out of three of the following conditions: a) at least one of the two scintillator 
layers of each hodoscopes has fired (STOP = SI AND S2). b) at least three of the 
four scintillator layers of both hodoscopes have fired (SCIN). c) There is a (PRLO) 
signal from calorimeter. 

2) A high — level electron trigger (ELHI) requires that all of the following signals 
are present: a) The (SCIN) signal, b) The (PRHI) signal from the calorimeter, c) 
The (SHLO) signal from the calorimeter. The high — level electron trigger (ELHI) 
does not use the Cerenkov signal. 

The electron trigger (ELREAL) is (ELLO) OR (ELHI). The reason for using 
two electron triggers (ELLO) and (ELHI) is to reduce the trigger inefficiency and 
to provide the most efficient electron-hadron separation by the Cerenkov and the 
calorimeter detectors. 

There is an additional pion trigger (PION) which requires the (SCIN) signal, and 
no Cerenkov signal (C). The (PION) signal is then combined with a prescaling circuit 
to form (PIPRE), which ensures that a low-rate sample of pions is sent along to the 
data acquisition system. Two copies of (ELREAL) are generated: one is fanned out 
into four logic units with dead times set between 30 and 120 ns in order to allow a 
measurement of the electronic dead time and another (ELREAL) OR (PIPRE) forms 
the signal (PRETRIG), which is forwarded to the Trigger Supervisor (TS) 
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When a run is started first 1000 pedestal triggers are generated by (PED PRE- 
TRIG) and the data acquisition system records the read-outs from the ADCs. Then 
the data acquisition system begins to record physics events. When a physics event 
is signaled by (PRETRIG), the data acquisition system records the read-outs of the 
detectors ADCs and TDCs. The TS sets TS BUSY while the data acquisition system 
records the read-outs. After the data acquisition system has completed recording 
the event information, TS BUSY is set off to allow for the next physics event. Both 
trigger and pretrigger signals are written into scalers providing information about the 
computer dead time. 
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Chapter 3 



Data Analysis 



The goal of data analysis is to perform and apply detector calibrations and extract 
cross sections from experimental data. In this chapter the calibration of detec- 
tors, along with background subtraction and cross section extraction methods are 
described. 

3.1 BCM Calibrations 



The BCMs measure the integrated beam current in two second intervals, an example 



is shown in Fig. [23| The Unser monitor is used to calibrate the BCMs since it has 
very well measured and stable gain. Despite its stable gain, the Unser's offset is very 
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Figure 23: Unser and BCM frequency versus time. 



unstable and can change depending on temperature and other physical conditions. 
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This is the reason it is not used for an absolute beam current measurement. The 
BCMs have very stable offsets but their gains can change over time and need to be 
calibrated a few times during this experiment. 

During a BCM calibration run the accelerator delivers beam to the hall in current 
steps lasting 2 minutes, first increasing up to the current at which data are taken, then 
decreasing, with equal time intervals of no beam between the steps. During the beam 
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Figure 24: BCM calibration results for the entire period of the experiment. The 
time interval between first and last measurement is two months. This shows that 
BCM gains and offsets were very stable. 



off intervals the zero offsets of the Unser monitor are measured by calculating the 
average frequency of the Unser before beam was on and after beam was off (before 
2 minutes interval and after it). This allows a very accurate determination of the 

Unser monitor's offset during the 2 minutes beam on interval. The beam current for 
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each beam on period is calculated using the measured Unser offset and gain. Plotting 
the BCMs average frequency, during the beam on period, versus the Unser current 
and doing a linear fit one can extract the BCMs' gains and offsets. These gains and 
offsets are used in the replay engine to calculate the beam current. 
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Figure 25: BCM cahbration residuals. 



In the E04-001 experiment the beam current for the iron target was limited to 
maximum of 40 //A, while for carbon, aluminum and copper up to 80 //A. For produc- 
tion runs the beam current was always higher than 35 /xA to minimize the uncertainty 
in the charge. 

The systematic error from the cahbration procedure, due to noise of the Unser 
monitor signal and uncertainty of the gain of the power meter signal used to mea- 
sure the beam current signal from the BCM, yields on overall absolute systematic 
uncertainty on the current measurement of 0.3 /iA. This was estimated by analyzing 
carbon runs taken with same kinematics but with beam currents ranging from 10 //A 
to 100 iiK. All runs were corrected by a factor 1.0/(1 -|- 0.3//) to take account this 



71 



offset. 

3.2 Detector Calibrations 

For all previous Hall C experiments, the software package HALL C Engine has been 
used for data analysis. It is a code based on GNU FORTRAN 77, which utihzes 
many of the data manipulation and display tools that are part of the widely-used 
CERNLIB package. It involves techniques for the HMS data analysis which, as a 
result of copious use in recent years, are stable and well-tested. The principal element 
is drift chamber (DC) tracking and the closely related magnetic reconstruction of the 
scattered electron momentum, direction and its reaction vertex at the target. In this 
chapter calibration of the DC, the Cerenkov, and the Calorimeter will be discussed. 

3.2.1 Drift Chamber Calibration 

Drift chambers provide the tracking information, coordinates and angles, for particles 
entering the HMS. The time difference between the fast START signal from the 
hodoscope counter and the STOP signal from the drift chamber is used to calculate 
the track position. In order to avoid negative drift times the overall offset between 
the times measured by the drift chamber and the times measured by the hodoscope 
is removed by doing hodoscope timing calibration. The time taken for electrons 
produced by ionization along a electron track to reach the anode wire is measured 
by TDC's and converted into a perpendicular distance from the wire. Since the hit 
time in a particular TDC depends on specific cable lengths and signal processing 
times, that may differ from channel to channel, a reference time (to) is found for each 
wire. This is used to extract the drift time from the actual TDC time. In order to 
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Figure 26: Drift time and distance distributions for the HMS drift chamber plane 
XI. These distributions are the sum of all production runs. 

determine how far the track was from the wire we generated a time-to-distance map 
using the following procedure. About 200000 events are analyzed to obtain the drift 
time spectrum. Assuming that hits are distributed evenly around the sense wire one 
can obtain the drift time by the following formula: 



J F{r)d7 
D(t) = Do^f^ 



(3.61) 



F{T)dT 



where D is the calculated distance from the wire, Do is the drift cell size (5 mm), T 
is the TDC time, tmim tmax are time limits corresponding START to STOP signals 
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(250 ns), F{t) is the measured drift time distribution. 

At high rates there are many accidental events which can force the tracking algo- 
rithm to fail and therefore introduce a tracking inefficiency. In order to reduce the 
number of these accidental events a narrow HMS drift chamber TDC time window 
from 2400 to 2900 channels (before it was 1800-3300 channels, 1 channel 0.5 ns) was 
used. 

3.2.2 Cerenkov Calibration 



The Cerenkov calibration is accomplished by finding the minimum ADC signal cor- 
responding to the production of one photoelectron in the PMT. To see the one pho- 
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Figure 27: ADC spectrum of the Cerenkov detector, PMTl and PMT2. 

toelectron peak one should select events that originate from less populated regions of 
the HMS acceptance to allow a minimal amount of light into the PMT. By doing so, 
it is possible to see the one, the two and even the three photoelectron peaks. The 



results of the calibration are shown in Fig. 27 where three different peaks are clearly 
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visible. For PMTl the one photoelectron peak is located at ADC channel 150, the 
two at 300 and a third peak three at 448 indicating that the calibration is done cor- 
rectly (the pedestal peak is centered at the zero ADC channel), since two and three 
photoelectron peaks are multiples of the one photoelectron peak. The same is true 



for PMT2. The ADC spectrum shown in Fig. |27| is the sum of all production runs. 
During the entire experiment the PMT pedestals had widths about 30—40 channels 
and were very stable. The calibrated signals from two PMTs are summed to get the 
total number of photoelectrons produced for each event. 



As described in Section 12.7.31 the Cerenkov detector has two mirrors. Due to 
engineering difficulties there is a gap at the region where the mirrors meet. The gap 
causes photon detection inefficiency around the central part of the 5p acceptance of 
the HMS. In order to estimate the size of the inefficiency and correct the cross sections 
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Figure 28: Inefficiency of the central part of the Cerenkov due to a gap between 
the mirrors. Inefficiency is parametrized in 5p and applied to all cross sections on 
bin-by-bin basis. All production runs are used to calculate the inefficiency. 
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the following procedure is done: the ELHI (see Section 3.4.7) trigger is required to 
fire in addition to all cross section cuts, the experimental yield is calculated with a 
Cerenkov cut with the number of photoelectrons greater than and greater than 2 
(see Sec. 



3.4.2 



for Cerenkov cuts) and the ratio of the yields is plotted versus 5p. The 



ratio is scaled by a constant to remove the effect of the Cerenkov cut 2 on pions (this 



exposes the inefficiency only due to the gap). As it can be seen from the Fig. 28 the 
Cerenkov cut 2 causes up to 2% inefficiency at the central part of the 5p spectrum. 
The ratio is parametrized in 5p and is applied to the cross sections (for \5p\ < 2% 
range) on bin-by-bin basis using the following formula 



a{W^,eUrr = (T{W^,0)/ / f{Sp)d6p. (3.62) 



It is important to mention that except for the inefficiency caused by the gap no 
position dependent ( 5p ) inefficiency is found. The point-to-point uncertainty (fit 
residuals) due to this parametrization is estimated to be 0.15%. 

3.2.3 Calorimeter Calibration 

Shower detector counters measure the energy deposited by the incoming particle. 
While passing through the lead glass of the shower counter, photons and electrons 
produce secondary photons and electrons, thus leading to an electromagnetic shower. 
The secondary electrons emit Cerenkov light, which after being reflected from the 
inner surface of the shower blocks, is collected on the cathodes of the photomultiplier 
tubes(PMT), mounted at the end of each block. The registered light is linearly pro- 
portional to the energy deposited by the incoming particle. Light flashes on the PMT 
cathode produce electrical signals, which are sent to ADC's and are stored in the 



76 




Figure 29: Calorimeter resolution. The top plot shows calorimeter resolution for 
all production runs versus HMS momentum. In order to find the energy resolution 
of the calorimeter the distribution of ratio of calorimeter track energy to the HMS 
momentum is fitted by a gaussian function. The width of the gaussian function gives 
the calorimeter resolution. The bottom plot shows the center of the gaussian function 
versus run number. 
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data stream for further analysis. The purpose of the shower cluster reconstruction in 
a shower detector is to calculate energy and in some cases the position of the particle. 
For this, at first shower cluster is identified. The cluster in a shower detector is deter- 
mined as a group of adjacent shower counters, where an electromagnetic shower has 
developed, i.e. an energy deposition was detected. The block that has the maximum 
energy deposition is called the central counter of the cluster. 

The parameters of the shower — the energy and the coordinates — are calculated. 
The energy E is calculated as a sum of energy deposition in all of counters included 
into the cluster by the formula 



i — number of shower counter, included into the cluster; 
M — set of counters numbers, included into the cluster; 
Ei — energy deposition in the i-th counter; 

For the shower parameters reconstruction, the shower counters gain calibration 
is required. The first step involves using the HMS electron energy Eg. 

The purpose of calibration is to define a coefficient for transformation of the ADC 
amplitude to the energy deposition for every shower counter, i.e., define coefficients 

Cj such that. 




(3.63) 



where: 



Ei — Ci ■ {Ai — Pj) , 



(3.64) 
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where: 



Ai — ADC amplitude; 

Pi — pedestal of the amplitude; 

Ci — calibration coefficient; 

Ei — energy deposition in the i-th counter. 
The calibration coefficients Ci are calculated by minimization of the functional 



n = 1 -T- — Number of event; 

i — Number of blocks, included in the cluster; 

— Set of counters numbers in the cluster; 
A" — Amplitude in the i-th counter; 

Pi — Pedestal of the i-th channel; 

— Known energy of a particle; 

Ci — Shower counters calibration coefficients to be fitted. 




(3.65) 



n=l ieM" 



Here: 



After obtaining the calibration coefficients, the energy deposition E is calculated 
by the formula (Eq. 3.63). 
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3.3 Target Coordinate Reconstruction 

In the transport coordinate system the trajectory of the particle at the focal plane, 
at the target and through the HMS magnetic elements is described by a vector (t), 
which expresses the track relative to the central reference trajectory. This vector is 
characterized by five components: 



(3.66) 



where: 

• X is the dispersive (vertical) displacement from the central trajectory expressed 
in meters. 

• x' is the angle the trajectory makes in the dispersive (vertical) plane relative to 
the central trajectory (dx/dz), expressed in radians. 

• y is the displacement in the non-dispersive (horizontal) plane. 

• y' is the angle the trajectory makes in the non-dispersive (horizontal) plane 
relative to the central trajectory (dy/dz), expressed in radians. 

• 5 is the fractional deviation of the momentum from the central value {p — Po)/Po, 
expressed in percent. 



X 

x' 

y 
y' 

5 



80 

In a first-order approximation, a transport matrix can be defined to relate the mea- 
sured focal plane coordinates to their counterparts at the target: 
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(3.67) 



This matrix involves only four unknown parameters since Xtar is known from the 
BPM data. In practice, the expansion of the focal plane coordinates is performed up 
to the fifth order. The reconstruction is performed with this formula: 

<ar = E"^, , Mi,^^(^f,yiyf,)\x'J{y'f^r for {l<j + k + l + m<N) (3.68) 



where M'^-f,^^ denote the elements of the reconstruction matrix. 
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3.4 Extraction of the Differential Cross Section 

The inclusive electron-nucleus cross section is extracted using the following formula: 



(PQ- _ PS X Nevents . ^ 

dVLdW^ ~ Acc X Anx AW^ x Nin x Eff x LiveTime x Ntg ^ ' ' 



where 

• PS is the prescale factor. 

• '^events IS the numbcr of scattered electron observed in solid angle dfl and dW^ 
range. 

• Acc is the acceptance of the dQ and dW^ bin. 

• AW^ is the width of bin. 

• AQ is the solid angle. 

• Nin = Charge/Q electron is the number of incident electrons. 

• Eff is the product of all efficiencies, trigger, particle identification, and tracking, 

• LiveTime is the product of electronic and computer live times. 

• Ntg=pNAd/A is the number of scattering centers per cm^. 

3.4.1 Acceptance Cuts 

The Hall C event reconstruction program reconstructs track parameters at the focal 



plane and using an optic model of HMS finds the target parameters using the Eq. 3.68 
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At the focal plane the track parameters are the x and y coordinates, and slopes of the 



track are and Xj^. These parameters are defined in Section 3.3 Cuts are applied 
to the reconstructed target quantities in order to eliminate events that are outside 
of the spectrometer acceptance but end up in the detectors after multiple scattering 
in the magnets or shielding. The acceptance cuts are given in the Table [6} The 



HMS acceptance cuts 




< 0.08 rad 


Wtar 


< 0.04 rad 


\6\ < 8% 



Table 6: Cuts on HMS reconstructed tracks. 

angle cuts are applied to eliminate events coming from outside the HMS acceptance. 
These events are the result of multiple scattering in HMS magnets and shielding 
material. The angle cuts are selected to be large enough to allow most events. The 
HMS momentum cut is applied to limit the HMS momentum acceptance since the 
reconstruction matrix elements could provide reliable tracking reconstruction only 
within these limits. 

For the HMS, the x'^^^ and y[^^ cuts typically rejected less than 0.3% of the total 
tracked events, and never more than 1.5%. For those rejected events, most come from 
events that are outside of the spectrometer acceptance but end up in the detectors 
after multiple scattering. The loss of these events is compensated by including mul- 
tiple scattering effects in the acceptance calculation. Also, these events are mostly 
lost from the very edge of the 6 bins. In the cross section calculation this edge bin is 
eliminated. The percentage of events lost due to not using the edge bins depends on 
kinematic setting and is in 0.1-0.5% range. 
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3.4.2 Particle Identification Cuts 

As in all experiments there are events coming from reactions other then the one the 
experiment wants to study. These events are regarded as a background. Different type 
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Figure 30: HMS calorimeter track energy hsshtrk/E' vs number of Cerenkov photo- 
electrons. The magenta line shows the Cerenkov cut, the red line shows the Calorime- 
ter cut. 



of detectors are constructed to decrease the number of background events in order to 
achieve the experiment's specific precision goals. For this experiment there are two 
detectors used for particle identification, a threshold Cerenkov detector described in 



Section 2.7.3 and an electromagnetic lead-glass calorimeter described in Section 2.7.4 



The main background in this experiment is from negative pions produced by charge 
exchange reactions. The ratio of pions to electrons varied from 0.1 to 30 for all 
runs. In this experiment the Cerenkov cut required that number of photoelectrons be 
bigger than 2 {hcerjnpe > 2), and the ratio of the Calorimeter track energy divided 
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to particle energy (determined from HMS) be greater than 0.7 [hsshtrk/hse > 0.7). 



Fig. 30 gives event distribution versus the number of photoelectrons on one axis and 
the track energy against the other. The pions and electrons are clearly identified. 



Cerenkov Cut EfRciency: Fig. 19 shows a HMS Cerenkov spectrum. The mean 
HMS signal is around 10 photoelectrons and most of the pions have zero photoelec- 
trons. The majority of the pions that have a signal with more than 2 photoelec- 
trons are pions that produce delta electrons at the front window or in the gas of the 
Cerenkov. These delta electrons emit Cerenkov light and the pion is misidentified as 
electron. In order to reduce the number of pions, a cut on Cerenkov signal is applied. 
The cut requires that number of photoelectrons to be greater than 2. This number 
is chosen to eliminate most of the pions and at the same time allows the majority 
of the electrons to pass. Not all electrons can be eliminated due to this cut so it is 
necessary to estimate how many. It is possible to estimate the number of electrons 
lost to this cut by looking at the total energy distribution in the calorimeter when 
applying all cross section cuts except the Cerenkov cut. This Cerenkov cut requires 
that number of photoelectrons to be in the range to 2, where is excluded to 



reduce number of pions in the energy distribution. In the top plot of Fig. [31] the 
ratio (hcaLet/hsshtrk) of total Calorimeter energy to HMS momentum is shown af- 
ter applying the aforementioned cut. The blue shaded area is the pion background 
estimated by a polynomial function. The parameters of the polynomial are obtained 
by fitting the distribution of hcaLet /hsshtrk ratio with a Cerenkov photoelectron cut 
equal to zero. The red hatched area is the estimated number of electrons that were 
cut by the requirement that the number of photoelectrons was greater than 2 and 
are the source of the Cerenkov inefficiency. The number of these events divided by 
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Figure 31: First plot: HMS calorimeter total energy hcaLet/E' distribution when 
number of photoelectrons are higher than but less than 2. Second plot: The 
Cerenkov cut efficiency as a function of scattered energy. 

the total number of electrons is the inefficiency of the Cerenkov cut and is shown in 



the bottom plot of Fig. 32 The Cerenkov cut efficiency is 99.7% and drops a little 
at the lowest scattered energy (most likely caused by high pion/electron ratio). The 
electron cut efficiency is parametrized in E' and used to correct the cross section. The 
systematic error of the Cerenkov cut is estimated to be about 0.1%. 
Calorimeter Cut EfRciency: The Cerenkov cut alone is not enough to achieve the 
necessary pion rejection power at low energies ( 0.4 — 1.0 GeV). An additional cut, 
hsshtrk/hse > 0.7, on the Calorimeter energy is applied to reduce the number of 
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Figure 32: HMS central momentum is 0.71 GeV. Top plot: HMS calorimeter track 
energy Etrack/E' (hsshtrk/hse) distribution without Cerenkov cut (the blue line) and 
with Cerenkov cut > 2 (the red line). Bottom plot: The Etrack/E' distribution after 
Cerenkov cut > 2 and Etrack/E' > 0.7 cut (the red hatched area). The solid blue 
area is the pion contamination. 

pions that pass the Cerenkov cut. The distribution of hsshtrk/hse, where hsshtrk is 
the track energy in the calorimeter and hse is the energy of the particle measured in 
the HMS, has electrons peaked at 1 and pions peaked at 0.3/ Etrack, see the bottom 



plot of Fig. 32 The hsshtrk/hse > 0.7 cut will remove most of the pions, but it 
will also remove some good electrons, which for some reason have less energy in the 
calorimeter than expected. It will cause an inefficiency (Calorimeter cut inefficiency) 
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and needs to be estimated. There is another source of inefficiency associated with the 
calorimeter cut called Wrong Tracking Efficiency (WTE). The WTE is not specific 
to calorimeters, but is rather a result of imperfect tracking algorithm. The WTE is 
absorbed into the calorimeter inefficiency since it is the result of the calorimeter cut. 
At the lowest scattered electron energy as this experiment, 0.44 GeV, the 



Calorimeter resolution is the worst and is equal to 7% (see Fig. 29). Even in this 
case the hsshtrk/hse > 0.7 cut is 3cr away from the peak and is more than 99.9% 
efficient. However, contribution from calibration errors, light collection inefficiency, 
and pedestal drift can leave some electrons below the the 3a threshold and will be 
cut by the hsshtrk/hse > 0.7 cut. 

To estimate the number of this electrons the following procedure is followed: 

• First, it is necessary to find the lowest value of hsshtrk/hse for good electrons. 
A few elastic runs with several E' are analyzed to estimate this value. It was 
found that small number of electrons, after a strict cut ( > 10 ) on the Cerenkov, 
can have hsshtrk/hse as low as 0.3. 

• Second, the cut efficiency is calculated with this formula 



cut _ Ar^^^"*^(E*°*^7E' > 0.7, hcaLel/hse > 0.2, Cerenkov > 10) 
~ N''vents(^El°\''J/E' > 0.3), hcaLel/hse > 0.2, Cerenkov > 10) 

where hcal_el is the total energy deposited in the first layer of the calorimeter. 
Since pions are less likely to deposit significant energy in the first layer (the 
first layer acts as preshower) this cut further reduces their number. In addition 
to the cuts in the above equation, acceptance cuts are also applied both in 
numerator and denominator. 



The procedure described above will count some pions as electrons even with strict 
cuts on the Cerenkov and on the first layer of the calorimeter. This means that it will 
underestimate the calorimeter cut efficiency. For each momentum setting there are a 
few runs with different angles, and therefore with different vr/e ratio. If there are no 
pions after El'^if/E' > 0.3 cut one would expect that the Calorimeter cut efficiency 
to be the same. In order to estimate the systematic error these pions introduce, the 
efficiencies from runs with the same momentum but different angles are averaged and 
the average is used as the cut efficiency for that momentum. The error bar is the 
standard deviation from the average. 

There are some good events (having passed the acceptance and the 
Cerenkov > 2 cuts) for which a cluster in the calorimeter is not found. These events 
can be seen in the spectrum of hsshtrk/hse around zero, as seen in the bottom plot 



of Fig. 32 For these events the tracking algorithm was unable to predict the correct 
position of the track on the front face of the Calorimeter, that is a wrong track energy 
reconstruction. However, the value of the total deposited energy in the calorimeter 
divided by the energy of the detected particle is still around 1 {E*°^i^ / E' ^ 1). The 
loss of these good events is compensated by calculating the efficiency due to the wrong 
tracking, denoted as e^^°'"-9 , The is calculated using this formula 



^wrong 



prevents (^^track^^f > 0.7, Cereukov > 10, utracks = 1) 



cai ^events (^E^otai I ^, > Q_7^ Cercukov > 10, fitracks = 1) ■ 

The requirement of having one track is essential to ensure that there is only one cluster 
in the calorimeter. At high counting rates it is possible to have two clusters produced 
by two high energy electrons. These two electrons will be counted as one and the lost 
electron will contribute as electronic dead time. In addition to the cuts in the above 
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Figure 33: Upper plot: Calorimeter efficiency versus E'. Second and third plots: 
Residuals of the fit when E' <1 GeV (left) and E' >1 GeV (right). 

equation, acceptance cuts are also applied both in numerator and denominator. 
The total Calorimeter efficiency is the product of WTE and Calorimeter cut ef- 



ficiency and is shown in Fig. 33 The top plot is the Calorimeter total efficiency vs 
E'. It is parametrized as a function of E' and used in the cross section analysis. The 
second and third plots are the residuals of the parametrization for E' < 1 GeV and 
E' > 1 GeV. The systematic uncertainty is estimated as the a of the residuals (left 
and right plots) and is used in cross section analysis. Systematic uncertainty is about 
0.13% for E' < 1 GeV and 0.03% for E' > 1 GeV. The difference is the result of higher 
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vr/e ratio at low energies. 

3.4.3 Background from the target walls 

Backgrounds are unfortunate part of all nuclear physics experiments. Depending on 
the goals of the experiment they must be taken into account by subtracting back- 
ground events or by estimating them and assigning a systematic uncertainty to the 
final measurement. In this experiment all backgrounds that can contribute more than 
half a percent are estimated and subtracted to obtain the final result. 




-5 5 

Y target (cm) 



Figure 34: Normalized yield distributions for deuterium target. The solid red line is 
the total normalized yield, the blue solid area is the contribution of estimated events 
(from dummy target) produced by aluminum walls of the liquid deuterium target, 
the cyan shaded area is the corrected deuterium yield. 

The cryotarget liquid targets have aluminum walls which can be source of back- 
ground induced by scattering on nucleons of the aluminum nuclei. In order to deter- 
mine the magnitude of this background, data on the aluminum dummy target have 
been taken at exactly the same kinematic settings as the deuterium data. The size 
of this background is estimated by the following formula 
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— V —V -'-Walls J^Alum. (n ^rV\ 

^corrected ^ Cryot. ^ Alum, rp -next ' yO.IUj 

J- Alum. J^Walls 



where, Y are the yields for the cryo and aluminum dummy targets, T is the thickness 
of the aluminum walls in the cryotarget and aluminum dummy target, and R is the 
external radiative correction for the aluminum dummy target or the aluminum walls 
in the cryotarget. Yield for one of the deuteron runs without aluminum background 



subtraction is shown in the Fig. |34j In the same figure the yield for the target wall 
background and the corrected deuterium yield are shown. 

The size of the background from the aluminum target walls is generally about 
10%. The background has been calculated for each run and the subtraction has been 
done on a bin-by-bin basis in (W^^, 9). 

3.4.4 Charge Symmetric Background 

In this experiment inclusive electrons are scattered from nuclear targets. The prob- 
ability of production of 7 and 7r° particles in the target is comparable to the cross 
section of inclusive electron scattering cross section in the kinematic range of this 
experiment. 

The 7r° decays predominantly to two photons which produce electron-positron 
pairs when passing through the target or any material between target and the drift 
chambers. In the beam energy range of the current experiment bremsstrahlung of 7 
particles (virtual and real) is the dominant electron energy loss mechanism. The 7 
particles will produce electron-positron pairs in the field of a nucleon (Bethe-Heitler 
process). The resulting electrons produced through the reaction mechanisms de- 
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scribed above is kinematically indistinguishable from inclusive electrons and passes 
all particle identification cuts. These electrons are a background and need to be 
subtracted from real inclusive electrons scattered from the target. Since the normal 
configuration of the experiment is unable to remove these electrons, another approach 
is used. From the production mechanism of background electrons it can be seen that 
there should be equal number of positrons and electrons (that is why this is called 
charge symmetric background (CSB)). Thus, the polarity of the HMS spectrometer 
is reversed and positron data are taken at the spectrometer settings where the charge 
symmetric background is large. This background is significant for larger scattering 
angles and small final electron energies of the electrons. 




,= 2.095 GeV 
,= 3.116 GeV 
,= 4.074 GeV 
,= 4.134 GeV 



0.5 1 1.5 

E' (GeV) 

Figure 35: Cross section of Charge Symmetric Background for the carbon target. 

The background contribution from kaon decay is negligible and is not taken into 
account, for more details see Ref. |75]. Runs with reversed HMS polarity were taken 
where the CSB was considered to have significant contribution. The CSB cross section 
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Figure 36: First plot: Ratio of CSB cross section to raw electron cross section versus 
E' for all carbon runs. Second plot: The residuals of CSB fit for ^=25,35,45 degrees 
(left) and 6'=60,75 degrees (right). 

is calculated for all these runs and then extrapolated by the following function: 



where E' is the electron energy. C{6) and S{6) are free parameters determined from 



fitting the cross sections for each beam energy and angle as shown in Fig. 35 The 
form of the function is chosen to go to zero at E' equal to the beam energy and 
have smallest possible number of parameters. After fitting the CSB cross sections at 
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each energy and angle, the C{6) and S{6) parameters are interpolated in theta by a 
polynomial function of second degree. The parameters of the polynomials are stored 
in a lookup table and are used in the cross section analysis. The CSB is subtracted 
from the total electron cross section for each bin {6, W"^). 

The CSB reaches 60% at the smallest values of E' = 0.44 GeV and the largest 



angles, 9 =76° and about 50% for 6 =75°, as can be seen in the top plot of Fig. |35 
The point to point uncertainty has been estimated as the deviation of the measured 
cross section from the parametrization and is equal to 0.1% for angles less than 45° 



and 0.4% for angles 60° and 75°. See the bottom plots of Fig. 36 



3.4.5 Pion Contamination 

The Cerenkov and the calorimeter cuts can't fully remove all pions. Most of the pions 
do not give a signal in the Cerenkov detector, but some pions can produce knock-on 
electrons of high enough energy ( more than 20 MeV) to emit Cerenkov light. These 
pions will pass the Cerenkov cut. In the calorimeter, pions can be registered due to 
ionization losses. This is visible as a wide peak around 0.3/i?'(GeV) in the normalized 
energy spectrum of the calorimeter. Also pions can undergo charge-exchange reactions 
and produce a neutral pion which decays into two photons. These photons will deposit 
their full energy in the calorimeter which leads to a high energy tail in the normalized 
energy spectrum of the calorimeter and can go beyond Ecai/E' = 0.7 ( the calorimeter 
cut). As a result both Cerenkov and calorimeter cuts will fail to reject some of 
the pions and they will be accepted as electrons, resulting in pion contamination 
in the electron sample. In the kinematic range of the current experiment the pion 
contamination is estimated to be up to 3%. Investigation showed no rate dependence 
of pion contamination. 
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3.4.6 Electronic and Computer Dead Times 



In this experiment electronic dead time was low (less than 3% ) since the highest 
counting rate is less than 500 kHz. Computer dead time reached up to 60% for very 



few runs and is generally lower than 30%; see Fig. 37 
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Figure 37: First plot: Computer live time as a function of rate. Second plot: 
Electronic live time as a function of rate. 



Electronic Deadtime: Electronic dead time (EDT) is caused when a trigger is 
missed because the hardware is busy when an event that should generate a trigger 
comes in. Detector dead time occurs when a detector is unable to respond to an event 
because it is still responding to a previous event. 

For an average event rate R the probability of finding n counts in a time interval 
T is given by the Poisson distribution, 
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-Rt 




(3.71) 



The live time is the probabihty of no events occurring in the interval r is P(0) = e"^"^. 
For small Rt this can be approximated by P(0) ft! 1 — Rt. An event will be missed 
when it arrives within a time r of an event accepted by the gate, where r is the gate 
width of the logic signal. Therefore, the fraction of the measured events is equal to 
the probability that the time between events will be greater than r. 



In order to determine NtotaU scalers with four different gate widths (i=50, 100, 150, 
200 ns) have been used and the number of events recorded in each of them have been 
measured (A^5o, -/VioO) -^i50) -^200)- A linear extrapolation back to zero gate width 
gives Ntotai- In Hall-C ENGINE the electronic dead time is calculated using only 



measured 



1-Rt . 



(3.72) 



X 

^ 250 



260 



Pretriggers(Ax=0) = 248.2 
Pretriggers(Ax=60) = 242.6 
Elec dead time = 2.3% 



0) 





220 







50 



100 150 



200 



Gate width (ns) 



Figure 38: Rates with different gate widths for run 63915. 
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two pretriggers, PREIOO and PRE150. The following formula is used to get the final 
electronic dead time (EDT) for each run. 



^^^^ (PREIOO -PREm At 

PREIOO 50 ^ ^ 



where r is the real electron trigger gate width (60 ns). 

All runs are corrected on a run by run basis, and no systematic uncertainty is 
assigned to the electronic dead time. For more information about electronic dead 
time see Ref. [7T] . 

Computer Deadtime: The computer dead time occurs as a result of the data 
acquisition computers being busy processing an event and not being available to 
process new events. In that case a new event is lost. Events recorded by the electronics 
have been logged as pretriggers. The pretriggers that have been processed successfully 
by the trigger supervisor are recorded as triggers. The ratio of triggers to pretriggers 



gives the computer live time (see Fig. 37) 



L 



trigger 



^pretrigger 



The computer live time is calculated for each run and applied on a run-by-run basis. 
The systematic uncertainty for the computer live time has been studied by compar- 
ing experimental yields at the same spectrometer kinematic setting but varying the 
prescale factors and therefore, computer live times [75]. These studies show that the 
yields agree within 0.2%. 



3.4.7 Trigger Efficiency 

Triggers are used to reduce the rate of useful events into a range manageable by 
the data acquisition equipment. Triggers also provide timing signals to the various 
detector parts. 

Some events can be lost due to inefficiencies of the detectors which are used to 
form different triggers of the HMS spectrometer. The schematics of the single arm 



trigger logic is shown in Fig. 22 It can be seen that the two electron triggers (low - 
level and high - level) are strongly correlated. For example, if the (SCIN) signal (at 
least three of the four scintillator layers of both hodoscopes have fired) is present for 
an event, there must also be a (STOF) signal (which requires one front panel and 
one back panel). If (PRHI) is present then there must also be a (PRLO) signal (high 
(PRHI) and low (PRLO) threshold on the energy in the first layer of the calorimeter). 
Thus, assuming that the Cerenkov signal (C) is always present, the efficiency for the 
low- level electron trigger can be calculated as epRio x ^stof- The small inefficiency 
of the Cerenkov is covered by the ELHI trigger which does not require the Cerenkov 
signal. If the calorimeter has an inefficiency, an electron trigger still could be produced 
if there is a (SCIN) signal. Accordingly the electron trigger efficiency for the HMS 
can be calculated (estimated) by the following formula, 

^trg = ^PRLO X ^STOF + (1 " ^PRLo) X ^3/4 , (3.74) 



where €3/4 is defined in Eq. 3.75 

The scintillator 3/4 efficiency (63/4) has been calculated by summing over all pos- 
sible combination that would satisfy the 3/4 trigger 
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Figure 39: HMS electron trigger efficiency as a function of the scattered electron 
energy. 



£3/4 = n + 5Z ~ n ' 

1=1,4 j=l,4 i^j 

where i, j are the plane numbers and e is the corresponding efficiency for each plane. 
The average scintillator 3/4 efficiency (63/4) is about 0.983. 

The (PRLO) efficiency {epRLo) has been calculated as the ratio of the events that 
have a (PRLO) signal and a signal from the Cerenkov over the events that have a 
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signal from the Cerenkov, 



PRLO & CER > 0.5 , ^ , 



The average (PRLO) efficiency is about 0.98. 

The (STOF) efficiency (cstof) has been calculated by using the efficiency of each 
scintillator (ei, €2, 63, €4,) by the following formula 

esTOF = fl - (1 - ei) * (1 - £2)) X fl - (1 - £3) * (1 - 64)) . (3.77) 



The average (STOF) efficiency is about 0.999. 



Finally, the total electron trigger efficiency for the HMS is calculated by Eq. 3.74 



and the result is shown in Fig. 39 The trigger efficiency is always higher than 
0.999. The trigger efficiency has been parametrized as a function of the scattered 
electron energy and during the data analysis this parametrization has been used 
on a run-by-run basis. The systematic uncertainty is estimated as the spread of 
the measured trigger efficiency from the parametrization and is equal to 0.007%. A 
detailed description of trigger efficiency calculation can be found in Ref. [76] . 

3.4.8 Tracking Efficiency 

Sometimes track reconstruction can fail even when there is a legitimate track passing 
through the detector system. These lost events should be taken into account during 
cross section calculations. Lost tracks happen when the tracking algorithm fails to 
reconstruct a track. There are reasons for this: ffist the ability of hardware to handle 
high rates and work effectively, second, the prohibitively high computer time needed 
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Fiducial Areas 




Figure 40: Illustration of the fiducial areas on the scintillator hodoscopes used for 
determination of the tracking efficiency. 



to reconstruct an event when there are multiple hits in drift chambers caused by 
background, and third, at high rates it is possible to have two legitimate tracks and 
the tracking algorithm is not able to reconstruct either of these events. In order to 
minimize the track reconstruction time the Drift Chamber TDC window is chosen 
to be 250 ns, it is required that each chamber have no more than 25 hits, and the 
maximum number of focal plane tracks should be less than 10. In order to calculate 
the tracking efficiency, only the central part of the HMS acceptance is chosen (a 
fiducial region). The fiducial area consists of paddles 4-13 in the X planes and 4-7 in 



the Y planes, as shown in Fig. 40 This defines the fiducial area, and guarantees that 
only one particle passed through both drift chambers within the central area of the 
spectrometer acceptance and should have been tracked. 

In order to calculate the tracking efficiency a clean sample of electrons is selected 
using electron identification cuts. The tracking efficiency is calculated as the number 
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Figure 41: Tracking efficiency as a function of tlie scintillator (3/4) rate. 

of events for which a track has been found, divided by the number of total events 
passing the electron identification cuts as, 



Trigger & PID & Track 
Trigger & PID 



^(tracking) — rn„.- t->tt^ ' [0.(o) 



where Trigger indicates that there is a trigger. Track is that at least one track is found 
and PID represents the particle identification cut. 



PID — hcerjnpe > 2.0 & heal jet / hpcentral > 0.7. 
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The tracking efficiency calculated for each run as a function of the scintillator (3/4) 



rate is shown in Fig. 41 The tracking efficiency has been parametrized as a function 
of the scintillator (3/4) rate and during the data analysis this parametrization has 
been used on a run-by-run basis. The systematic uncertainty is estimated as the 
spread of the measured tracking efficiency from the parametrization and is equal to 
0.15%. For a detailed description of tracking see Ref. [76] . 

3.4.9 Acceptance Calculation 

Knowledge of the HMS acceptance is one of the dominant sources of uncertainties in 
determining the cross section. For the HMS the acceptance is normally defined by a 
collimator and subsequent magnet apertures. As the magnet apertures partially define 
the acceptance, the magnetic model must be known in order to find the acceptance 
function. 

For a given angle and momentum the HMS can detect particles which have angles 
and momentum around the given values. The HMS acceptance is a function of the 
three target coordinates X, Y, Z and three spectrometer coordinates S, X', Y'. For 
thin targets the cross section is independent of X, Y, Z and therefore the acceptance 



is a function of only 6, A', Y' . Spectrometer angles A', Y', defined in Section 3^ 
can be related to the polar angle 6 by the following formula 

9 = arccos(cos(A') cos{9hms — ^')) (3.79) 



where Ohms is the central angle of the HMS. 

Since the inclusive cross section is independent of azimuthal angle (f), the accep- 
tance is a function of only two variables, A{6, 6) where 6 is the momentum fraction 
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9 and is the polar angle . Before calculating the acceptance the solid angle of each 
9 bin is calculated. During the cross section analysis events are binned in twenty 9 
bins from —0.035 to +0.035 radian. For each bin the solid angle is calculated from a 
Monte Carlo. First, events are generated within the range \X'\ < 0.1 and \Y'\ < 0.1 
and than the number of events that passed through the acceptance of the HMS are 



compared to the number of initially generated events. In Fig. |42]the thickness of the 
blue region is the width of the 9 bin (0.0035 radian). If the aperture is determined 
by the collimator only, the solid angle will be the blue region, but since the HMS has 
other apertures, the real solid angle is determined in the following way. The solid 
angle for that given 9 bin is calculated as the ratio of number of events in red region 
to number of events in the blue region (events are sampled within \X'\ < 0.08 and 
\Y'\ < 0.04) multiplied by total solid angle of AX'Y'. 

To calculate the acceptance function A{E', 9) events are generated in X, Y , Z, 5, 
X', Y' and binned in 9. In general it is convenient to use W"^ bins instead of 5{E') 
bins. There are several advantages in using bins over 5 bins. For example, in the 
kinematic range where the cross section varies significantly the density of W'^ bins is 
higher than the density of 5 bins. Also the radiative corrections are calculated at the 
centers of the W'^ bins, reducing the model cross section interpolation errors. 

In this analysis W'^ binning is not done in a direct way. Instead of binning the 
data in W'^ bins, data are binned in E' bins in such a manner that the lowest and 
the highest E' bins are within the known range of the HMS momentum acceptance 
of < 8 and at the same time correspond to a point in the W"^ grid. Here the W"^ 
grid points are given by this formula Wf = 0.05 + dW"^ x i, where z = 0, 150 and 
dW'^ = 0.04 GeV^. This binning method is chosen over the direct binning based 
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Figure 42: Solid angle calculation for each theta bin. The blue region is the solid 
angle of a theta bin when there is no aperture, the red region is when the aperture of 
HMS is applied. The solid angle is equal to the ratio of events with aperture on and 
off multiplied by the total solid angle of AX'Y'. 

on the following consideration: the acceptance does not always include a symmetric 
region in 6 about the central value in W'^, for bins corresponding to high or 
low values of 6 at the central angle only part of the 6 acceptance lies within the 
spectrometer acceptance. The result is the maximum bin centering corrections are 
at the edge of the acceptance, where the acceptance is more poorly known than in 
the central part. Therefore, for W"^ bins at the edge of the momentum acceptance, 
the systematic error is larger than for W"^ bins in the central part. Also, taking into 
account the fact that there are less events available at the edge of the acceptance, the 
edge bins will have bigger statistical errors. 

There are three main elements in the Monte Carlo simulation code: the event gen- 
erator, the transport of the particles through the magnets, and the list of materials 
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and apertures that cause multiple scattering or stop the particles. The initial coordi- 
nates are randomly generated along the target length, while the quantities 6, 9 and 
are chosen randomly within their allowed limits. The particle is transported through 
the spectrometer to the detector hut using the computer program COSY Infinity ^78j , 
which models the magnetic part of the spectrometer, the magnet positions, their in- 
ternal dimensions and their magnetic field maps. If the particle successfully traversed 
the spectrometer and the detector stack, it is considered a success and contributes to 
N success- The spectrometer acceptance is calculated by the following formula. 



Acceptance{e. - B^. E\) = ^!7'7f ^ /"f^'^ , (3.80) 

^^gen[oi — do, ^j) 



where N^enl^'i — do,Ej) and ^success{(^i — (^o,Ej) are the number of events generated 
and detected in a given {9i — 9o, E'-) bin, respectively. In the cross section analysis 
data are binned in {9i — 6o,Ej) bins and each bin is corrected by the corresponding 



Acceptance{6i — 6o,Ej). In Fig. 43 the HMS acceptance is shown for Phms = 0.44 
GeV and 9 =Tb°, obtained by generating five million Monte Carlo events. One can 
see that the momentum range used in this analysis {\5\ < 8%) is well within the 
fiat range, while the angular acceptance falls off quickly at the age of the angular 
acceptance. This fall off is due to the shape of the collimator. 

During the cross section analysis a systematic dependence of the experimental 



cross section divided by the model cross section was observed and is shown in Fig. 44 
After studying this effect under different kinematic conditions and offsets no target 
or kinematic E' , 6) dependence was found. This effect is believed to be caused 
by mis-calibration of the optics of the HMS. In order to correct for this effect a set 
of runs are selected at different HMS central momenta and angles. These runs are 
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Figure 43: HMS Acceptance function at Phms = 0.44 GeV and 6 = 75 degree, 
calculated by HMS Single-Arm Monte Carlo program. 

selected to have good agreement with the model cross section in order to rule out 
any model dependence in the extraction of the correction function. The cross section 
ratios are parametrized as function of 6p and are used in cross section analysis in 
bin-by-bin basis. The uncertainty associated with this parametrization is ~ 0.6%, 
estimated from the spread of data to the fit. 

In addition, the position uncertainties on the target, collimator, magnets, and 
detector package contribute to both the point-to-point and normalized acceptance 
uncertainties. The total point-to-point uncertainty due to the position uncertainties 
on the target, collimator, magnets, and detector package is ~ 0.3%. As a result, 
the total point-to-point (kinematic dependent) uncertainty on the acceptance correc- 
tion was ~ 0.7%. The normalized (kinematic independent) uncertainty on the HMS 
acceptance is ~ 0.6%, which is determined by changing the positions slightly of the 
collimator, magnets, and detector package respectively in the HMS Single- Arm Monte 
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6 p (%) 

Figure 44: Ratio of experimental to model cross sections versus 5p. The red line is 
the same ratio obtained in a previous experiment, E99-118. This effect is caused by 
optic mis-calibration. 

Carlo code and checking how this affected the extracted cross sections. 
3.4.10 Bin Centering Corrections 

The goal of the analysis is to extract electron-nucleus inclusive differential cross sec- 
tion for a range of E' at a fixed scattering angles. During data acquisition each run 
is taken at a particular HMS angle and momentum. Since the HMS spectrometer 
has a relatively large angular acceptance of ± 2 degrees, the cross section can vary 
greatly across the angular acceptance. In order to calculate the cross section for each 
E' bin at fixed 9 value, a 9 bin centering correction must be done. If the cross section 
varied linearly across a symmetric 9 acceptance there will be no need for bin centering 
correction. 

This correction removes the cross section dependence in the angular acceptance 
range using a known model cross section. The bin centering correction depends on 
the model cross section used, and the model dependence in this correction can be a 
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large systematic uncertainty in the analysis. The correction is applied by rescaling 
each {E', 6) bin by the ratio of its cross section at the center of the bin to the central 
cross section {E',6o). The bin centering correction for {E[^9o) bin is calculated with 
this formula: 

BC{E, E\ 6) = aZTiE, E[, ^-"*™') / aZt\E, E^ 9), (3.81) 



where E'^ is the center of the E' bin, cr^^'^^' is the total radiated cross section. 



3.4.11 Radiative Corrections 

The differential cross sections measured in the resonance region for electron-nucleon 
scattering may have large contributions from processes other than the one photon 
exchange approximation. This is also called Born approximation and is shown in 



Fig. 45 The Feynman diagrams for higher order electromagnetic processes in a are 
shown in ( 2, 3, 4, 5, ) which include vacuum polarization (creation and annihilation 
of particle-antiparticle pairs), vertex processes (emission and reabsorption of virtual 
photons), and bremsstrahlung (emission of real photons in the field of the nucleon 
during interaction). 

The cross section measured in the kinematic range of this experiment may have 
up to a 30% contribution from those processes. In order to determine the differential 
cross section for the one— photon exchange process, all the other contributions from 
the higher order processes in a have to be estimated and corrected for in the measured 
cross section. 

The impact of the radiative process to experimental data has the following three 
effects: 
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1) Born 2) Vacuum polarization 3) Vertex correction 




4) Bremsstrahlung 5) Multi-photon emission 



Figure 45: 1) is the lowest-order diagram for charged-lepton-nucleon scattering. 
2)-4) are the lowest-order electromagnetic radiative corrections. 

• The incoming particle's momentum may be reduced, thereby changing the kine- 
matic configuration for the event, and also the probability for the interaction 
(through the energy dependence of the cross section). 

• The outgoing particle's momentum may be reduced, moving events from one 
kinematical configuration to another. 

• The overall cross section for the process in a particular kinematical configuration 
will change. 

These effects are taken into account by applying radiative corrections to the mea- 
sured cross section. Electromagnetic radiative effects in electron scattering are divided 
into two categories: external and internal. The external radiative processes, shown 



in Fig. 46 , can take place before and after the scattering in the material the electrons 



passes through (the largest contribution happens in the target material). Both cases 



Ill 

have to be corrected in order to calculate the Born cross section. Corrections for 
ionization energy losses are also applied, but have a much smaller contribution than 
the radiative processes. 





Figure 46: The Feynman diagrams for external radiative processes. 

Internal effects occur at the scattering vertex and are calculable in Quantum 
Electrodynamics. The internal effects include internal bremsstrahlung (emission of 
photons in the field of the nucleon during the scattering process), vacuum polarization, 
vertex processes, and multiple photon exchange (important at low Q"^). In the first 
order in a there is also contribution from the diagrams of vacuum polarization by 
hadrons, which is not shown here but is taken into account in Ref. [79] . 

The program used to calculate radiative corrections for this experiment, including 
both the internal and external effects and the elastic tail, are based on the program 
developed at SLAC which is described in Ref. [80] . 

The measured cross section is the combination of radiated elastic, quasi-elastic 



and inelastic cross sections, are depicted in Fig. 47, and the total measured cross 
section is the sum of all three, 



and in inclusive experiments the inelastic and quasi-elastic cross sections are not 
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1) Elastic 



P_ 

A A ^ ^ A-1 X 

2) Quasi-elastic 3) Inelastic 




Figure 47: Feynman diagrams of processes present in the kinematic range of this 
experiment. 

separated. 

The total radiative cross section (quasi-elastic + inelastic) is proportional to the 
Born cross section, and can be calculated by the following formula: 



inel+qe 
Born 



rad 



a. 



mod 
rad 



a 



rad 
el 



X a 



mod 
Born 



(3.83) 



where cr™^'^ is the total radiated model cross section (internal + external ), cr™'^ is the 
elastic radiated model cross section (internal + external ). The radiated model cross 
section is the convolution of internal and external radiative cross sections: 



= internal ® external ® Born. (3.84) 

For this experiment, the external radiative corrections are computed using a complete 
calculation of Mo-Tsai |HI] with a few approximations. This approach, MTEQUI, uses 
the equivalent radiator approximation [80]. In the equivalent radiator method, the 
effect of internal Bremsstrahlung is calculated using two hypothetical radiators of 
equal radiation length, one placed before and one after the scattering. It is important 
to note that the energy-peaking approximation is not used for the computation of 
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external contributions. The internal contribution in MTEQUI method is evaluated 
by setting the radiation length of the material before and after the scattering point 
to zero, and ignoring the target length integral (see Eqn. CI in Ref. [81j). 

The internal radiative corrections a^^f have been calculated using the formalism 



of Bardin [79]. The calculations were done for all diagrams in Fig. 45 without any 
approximation. 

In order to reduce the effects of any approximations in treating the external ra- 
diative effects we use the following expression to calculate a Born, 



ii=el,qe,inel 



ii=el,qe,inel 



^rad 
ii,int 



(3.85) 



\ ii=el,qe,inel / Mo,Tsai 



It was found that using Mo and Tsai's MTEQUI method to treat external radiative 
correction is consistent for targets with different radiation lengths, see Ref. |80j . 

The uncertainty on the cross sections due to the radiative correction is estimated 
at 1% for the e dependence, and the normalized uncertainty is also ~ 1%, according 
to the radiative correction studies done at SLAG [80j. The e dependent uncertainty 
and normalization uncertainty caused by the two-photon corrections to the one pho- 
ton exchange diagram are estimated to be ~ 0.3% respectively, according to the 
results of an early SLAG experiment |90]. That experiment measured the difference 
of electron/positron cross sections at Q^=1.2-3.3 GeV^ and W^=l. 3-17.0 GeV^ range 
(comparable to the kinematic range of this experiment) and found (Te+/o"e-=1.0027± 
0.0035 with no significant dependence on or W^. Since no significant difference 
was found between electron/positron cross sections in that experiment, one can con- 
clude that two photon exchange contributions are small in the kinematic range of this 
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(GeV) 



Figure 48: Systematic errors due to contribution from o? processes are shown in 
four different beam energies. The error is set to be 3% of the size of the quasi-elastic 
tail, relative to the inelastic [H2]- The error has a slight e (beam energy) dependence 
following Ref. [80]. 
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experiment. 

In the kinematic range of this experiment the radiative corrections are a maximum 
of 30%. Systematic errors due to contributions from second order in a processes are 



set at 3% of the size of the quasi-elastic tail, relative to the inelastic [H2], see Fig. 48 
The error is less than 1% for W"^ > 1.4 and has slight e (beam energy) dependence. 



The Born model cross section used for this analysis is discussed in Sec. |3.5[ The 
radiative correction is applied to measured cross sections on bin-by-bin basis to get 
the corrected cross sections. 

The recoil polarization measurements of the form factor ratio /^pC^/C^^ |83t l8l] 
contradict the Rosenbluth measurements (See Ref. |H5] for a compilation and refer- 
ences) and it has been suggested that the earlier experiments might have not have 
fully understood their systematic errors or had normalization problems. The Rosen- 
bluth measurements have been reexamined [BS] and this global reanalysis could find 
no systematic or normalization problems that could account for the discrepancy. The 
author of Ref. [86] concluded that a modest linear e-dependence correction (of origin 
yet unknown) to the cross section measurements might explain the difference. Sev- 
eral investigators [871 EHl EH] have explored the possibility of two-photon exchange 
corrections (which would be less important in the direct ratio measurement of recoil 
polarization) to explain the discrepancy. While only incomplete calculations exist, 
the results of Ref. [HH [89] account for part of the difference. 

As this experiment exploits the Rosenbluth technique the question naturally arises 
whether two-photon effects might play role. Unfortunately no calculations have been 
done specifically to answer this question. Nonetheless we have attempted to gauge 
what effect an e dependence might have in our results. To do so we artificially 
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introduced a 2% e dependence in the radiative corrections which resulted in ±0.028 



uncertainty in R, see Sec. AA. Hence we can conclude that for our specific conditions 
the neglect of two photon effects (as manifested in a e-dependence) has only modest 
consequences. 

3.4.12 Coulomb Correction 

In the field of high Z nuclear targets the incoming and outgoing electron wave func- 
tions are affected by the Coulomb field. From the classical point of view the incoming 
electron accelerates in the field of the positive charged nuclei and the scattered outgo- 
ing electron decelerates. This effect causes an increase in the momentum of incoming 
beam electron and a decrease in the momentum of the scattered electron relative to 
the vertex values. At high enough beam energies this effect is negligible while at the 
beam energies of this experiment it is not and has to be taken into account. The 
change of the vertex kinematic variables can have a big impact on the measured cross 
sections. The change of the energy of the incoming electron will change the cross 
section of the process of interest while the change of the scattered electrons energy 
change will cause it to populate a different kinematic bin and change the measured 
cross section. 

Coulomb corrections have the result that the plane wave Born approximation 
loses its validity and hence requires a correction. For this experiment the Effective 
Momentum Approximation (EMA) is used as described in Ref. [91] and discussed 
below. 

Assuming a spherical charge distribution in the nucleus, the electrostatic potential 
inside the charged sphere can be defined as followed: 
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V{r) 



3a{Z-l) a{Z-l)r 



(3.86) 



2R 2R R 



where the radius of A is given by this formula: 



(3.87) 



Because most the nucleons of heavy nuclei are located in the periphery of the nucleus, 
taking the electrostatic potential at the center of the nucleus will be an overestimate 
of the Coulomb effect. In the effective momentum approximation (EMA) [91j, the 
effective potential used is Veff ~ (0.75 — 0.8)V{r = 0). This value for Vg// agrees 
with the extracted effective potentials from positron and electron inclusive scattering 
experiments, see Ref. [92]. In order to take account the K//, is replaced by Qlff 



There are two methods to do Coulomb correction: one assumes that the Mott 
cross section is unchanged while the nuclear spectral function is subjected to a trans- 
formation when is replaced by Qe// defined above. The second method changes 
both the Mott cross section (Q^ — > Qlff and E' — > E' + Veff) and the nuclear spectral 
function. Also, for the second method a focusing factor for the incoming electron is 
given by the following formula 



where 




(3.88) 
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Figure 49: The Coulomb correction factor as a function of W^, for low to high e 
range and for carbon, aluminum and iron targets. For heavy nuclei, the Coulomb 
correction factor is significant, reaching a maximum of ~ 15% at low W"^, near the 
quasi-elastic peak. The Coulomb correction is smaller for high beam energy (high e 
range) than for low energy (low e range). 
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Ffoc = t (3.89) 



is applied. In this experiment the average potential of K// = 0.775 x V{r = 0, Z, A) 
is used. The Coulomb correction is applied to extracted cross section in the following 

way 



E') = a„^.(E, E') X .I'Tv':;''^ +V „) 

'^mode/l-^ ^ ^eff^-^ ^ ^eff) ^ foe 



The Coulomb correction factors versus are shown in Fig. 49 for beam energies 
2.1, 3.1, 5.1 GeV. For heavy nuclei, the Coulomb correction factor is significant, 
reaching a maximum of ~ 15% at low W'^, near the quasi-elastic peak. Different 
beam energies indicate a different range in e. The Coulomb correction is bigger for low 
epsilon (low beam energy) data than for high epsilon data (high beam energy). It is 
necessary to apply the Coulomb correction to properly extract the nuclear R=(Ti/aT. 

3.5 Cross Section Model 

The purpose of this experiment is to measure the Born cross section for inclusive 
electron nucleus scattering in resonance region. Since the exact Born cross section 
for inelastic electron-nuclei scattering is impossible to measure due to superposition 



of different reactions (see Fig. 47) present in the kinematic range of this experiment, 
a model cross section is necessary. Also a cross section model is necessary for the bin 
centering correction, radiative corrections, and the Coulomb corrections. 
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3.5.1 Quasi-elastic model 

In earlier sections the scaling of the deep inelastic cross section has been discussed. 
This refers to the observation that the cross section for scattering from a nucleon, 
normally of function of two kinematic variables, reduces to a function of a single 
variable, x = and reflects that the scattering takes place from a structureless, 
point like object that is essentially free within the nucleon. 

There exists a analogous scaling behavior in scattering from quasifree nucleons 
in the nucleus |93]. In the PWIA quasi-elastic scattering, where the nucleons are 
moving independently in the mean field of the nucleus, the cross section as a function 
of final electron energy E' is a convolution of the spectral function S{k, E) with the 
elementary electron-nucleon cross section cjei- The spectral function S{k, E) is the 
joint probability to find a nucleon with momentum k and separation energy E in the 
nucleus. Schematically then. 



where a summation over all the protons and nucleons is implied and the delta function 
argument conserves energy and momentum. 

At large momentum transfer the expression can be rewritten (with a few important 
assumptions, see Ref. [HS]) as 



with K a kinematic factor and F{y) is the longitudinal momentum distribution, y is 
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the scaling variable and is determined from energy conservation with v — E — E', 



u + Ma^ [{Ma-1 + Ea-iY + k^Y'^ + [M^ + (k + q)''Y'\ 



y is the minimum momentum of the struck nucleon satisfying the previous expression 

- the longitudinal momentum of the struck nucleon. 

Hence the quasi-elastic cross section from a particular nucleus is a function F{y) of 
a single variable, y, itself a function of q and u, and is independent of them separately 

- the cross sections scale. This is called scahng of the first kind. 

The difference between x and y scaling is that in y scaling the struck object has 
structure and must be accounted for through the division of the cross section by 
the elementary elastic electron-nucleon cross section (which is a strong function of 
momentum transfer). While useful for the modeling of the cross section, y scaling 
can provide access to the nucleon momentum distribution n{k) through F{y), the 
scaling function 



with F{y) the longitudinal momentum distribution. Hence quasi-elastic scattering 
would provide a direct measure of n{k). 




and if E, 



'max 



= OO then 




n{k) 



1 dFjy) 
2'Ky dy 
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In a series of papers by Donnelly and Sick [Oil ESI EHl EH] and collaborators it was 
found that scaling of the quasi-elastic cross section could be extended to account for 
the difference in atomic number A. This is scaling of the second kind. The scaling 
variable ip is given by the following formula 



where A, k, r are dimensionless energy, 3-momentum and 4-momentum transfers re- 
spectively (these variables are normalized by a factor l/rriN, where is the nucleon 
mass.). The new scaling variable ip is related to y but lacks one important property, it 
does not take account the small energy shift EsUft which is included in separation en- 
ergy Eg. An improved phenomenological dimensionless scaling variable is employed 
in treatments of superscaling [97], which included the empirical shift Eghift- One 
should note that scaling function is very sensitive to kp and yet it is possible to find 
a value for it to line up all data on one curve. The dependence from 'EisM ft is not that 
pronounced but it needs to be taken into account to put the quasi-elastic peak in the 
place where the scaling variable if) is zero. 

In order to fix this a new variable ip' is introduced, which has the same form as 
but the dimensionless variables A and r are shifted. The new scaling variable is given 
by the following formula 



1 



A -r 



(3.91) 




1 



A'-r' 



(3.92) 
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A 


EsMft (GeV) 


kpiGeV) 


Carbon 


0.020 


0.228 


Aluminum 


0.018 


0.236 


Iron 


0.018 


0.241 


Copper 


0.018 


0.245 



Table 7: Values of Eghift and kp for nuclear targets used in this experiment. 

where Xshift = Eshiftf^mN, X' = X - Xshifu t' = \/ {k^ - A'^)- The scaling function 
in ip' is given by the following formula: 



1.5576 



(1 + 1.7722(7/;' + 0.3014)2)(1 + e(-2-429iV''))A;^ ^ 



(3.93) 



see Ref. p8] for details. The values of E shift and kp are given in the Table (T) The 
structure functions Fi and F2 are calculated using the following formulas 



Fi = MpF(^')G'T/2 

(3.94) 

F2 = z/F(V'')KG'l + z/tG't) 

In the above equation the Gl and Gt are related to the nucleon elastic form-factors 
and their values are taken from the Bosted fit for the nucleon form factors |85]. Pauli 
suppression is taken into account according to Eq. B54 of Ref. • From Table [T] 
one can see that the energy shift does not vary too much for different nuclei. The 
values of kp don't vary too much after carbon. 

3.5.2 Inelastic Model 

The kinematic range of this experiment covers W"^ = 0.0-4.5 GeV^ and 
Q2 = 0.5-4.5 GeV^. The inelastic model cross section used in this experiment 
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is based on the structure functions of proton and deuterium, as well as on the 
EMC |100] effect measured in the DIS region. The proton structure functions are 
extracted from an inclusive inelastic electron-proton cross sections fit described in 
Ref. jlUlj . The fit is constrained by photoproduction data at Q^=0 and makes a 
smooth transition to DIS at higher values. The other important aspect of this 
fit is that it is based on R = ai/cTT obtained by Rosenbluth separation. Since it 
is impossible to extract the on-mass shell neutron structure functions directly by 
measuring them using a neutron target, a model based on the deuteron structure 
function is used. The deuteron fit, described in Ref. |102j . relies on a fit of the 
ratio Rp = (Tl/'^t and the assumption Rp = The fit includes photoproduction 
and low data points and data from several other experiments. Since there are 
not enough data at low and no Rosenbluth separation has been done only the 
transverse portion of the cross section is fitted. Fermi motion is taken into account 
in a Plane Wave Impulse Approximation (PWIA). The region (dip) between the 
quasi-elastic peak and the A (1232) resonance is systematically under- predicted 
at low Q^, possibly because Meson Exchange Currents (MEC) and Final State 
Interactions (FSI) have been ignored. An additional empirical function was used to 
fill in the missing strength in the dip region. R£, for deuteron is evaluated by doing 
Fermi-smearing to proton ai and ctt described in Ref. |101] . The Fermi-motion of 
the nucleons in the deuterium is taken into account using a PWIA calculation and 
the Paris |1U3] deuterium wave function. After performing Fermi smearing of 
and F^' the neutron Ff is calculated as 



F^^ = 2F^ - F[ 



(3.95) 



where Ff^ is deuterium structure function per nucleon. 
the nuclei (A,Z) is calculated by the following formula, 



125 

The structure functions of 



N 



(3.96) 



^2 



F^{l + Rp)/{l + u'/Q^)/u. 



The structure functions given in the Eq. 3.96 are further corrected for the EMC effect. 
The EMC effect has been studied by comparing F2 measured on bound nucleons in 
nucleus A and deuterium {x = 0.0085 - 0.09) [m\, {x > 0.125) [TU5] . 



EMC{x, A) = {a^/a^)is = C(x)A'^(^) (3.97) 

where logC(x) = a + 61og(a;) + c(logx)^ and a{x) is an eight order polynomial 
function of a; = Q'^/{2Mv). The "is" subscript means the isoscalarity correction 
(taking into account the fact that for all nuclei). For x lower than 0.0085 the 

EMC correction is taken to be that for x = 0.0085 (no reliable data exists below this 
value). For x higher than 0.7 the EMC correction is taken to be that for x =0.7, and 
the rest is taken into account when Fermi smearing is done. An empirical fit form is 
added to the dip region between the quasi-elastic peak and the A (1232) resonance. 
The function has the same form as in the deuteron fit. 



3.6 Global Fit and Model Iteration 

In order to extract the Born cross section from experimental data a model cross 



section is necessary. This can be seen from Eq. 3.83 where the extracted cross section 



is proportional to the model cross section. Also, the model cross section is used to 
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calculate radiative corrections (Eq. 3.83), Coulomb corrections (Eq. 3.90) and apply 



bin centering corrections (Eq. 3.81). Since the model cross section is used to apply 
all the corrections mentioned above, the extracted cross section clearly depends on 
knowledge of the model cross section. In order to minimize the model dependence, 
an iteration procedure is used. The procedure is described below, step-by-step. 

1. The Born cross section is extracted using the model cross section described in 
Sec. 13.51 The same model is also used to calculate radiative corrections and bin 
centering correction. 

2. The extracted cross section is combined with the cross section results of other 
experiments |106[ 1107] . some of them done dating from 1970. These data are 
fitted with the model cross section described in Sec. 13.51 with some additional 
corrections. These corrections are listed below. The quasi-elastic structure 
function Ff'^ is multiplied by polynomial function of fourth degree in y. The 
parameters of the polynomial are constrained in such way that the value of the 
polynomial function is always positive in the kinematic range of quasi-elastic 
scattering. The inelastic structure function is multiplied by a function of 
W"^ and |108j . When — )• oo the value of this function goes to 1 and does 
not have any impact on A linear A (nuclear) dependence is introduced to 



Ra- The result of the first iteration is shown in Fig. 50 



3. The corrected model is used to calculate radiative corrections again and redo 
the cross section analysis. Steps 1 and 2 are repeated until the fitted model 
cross section converges. 

4. Using the corrected model the Ra — Rd is extracted (by means of a Rosenbluth 
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Figure 50: The agreement between experimental data and model cross section is 
shown for six bins after doing the first iteration of the global fit. 
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separation) for each W"^ bin. Since the dependence of Ra — -R_d is very weak 
two bins, 0.5—2.5 GeV^ and 2.5—4.5 GeV^, are used to parametrize the 
Ra — Rd versus W"^. The parametrizations of Ra — Rd are used in the global 
fit program to run steps 1,2,3,4 repetitively until no improvement in model cross 
section can be found. 



3.7 Extraction of R and F2 

The electron-nuclear inclusive scattering cross section can be written in terms of 
photoabsorption cross sections for transverse (helicity ± 1) photons and longitudinal 



(helicity 0) photons, see Eq. 3.98 



1 (fa 



In Eq. 3.98, where F is the flux of transverse virtual photons, the right side is a 
function of W"^, and e (beam energy). In order to determine ctt and ctl from this 
equation, measurements are done at the same {W^, Q"^) but different e (beam energy). 
After calculating ax and the structure functions R and F2 can be determined by 
the following formulas: 



R=^ and (3.99) 



F2 = [OT + ^l)t^^^, (3.100) 



where K is called "equivalent photon energy". In this experiment the Rosenbluth 
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Figure 51: An example of Rosenbluth separation for carbon target. The slope of the 
curve is Ra — Rd, the intercept is a^/a^. 



separation is done using the nuclear cross section ratio to the deuterium cross sec- 
tion ^80j. The ratio of the cross sections is given by 



^ = ^[l + e'{RA-RD)], (3.101) 



where Rd is the deuterium R and e' = e/(l + eRi)). This method has it's advantage 



compared to Rosenbluth separation method which uses Eq. 3.98 If Rosenbluth sep 



aration is done using Eq. 3.98 systematic uncertainties may have big impact. Since 
eRo is small, it's impact on e' is small, for a 20% uncertainty when R^, = 0.2, the 
uncertainty on e' is only 4%. Taking the ratios of cross sections aA/cyo cancels most 
systematic uncertainties: acceptance, beam charge, offsets in beam energy, HMS an- 
gle, HMS momentum, model dependence to some extent, radiative corrections. 
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An example of Rosenbluth separation is shown in Fig. 51 Cross sections are 
not measured at exactly same W^,Q^ values, instead, for each there are several 
Q'^s, which can be combined to have one average Q^ean value. In this analysis the 
values of Q^s that are combined to perform Rosenbluth separation are chosen to be 
within a AQ^ =0.5 GeV^ range {Qmean is calculated as the mean of these Q^s). After 
calculating Qfneany cross sections are centered to the same value of W'^, Q^ean using 
the model cross section, via 



^Qmean) = (^exp[W ,Q ) . (3.102) 



If the model cross section correction is greater than 50% that particular point is 
excluded from the fit. 

In order to have enough range in e, each Rosenbluth separation is only done if the 
Ae > 0.3 condition is satisfied. Each data point is used only once to avoid correlated 
uncertainties. For this experiment 280 Rosenbluth separations are done for carbon 
and aluminum, 120 for iron, 110 for copper. A program is written to search the cross 
section data for Rosenbluth separation candidate cross section values and perform 
the Rosenbluth separation. The program is designed in such a way that adding data 
from another experiment to the data set of this experiment requires no modification. 
An algorithm was developed to exclude a given cross section value being used more 
than once, and to find data within AQ^ =0.5 range while at same time have the 
minimum spread around the calculated mean Qmean value. This minimizes the model 



dependence introduced by Eq. 3.102 



Though R can be extracted in a model dependent way by using the expression 
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(Pa 2MxF2 1 + eR 

(^Mott^^ — , , „ , (3.103) 



dndE' Q^e \1 + R 

the experimental cross section and a model of the structure function F2, it was not 
done for this analysis. 

3.8 Systematic Uncertainties 

The total systematic uncertainty in the cross section extraction is taken as the sum 
in quadrature of all systematic uncertainties of the quantities that contribute to the 
cross section. 

Systematic uncertainty can be divided into two groups: point-to-point uncer- 
tainties and normalization uncertainties. Point-to-point uncertainties are caused by 
changes in experimental conditions during data acquisition, and therefore their effect 
is uncorrelated between different data points. These include uncertainties arising from 
a variation in the efficiencies of detectors and data acquisition systems, changes from 
one spectrometer setting to another. Point-to-point uncertainties can be removed 
if the same measurement is done more than once, while in contrast, normalization 
uncertainties can't be avoided by doing more measurements. An example of normal- 
ization uncertainty can be a systematic shift in the offset of a current measurement 
device, target thickness measurement, acceptance, etc. 

Kinematic uncertainties are caused by uncertainty in beam energy, spectrometer 
momentum, and spectrometer angle. These uncertainties are estimated from a study 
of elastic electron-proton scattering data. A procedure described in Ref. [76] provides 
a way to estimate the uncertainties in beam energy, spectrometer momentum, and 
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Section 



3.1 



3.4.6 



3.4.6 



3.4.9 



3.4.11 



Quantity 



Beam Energy 
Beam Charge 
Scattered e' Energy 
Scattered e' Angle 
Elect. Dead Time Correction 
Comp. Dead Time Correction 
Acceptance 
Model Dependence 
Radiative Correction 



Uncertainty 



0.05% 
0.3 /iA 
0.06% 
0.2 mrad 
27.0% 

0.2% 

0.6% 

0.6% 

1.0% 



0.25% 
0.37-0.75% 
0.025% 

0.3% 
0.04% 

0.2% 

0.6% 

0.6% 

1.0% 



Table 8: Normalization systematic uncertainties in the experimental parameters (col- 
umn 3) and the corresponding systematic uncertainties in the differential cross section 
(column 4). 

spectrometer angle, and their values are shown in the third column of Table [8j A 
model cross section is used to estimate their impact on measured cross section and is 
shown in the fourth column of Table HI 



The beam charge measurement is discussed in Sec. 3.1 The point-to-point un- 



certainty is estimated to be 0.3%. This is obtained by studying the residuals of the 
measured currents during the calibration procedure. An additional scale uncertainty 
of 0.3% is assumed for the charge measured, due to the UNSER calibration. This is 
estimated by examining runs taken with the same kinematics settings but different 
beam currents on a carbon target ( with 10 fiA steps starting from 20 fiA up to 100 
fiA). 

The scale uncertainty of the HMS acceptance correction is 0.6%. This is estimated 
by changing positions of the target, collimator, magnets and then calculating the 
acceptance. Cross sections are extracted using this shifted acceptance. The spread 
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Section 



3.1 



3.4.6 



3.4.6 



3.4.7 



3.4.8 



3.4.2 



3.4.2 



3.4.4 



3.4.9 



3.4.1 



Quantity 



Beam Energy 
Beam Charge 
Scattered e' Energy 
Scattered e' Angle 
Elect. Dead Time Correction 
Comp. Dead Time Correction 
Trigger Efficiency 
Tracking Efficiency 
Cerenkov Efficiency 
Calorimeter Efficiency 
Charge Symmetric Background 
Acceptance 
Radiative Correction 



Uncertainty 




0.05% 


0.25% 


0.3 /iA 


0.3% 


0.06% 


0.025% 


~ 0.2 mrad 


0.3% 


27.0% 


0.008% 


0.2% 


0.2% 


0.007% 


0.007% 


0.15% 


0.15% 


0.15% 


0.15% 


0.05-0.2% 


0.05-0.2% 


0.1-0.4% 


0.1-0.4% 


0.7% 


0.7% 


3.0% of QE 


< 1.0% 



Table 9: Point-to-point systematic uncertainties in the experimental parameters (col- 
umn 3) and the corresponding systematic uncertainties in the differential cross section 
(column 4). 

of distribution of the ratio of shifted cross sections to cross section with a known 
positions of the target, collimator, magnets is assigned as the systematic uncertainty 



for the acceptance correction. The momentum offset correction discussed in Sec. 3.4.9 
contributes about 0.7% point-to-point uncertainty. 

No normalization uncertainty is assigned to the tracking efficiency. At a fixed 
rate the tracking efficiency should scale linearly with the DC time window. This was 
studied by varying the DC TDC time window and plotting the tracking efficiency 
versus time window width. No nonlinearity is observed according to Ref. [77] . Also 
a point-to-point uncertainty of 0.15% is assigned to the tracking efficiency based on 



the spread of the points in tracking efficiency versus SCIN rate plot, as in Fig. 41 
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The trigger efficiency is better than 99.9%, and no scale uncertainty is assigned 



to it. The point-to-point uncertainty is 0.007%, see Fig. 39, and is neghgible. 



Electronic and computer dead times are discussed in Sec. |3.4.6] A scale uncertainty 
of 0.04% (20 kHz average rate) is found for the electronic dead time, mainly from the 
deviation of the measured value of r (from the plot of electronic dead time versus 
pretrigger rate, electronic dead time is found to be ~ 80 ns ) from the expected value 
of 60 ns. The point-to-point uncertainty (spread of electronic dead time from expected 
value defined by 1 — 60ns x rate) is 0.008% and also is negligible. The uncertainty in 
computer dead time is estimated by taking runs at the same kinematics but different 
prescale factors. Results are found to agree within 0.2%, see Ref. [75]. A point-to- 
point uncertainty of 0.2% is assigned. 

Pion contamination can be up to 3% at some kinematics. It should be the same 
for positron runs, so no explicit subtraction is done. It is automatically subtracted 
when the charge symmetric background is subtracted. The uncertainty is absorbed 
into the charge symmetric background uncertainty. 



The charge symmetric background is discussed in Sec. 3.4.4 At scattered electron 
angles smaller than 50° (low E') point-to-point uncertainty is estimated to be 0.1%, 
for angles greater than 50° (high E') the point-to-point uncertainty is estimated to 



be 0.4%, (see Fig. [36]). 

The effect of the model on the bin centering corrections is studied by varying the 
shape of the model. This is done by supplying an artificial dependence as input 
to the individual DIS and QE cross sections. The point-to-point of 0.5% uncertainty 
due to model dependence is assigned as a result of the study. 

In the kinematic range of this experiment radiative correction are maximum 30%. 
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The systematic error for higher order a contributions are set to be 3% of the size 
of the quasi-elastic tail, relative to the inelastic [H2]. The uncertainty on the cross 
sections due to the radiative correction is estimated at 1% for e dependence, and 
the normalized uncertainty is also according to the radiative correction studies 

done at SLAC [50]. 
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Chapter 4 



Results and Discussions 



In the last chapters we have discussed the cahbration of the detectors and the correc- 
tions apphed to the experimental data to extract the cross sections. In this chapter we 
will present the results of the analysis. First, the extracted differential cross sections 
for the aluminum target will be presented. Next, the results of nuclear dependence 
oi R = ai/ cft will be shown and nuclear dependence of Ra — Rd and Ra — Rc will 
be discussed. The Rosenbluth separated structure functions R — oi^juT and Fi, 
and Ft will be compared with a model obtained from an empirical fit performed in 
the nuclear resonance region. Next, the extracted F2 structure function using model 
cross sections will be discussed in terms of quark-hadron duality. Finally, Rosenbluth 
separated structure functions will be used to study quark-hadron duality in nuclei. 
Before beginning this presentation of the results I will take a short detour to discuss 
how elastic electron-proton scattering is used to gauge our estimate of systematic 
errors. 

4.1 Elastic Electron-Proton Cross Section 

Elastic electron-proton scattering can be used as a cross check of the systematic un- 
certainties present in the setup of an experiment. Since the elastic electron- proton 
cross section has been measured in the past with great precision, any deviation from 
that measurement would indicate a problem present in the collected data or in cali- 
brations. In addition, clastic electron-proton scattering can be used to determine if 
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normalization uncertainties like charge, spectrometer acceptance and detector ineffi- 
ciencies and dead times are taken into account correctly. 
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Figure 52: Ratio of extracted elastic cross sections to a fit of the world's elastic 
data [86j. Point-to-point systematic errors are also shown and are the dominant 
uncertainty. The dashed lines indicate ± 2% range. 



In this experiment elastic data are taken at twelve different kinematic settings. 
The elastic data were analyzed using the same technique as the data in this experiment 
with a few exceptions. In addition to the cuts used in resonance cross section calcu- 
lation, a cut 0.8 < < 1.15 GeV^ is applied to avoid the pion electro-production 
region. In order to compensate for radiated elastic events with W"^ > 1.15 GeV^ (low 
W'^ bins radiate to higher W"^ bins) a correction factor is calculated and applied to 
each 6 bin. For each run the elastic cross sections is averaged over 6 bins using an 
elastic cross section model and the results are shown at a given central value in 



Fig. 52 These elastic cross sections agree with the world's elastic cross section data 
within better than 2%. Since the measurement of elastic cross sections are more sen- 
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sitive to kinematic offsets than inelastic resonance cross section, it is expected that 
systematic uncertainties in the resonance region are, in the worst case, the same as 
in the elastic case. For more information about elastic electron-proton cross section 
extraction method see Ref. |lU9j . 



4.2 Differential Cross Sections 



The method of extracting the cross sections was discussed in the Section |3.4[ The 
differential cross sections were extracted for four different targets: carbon, aluminum, 
iron and copper. Carbon and aluminum data were taken in the range 0.0 < W"^ < 4.5 
GeV^ and 0.5 < Q'^ < 4.5 GeV^. For iron data were taken in the range 0.0 < W'^ < 
4.5 GeV^ and 0.5 < < 2.5 GeV^, and for copper 0.0 < < 4.5 GeV^ and 
2.5 < < 4.5 GeV^. 



In the analysis process, the data were stored in equally spaced W"^ = 0.04 GeV^ 



bins. After iterating the model described in the Section 3.6 , double differential electro- 
production cross sections were extracted. 



Figures 53 to 62 show examples cross section spectra for aluminum (cross sections 



for C, Fe, and Cu are in Appendix 1.3). The blue solid curve in these figures is the 



fit to the extracted cross sections of this experiment. The red curve is the fit based 



on proton and deuteron fits described in the Section 3.6 and was used as the starting 
input cross section model in this analysis. Before starting the iteration procedure, 
77% of the cross section points agreed within 5% of the corresponding model value, 
after the iteration procedure 90% agreed. The iterative procedure is stopped when 
convergence is achieved. Convergence is confirmed by taking the difference with the 
extracted cross sections between successive iterations and making sure that the differ- 
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ence is less than 0.5%. In total three iterations were necessary to achieve convergence. 
Only the statistical uncertainty of the data is shown in the figures, and it is typically 
smaller than the symbol size. 
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Figure 53: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 54: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 55: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 56: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 57: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 58: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 59: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 60: Extracted differential cross section for aluminum compared to the model 
cross section. 
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Figure 61: Extracted differential cross section for aluminum compared to the model 
cross section. 



148 



150 
100 
50 

^ 

> 100 
CD 



50- 



40 



: E = 5.15GeV,e = 30 
: = 2.3-4.1 GeV^ 












































r r 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 r-v 


Old model 

" New model 

- • Data 






















- E = 5.15GeV,e = 33 
= 2.8-4.2 GeV^ 


















r 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 r-v 


- E = 5.15GeV,e = 38 
= 3.2-4.2 GeV^ 















































(GeV^) 



Figure 62: Extracted differential cross section for aluminum compared to the model 
cross section. 
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4.3 Cross Section Ratios and Extraction of Ra — Rd 



During the running period of this experiment the experiment E06-009 |69] took data 
on a deuterium target at the same kinematic points as the current experiment. In 
order to study the nuclear dependence of i?, the deuteron cross sections were also 
extracted by the author in addition to the cross sections from heavier nuclei of this 
experiment. 



The ratio of cross sections a a /err) are shown in Figs. 63 to Fig. 67 The cross sec- 
tion ratios of nuclear targets to deuterium are indicated by different colors. The three 
dip regions (where enough W'^ range exists) correspond to the nucleon resonances P33 
(1^2 = 1.52 GeV^), 5ii {W^ = 2.28 GeV^), and F15 {W^ = 2.82 GeV^). Since the 
Fermi momentum of deuteron is the smallest of all nuclei, the resonance structure is 
not washed out, and the cross sections for deuteron at the positions of resonances are 
larger than that of other nuclei. Therefore, resonances are seen as dips in the cross 
section ratios. Coulomb corrections were applied to the cross section ratios shown in 



these figures as discussed in Section 3.4.12 and the strength of the correction is shown 



in Fig. |49j Since these are ratios of cross sections and not of separated structure func- 
tions, therefore it is necessary to assume that there is no nuclear dependence of R in 



order to interpret these data. This follows from Eq. 4.104 below. 



The nuclear to deuteron cross section ratio can be written in terms of structure 
function R as shown in the following equation 



0^4 



(J- 



l + eRA 
l + eRo 



a. 



D 



[1 + e' {Ra - Rd)] 



(4.104) 



where e' = t/ {l + tRo)- At fixed values of W'^ and the right side of this equation is 
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Figure 63: Cross section ratios of nuclear targets to deuteron versus is shown. 
Colors indicate different nuclear targets. The three dip regions (where enough 
range exists) correspond to the nucleon resonances P33 {W"^ = 1.52 GeV^), 5*11 
(iy2 ^ 2.28 GeV^), F15 {VV^ = 2.82 GeV^) and are indicated by arrows. Since the 
Fermi momentum of deuteron is the smallest of all nuclei, the resonance structure is 
not washed out for it, and the cross sections for deuteron at the positions of reso- 
nances are larger than that of other nuclei. Therefore, resonances are seen as dips in 
the cross section ratios. 
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Figure 64: As in Fig. 
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Figure 65: As in Fig. 
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Figure 66: As in Fig. 
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Figure 67: As in Fig. 
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function of e' only and represents the equation of a straight hne with constant / 
and slope {a^/a^){RA — Rd)- 

Since eRo is small, the uncertainty of e' due to Rd is negligible. The Ra — Rd 
results are consistent with zero for W"^ > 2 GeV^, in agreement with DIS data where 
no significant A dependence of R was found. These results indicate that contributions 
to R from nuclear effects are small. The range of this data is not low enough to 
investigate the predictions of Miller pT] who suggests a significant enhancement of 
longitudinal cross section at x = 0.4 and = 0.3 GeV^ for heavy nuclei. This will be 
possible when the results of E02-109 [18] experiment are available. Combining these 
data with the data of E02-109 will allow to perform a precise Rosenbluth separation 
in wide range of for carbon, aluminum and iron. 



The observed enhancement of Ra — Rd near W =1.5 GeV (Figs. 68 and 69) 
indicates some nuclear dependence, which vanishes with increase of and with 
nuclear number A. A final conclusion about the origins of the enhancement of Ra—Rd 
can not be drawn until the analysis of experiment E06-009 [69] is complete and 
precision deuteron data become available which will allow a better representation of 



(T£) in Eq. |4.104[ These data will allow improvements to the deuteron model and 
precise cross sections in the same kinematic region as the current experiment. 



In addition to extracting Ra — Rd, Ra — Rc was also extracted, shown in Fig. 70 



The model for Rc was obtained from the current fit to the data as described in 
Section 3^ No nuclear dependence of Ra — Rc was found over the range of W"^ and 
in the four bins for aluminum, iron and copper. 
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Figure 68: The results for Ra — Rd are plotted versus W'^ for two Q"^ bins. Top two 
plots are for carbon, and bottom two are for aluminum. 
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Figure 69: The results for Ra — Rd are plotted versus W"^ for two bins. Top plot 
is for iron, and bottom one is for copper. 

4.4 Rosenbluth Separated R 

Once the cross section model has converged, the final cross section can be extracted 
and used to perform Rosenbluth separation using the following formula 



dndE' 



r = ((Tr + eo-i) 



(4.105) 



where T is the flux of virtual photons and is given by Eq. 1.13 In order to perform a 



Rosenbluth separation using the Eq. |4.105| it is necessary to measure cross sections at 
the same values of and but at different values of e. In this analysis the cross 
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Figure 70: The results for Ra — Rc are plotted versus W"^ for two Q'^ bins. Top two 
plots are for aluminum, then iron and copper. 
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sections are calculated at fixed W'^ bins, as described in Section 3.4.9, where for each 
W"^ there are several values of corresponding to different e's. In order to reduce 
systematic uncertainties of R, each Rosenbluth separation was performed when there 
were more then three points with an e range larger than 0.3, one of the e points must 
be greater or smaller than 0.5 and all points were within AQ^ < 0.5 GeV^ range. The 
average e range for all Rosenbluth separations is 0.5±0.1. The model cross section 
obtained after the iteration procedure was used to interpolate all cross sections within 
the AQ^ < 0.5 range into an average value of Q^. In order to avoid having correlated 
errors no data point was used twice in the Rosenbluth separation. The value of per 
degree of freedom is about 1 indicating that point-to-point systematic uncertainties 



are taken into account correctly. A few Rosenbluth separations are shown in Fig. 71 



Rosenbluth separations were performed for carbon, aluminum iron and copper. 
The Coulomb field of these nuclei are not negligible and a correction is necessary to 
take into account its effect. This is very important for Rosenbluth separation since 
the Coulomb correction is e (beam energy) dependent and can change the value of R 
dramatically. 

Another source of systematic uncertainty is the e dependence of radiative correc- 
tions. Although the radiative corrections for inclusive electron-nuclei are well under- 
stood, at the level of ~ 1%, it is still possible to have a systematic error in R if the 
correction is large. In this experiment the maximum radiative correction was about 
30%. In order to study the effect of radiative corrections, an artificial e dependence is 
introduced to cross sections during the Rosenbluth separation and then the separated 
Rs are compared to the Rs separated without the artificial e dependence. The size 
of this artificial e dependence is set to be within ± 1% range. A systematic uncer- 
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5.15 18.0 1.91 0.91 

3.27 33.0 1.98 0.74 

2.10 75.0 2.06 0.30 

2.10 75.0 2.06 0.30 
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Figure 71: Some examples of Rosenbluth separated R for aluminum. The horizontal 
axis, e, is the relative polarization of longitudinal photons. The vertical axis is the 
reduced cross section defined in the left side of Eq. 4.105[ The intercept of the linear 
fit with the y axis defines ctt, while the slope defines ai- R is the ratio cri/aT- For 
each Rosenbluth separation the beam energy, electron scattering angle, and e are 
given in the legend. 



161 



Q' = 0.9 (GeV) 



O 0.5 

DC 



O 0.5 

DC 




— 



= 3.0 (GeV") 



O 0.5 

DC 



Q2 = 3.7 (GeV^) 



O 0.5 

DC 




(GeV^) 



Figure 72: Rosenbluth separated R of carbon as a function of PF^ at four different 
ranges. In order to take into account the dependence of R for each range, the 
model obtained from the global fit is used to move Rs to an average shown in the 
plot. 
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Figure 73: Rosenbluth separated R of aluminum as a function of at four different 
ranges. In order to take into account the dependence of R for each range, the 
model obtained from the global fit is used to move Rs to an average shown in the 
plot. 
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Figure 74: Rosenbluth separated R of iron and copper as a function of at four 
different ranges. In order to take into account the dependence of R for each 
range, the model obtained from the global fit is used to move Rs to an average 
shown in the plot. 
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tainty of AR = 0.028 is found and is assigned to all values of R attributed to the e 
dependence of the radiative corrections. 

The Rosenbluth separated R for carbon, aluminum, iron and copper are shown 



as a function of W"^ over the resonance region in Figs. 72, 73 and 74 respectively. 
Each figure has four plots indicating four ranges from top to bottom respectively: 
0.5-1.0 GeV^, 1.0-2.5 GeV^, 2.5-3.5 GeV^ and 3.5-4.5 GeV^. For each range 
the R is interpolated, using models of R, to the average value shown on top of each 
plot. The lines shown on all plots correspond to different models for R. The red line 
(Resonance new) indicates the model obtained after iterating the cross section model. 
The blue line is the initial model used to extract the cross sections. The green line is 
the DIS R fit "R1998" [25] and is based on a fit to world R measurements. In general 
all three models agree with the data, except for the lowest case. At lowest = 
0.9 GeV^ the agreement is better with the "Resonance new" model, which is natural, 
since this model is based on the first cross section measurements done in the nucleon 
resonance region, but it still needs improvement. This improvement will be done 
when the low data of E02-109 [12] is available to include in the iterative procedure 



discussed in Sec. 3^ At = 2.0 GeV^ for carbon and aluminum the agreement 
between data and DIS model for W"^ < 3 is lacking. This, again, may be caused by not 
having the low data of the E02-109 in the iterative procedure. For iron at = 
2.0 GeV^ the data have larger uncertainties and it is impossible to make a definite 
observation about the agreement of data and models. This is further complicated by 
the fact that the difference between the different models is not significant. At = 
3.0 and 3.7 GeV^ for carbon, aluminum and copper the agreement of the data and 
models is good within the statistical errors and the difference between the models 



165 



decreases with increasing Q^. 

4.5 F2 Structure Function and Duality in Nuclei 

The extraction of the F2 structure function from cross section data can only be 
accomphshed with some input for the ratio R of the longitudinal to transverse cross 
sections. In order to perform a model-independent extraction of the unpolarized 
structure function F2 from inclusive cross section data, R is obtained by performing a 
Rosenbluth separation, described earlier. The Rosenbluth separated Ra still contains 
some model dependence since the Ra = Rd + Ffu where Rd is known from a model 
and Ffit is found by fitting Ra — -Rd at two bins as shown in Figures 68 and 69 



However, it should be noted that the model of -R^ is based on a fit of Rosenbluth 
separated data [75] . 

After extracting Ra one can calculate using 



Mott i + en 



The results of extracted F^ are shown in Figs. 75 , 76 , 77 and 78 for carbon, aluminum, 
iron and copper respectively. In the top plot the F^ structure function is shown versus 
the Bjorken x variable, in the second plot, it is shown versus the Nachtmann variable 



^ defined in Eq. 1.35 The curves shown on the bottom plots of Figs. 75, 76, 77 



and [78] are the proton F2 based on different parametrizations and are corrected for 
the nuclear EMC effect. The parametrization ALLM97 |110] is F2 of the proton in 
the DIS region. It rests on a Regge motivated approach, similar to that used earlier 
by Donnachie and Landshoff |111] . extended into a large regime in a way that 
is compatible with QCD expectations. The parametrization of SLAG data is based 
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on the old SLAC structure functions and authored by Bodek and Atwood. The 
parametrization GJR08 |112] stems from the dynamical parton distributions of the 
nucleon. 



In the top plot of Fig. 75 (carbon) and Fig. 76 (aluminum), scaling in Bjorken x 
is not observed for the most part. This is expected since the F2 set includes values 
at momentum transfers ranging from Q^=0.5 GeV^ at low x to Q^=A.5 GeV^ at the 
higher x values. In addition the quasi-elastic peak and resonances are not completely 
washed out by the Fermi motion of the nucleons inside the nucleus. The scaling 
is violated the most at x = 1, which corresponds to the quasi-elastic peak. For 
aluminum and carbon the Fermi momentum is lower than for heavier targets and 
the quasi-elastic peak is not broadened enough to show scaling. In the bottom plot 



of Fig. 75 and Fig. 76 F2 is plotted versus the Nachtmann variable ^, which takes 
account the finite target mass, and seems to become independent of Q^. Note that 



the F2 of iron, shown in Fig. 77, has a range of 0.5-3.0 GeV^ in contrast to both 
carbon and aluminum where the is 0.5-4.5 GeV^. Even though Fermi momentum of 
iron is larger than that of carbon and aluminum the F2 scaling in Bjorken x is still 
not observed (top plot). In the bottom plot, where F2 is plotted versus Nachtmann C,, 
the situation is completely different and scaling is observed. The quasi-elastic peak is 
not visible. This is an important observation, since it shows that Fermi momentum 
of iron is enough to remove dependence of F2 structure function near the quasi- 



elastic peak even at values around 1 GeV^. In Fig. 78 F2 is shown for copper 
in the range between 1.3-4.3 GeV^ higher than that for other targets. As it can 
be seen from the figure that F2 scales versus both Bjorken x and Nachtmann ^. For 
copper it can be concluded that smearing effect of Fermi motion on the resonances 
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combined with high enough (> 1.5 GeV^), scahng of is achieved even without 
applying target mass corrections. 

For nuclei heavier than iron even the quasi-elastic peak is washed out by the 
Fermi smearing at higher Q^, and scaling is seen at all values of x and ^. Here 
the resonance region is essentially indistinguishable from the DIS scaling regime. The 



observed scaling is consistent with the duality arguments discussed in Sec. |1.37[ There 
local averaging of the structure function over ^ for the nucleon resonances seen at low 

and quasi-elastic peak, produces structure function consistent with the high Q'^ 
scaling limit of structure function. Here we see indications that the Fermi motion of 
the nucleons in the nucleus are performing local averaging and there is no need to 
use the finite energy sum rule in order to quantify the similarity of scaling functions 
in resonance and deep inelastic regimes. 

Qualitatively, the nuclear effects in the resonance region appear to be similar to 
those in the deep inelastic region. This is surprising since the nuclear dependence of 
the scaling structure functions is not expected to be the same as the nuclear depen- 
dence of resonance production. There is a priori no reason why these modifications 
would be the same as those for structure functions measured in deep inelastic scatter- 
ing. On the other hand, this may be viewed as another consequence of quark-hadron 
duality. 

The extracted F2 structure functions (per nucleon) are similar for all nuclear 



targets as it can be seen in Fig. 79 Figure 79 shows the structure function for 
carbon, aluminum, iron and copper at ^ =0.3, 0.4, 0.5, 0.6, 0.7 and 0.8. The dashed 
lines indicate fixed values of Nachtmann ^. The kinematic coverage for the iron and 
copper targets is less than for carbon and aluminum. For iron ranges from 0.5-2.5 
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Figure 75: Top plot: Structure function per nucleon versus Bjorken x for carbon 
(per nucleon) from the present measurement. Bottom plot: The i>W2 = F2 structure 
function for carbon (per nucleon) as a function of Nachtmann ^. The ranges are 
given at each 6 angle. The curves are the proton F2 parametrizations at = 10 
GeV^, corrected for the nuclear EMC effect. Arrows on the right side (^ ~ 0.8) of 
the bottom plot indicate positions of quasi-elastic peaks. 
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Figure 76: Top plot: Structure function per nucleon versus Bjorken x for aluminum 
(per nucleon) from the present measurement. Bottom plot: The = F2 structure 
function for aluminum (per nucleon) as a function of Nachtmann ^. The ranges 
are given at each 9 angle. The curves are the proton F2 parametrizations at = 
10 GeV^, corrected for the nuclear EMC effect. Arrows on the right side (^ ~ 0.8) of 
the bottom plot indicate positions of quasi-elastic peaks. 
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Figure 77: Top plot: Structure function per nucleon versus Bjorken x for iron (per 
nucleon) from the present measurement. Bottom plot: The VW2 = F2 structure 
function for iron (per nucleon) as a function of Nachtmann ^. The ranges are 
given at each 6 angle. The curves are the proton F2 parametrizations at = 10 
GeV^, corrected for the nuclear EMC effect. Arrows on the right side (^ ~ 0.8) of 
the bottom plot indicate positions of quasi-elastic peak. 
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Figure 78: Top plot: Structure function per nucleon versus Bjorken x for copper 
(per nucleon) from the present measurement. Bottom plot: The VW2 = F2 structure 
function for copper (per nucleon) as a function of Nachtmann ^. The ranges are 
given at each 6 angle. The curves are the proton F2 parametrizations at = 10 
GeV^, corrected for the nuclear EMC effect. The arrow on the right side (^ ~ 0.85) 
of the bottom plot indicate the position of quasi-elastic peak. 
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Figure 79: Structure function per nucleon for C, Al, Fe and Cu as a function of 
at fixed values of Nachtmann variable ^. The dashed line serves only to guide the 
eye. 

GeV^, for copper from 2.5-4.5 GeV^. At values of ^ corresponding to the top of the 
quasi-elastic peak, the structure function decreases slightly with A, as the increased 
Fermi momentum broadens and lowers the peak. At lower values of C,, far away from 
the quasi-elastic peak, the structure function per nucleon is nearly identical for all of 
the nuclei. 

In Fig. world data on the structure function F2 of the proton versus is shown 



at several values of x. In Fig. 79 a similar plot is depicted for nuclear structure 
function in a limited range. Overall the two plots follow the same pattern, at small 
^ (x) F2 (-Ff ) ^ises while at high ^ (x) decreases. As it is described in Sec. 



1.1.1 



this 



behavior is consistent with the QCD predictions indicating scaling violations. In the 
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F2 structure function at higher ^ and in a much smaller range of (compared to 
F2 of proton) observed scaling violation is significant and can not be considered as 
logarithmic in Q^. Higher values of ^ are near the quasi-elastic peak meaning higher 
twist contributions are large and the observed violation should have l/((5^)" (power 
correction) dependence. 

4.6 Structure Functions and Duality Studies 

The Rosenbluth separation allows the extraction of cxi and (Tt independently. Then 
the structure functions F2, Fi, and F^ can be calculated using the knowledge of ai 
and (Tt by the following formulas 

K 

Fi = -^Mar (4.107) 



' l + ^][aT + (JL] (4.108) 



where is a kinematic variable. From the definition of Fi one can see that it is 
related only to the transverse virtual photon coupling, while F2 is a combination of 
both transverse and longitudinal couplings. It is useful therefore to define a purely 
longitudinal structure function Fl 



1 + ^ ) F2 - 2xFi. (4.109) 



One of the goals of this experiment was to study quark-hadron duality for the 
transverse structure function Fi, the longitudinal structure function and for F2. 



In Section 4.4, where R is extracted in the resonance region and compared to R in 
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DIS region, no significant difference was found in these regimes, except at the lowest 
(0.9 GeV^). The similarity of R in the two kinematically distinct regions for high 
enough is expected from quark-hadron duality and for the proton has been studied 
extensively. The results of this experiment allow a similar comparison to be made 
in the resonance region for heavy nuclei. Studying quark-hadron duality in nuclear 
targets can provide information on what role the nuclear effects have. 

The results of structure functions Fi, and are presented below and compared 
to models of the structure functions in DIS region. The 2xF\ and are plotted as 
a function of x at different to provide a direct comparisons of the resonance 
transverse and longitudinal structure functions to those in the DIS region, and are 



seen in Figs. [80|-|88| In all plots the resonance structure is not visible because of 
averaging effect of Fermi motion of nucleons inside nuclei. The red curve shown in 
the plots is calculated from the fit to current model, the blue curve is calculated 
from MRST (NNLO) fit |113] . without target mass corrections. The green curve is 
calculated from the MRST (NNLO) fit |113] with target mass corrections. In MRST 
(NNLO) fit [1131 the parton distributions are determined by a global analysis of a wide 
range of deep inelastic and related hard scattering data. The Bjorken x dependence 
of the distributions are parametrized at some low momentum scale, and a fixed order 
(either LO or NLO or NNLO) DGLAP [331 [311 135] evolution performed to specify 
the distributions at the higher momentum scales where data exist. The target mass 
correction is done (calculated based on prescription of Georgi and Pohtzer [3H]) using 
a program provided by M.E. Christy [III]. The MRST and MRST (NNLO) are also 
corrected for EMC effect, and isospin weighting (see below) and nucleon structure 
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functions from Ref. |115] . The isospin correction is given as 

= {Z + {A- Z){F^/F^)) j A (4.110) 

and multiplies the structure functions of the proton. 

The comparison of the DIS data at high to resonance data at low values 



follows the original idea of Bloom and Gilman, discussed in Sec. 1.3.1 Bloom and 
Gilman compared low resonance data to high DIS data and argued that the 
observed similarity between the two structure functions indicates that resonances are 
not a separate entity but are an intrinsic part of the scaling behavior of vW^iy^Q'^^ 
structure function |?T]. The shortcoming of this approach is that the comparison is 
made at different values, although at the same ^, which is equivalent to making 
the comparison at different Bjorken x and sensitive therefore to different parton dis- 
tributions. In order to be more precise the comparison of DIS data should be made 
at the same where the resonance data is taken. This is achieved by utilizing par- 
ton distributions determined by a global analysis of a wide range of deep inelastic 
and related hard scattering data, such as MRST (NNLO) fit [113J (the green curve). 
Comparison of resonance region data with PDF-based global fits, DGLAP |33l [311 ES] 

evolution and target mass correction allows the resonance-scaling comparison to 
be made at the same values of x and . 

The transverse structure function 2xF\ versus Bjorken x for carbon, aluminum. 



iron and copper is shown in Figs. 80 , 83 and 86 the longitudinal structure function 



Fi^ versus Bjorken x is shown in Figs. 82, 85 and 88 The structure function F2 is 
shown in Figs. 81, 84, 87 The structure functions are compared at four different 



ranges for carbon and aluminum. For iron only the two lowest ranges are shown. 
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Figure 80: The purely transverse structure function 2xFi of carbon per nucleon, 
measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect and isospin weighting. Errors are smaller than the symbol 
size. 
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Figure 81: F2 structure function of carbon per nucleon, measured in the resonance 
region as a function of Bjorken x is compared with MRST(NNLO) with and without 
target mass correction. Both curves are corrected for the nuclear EMC effect and 
isospin weighting. Errors are smaller than the symbol size. 



178 



0.3 
0.2 



Q2 = 0.9 (GeV^) 




Figure 82: The purely longitudinal structure function Fl of carbon per nucleon, 
measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect. 
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Figure 83: The purely transverse structure function 2xFi of carbon per nucleon, 
measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect and isospin weighting. Errors are smaller than the symbol 
size. 
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Figure 84: F2 structure function of aluminum per nucleon, measured in the resonance 
region as a function of Bjorken x is compared with MRST(NNLO) with and without 
target mass correction. Both curves are corrected for the nuclear EMC effect and 
isospin weighting. Errors are smaller than the symbol size. 
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Figure 85: The purely longitudinal structure function Fl of aluminum per nu- 
cleon, measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect. 
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Figure 86: The purely transverse structure function 2xFi of iron and copper per 
nucleon, measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect and isospin weighting. Errors are smaller than the symbol 
size. 
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Figure 87: F2 structure function of iron and copper per nucleon, measured in the 
resonance region as a function of Bjorken x is compared with MRST(NNLO) with 
and without target mass correction. Both curves are corrected for the nuclear EMC 
effect and isospin weighting. Errors are smaller than the symbol size. 
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Figure 88: The purely longitudinal structure function Fl of iron and copper per 
nucleon, measured in the resonance region as a function of Bjorken x is compared with 
MRST(NNLO) with and without target mass correction. Both curves are corrected 
for the nuclear EMC effect. 
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for copper only the two highest since there is no data at high for iron and at low 
data for copper. The agreement of MRST (NNLO) parametrization with the data 
is the poorest as can be seen in all the plots, while MRST (NNLO) with target mass 
corrections does a better job at describing the data. This makes it clear that target 
mass effects are required to describe the data. 

These results confirm the validity of quark-hadron duality in the resonance region 
in the separated transverse Fi and longitudinal structure functions of the nuclei 
for larger than 2 GeV^ and x > 0.6. The duality is similarly observed in the F2 
structure functions for larger than 2 GeV^ and x > 0.6. For larger than 2 GeV^ 
and X < 0.6 duality is violated up to maximum 20% in all structure functions. This 
can be caused by nuclear effects. 

Qualitatively, the nuclear effects in the resonance region appear to be similar 
to those in the deep inelastic region. This similarity may be viewed as another 
consequence of quark-hadron duality. 
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Chapter 5 



Summary and Conclusions 



Inclusive electron-nuclear scattering cross sections were measured in the nucleon res- 
onance region at 0.0 < W"^ < 4.5 GeV^ for four-momentum transfer values between 
0.5 < < 4.5 GeV^. The data were taken in May- July 2007 at Jefferson Labora- 
tory (Newport News, Va, USA) at scattering angles between 12° and 76° degrees on 
i^C, 2^A1, s^Fe, ^'^Cu targets with electron beam energy 2.097, 3.116, 3.269, 4.074, 
4.134 and 5.150 GeV. The cross sections are estimated to have about 2% systematic 
uncertainties. 

An empirical fit to inelastic electron-nuclei scattering has been performed which 
describes available data reasonably well, within 3% to 5%, for 0.0 < W'^ < 4.5 GeV^ 
and 0.5 < < 5.0 GeV^. The fit is useful in the evaluation of radiative corrections 
to experimental data, for extraction of spin structure functions from asymmetry mea- 
surements, and for the evaluation of structure function moments. 

The cross sections extracted from a deuteron target are used to study the nuclear 
dependence of the structure function R. The cross section ratios of nuclear and 
deuteron targets showed no nuclear dependence of R. This was further confirmed after 
extracting Ra — Rd hj performing Rosenbluth separations. No significant nuclear 
dependence was found for > 2 GeV^ for all targets. The observed enhancement of 
Ra — Rd near W'^ = 1.5 GeV^ indicates some nuclear dependence which vanishes with 
increase of and with nuclear number A. A final conclusion about the origins of the 
enhancement of Ra — Rd can only be made after final analysis of experiment E06- 
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009 [Sn] which will provide an iterated model and precise cross sections for deuteron 
in the same kinematic region as the current experiment. The results of Ra — Rc 
indicate that there is no nuclear dependence of R for heavy nuclei. 

The F2 structure function is extracted for all measured cross sections in a model 
dependent way by interpolating Ra — Rd with a function and using a model for 
deuteron Rj:,. The extracted F2 shows that resonance structure is washed out by the 
Fermi motion of nucleons in nuclei for > 1 GeV^ for nearly all targets and the 
scaling regime of F2 versus Nachtmann ^ is reached. Since the Fermi momentum of 
copper is larger than that of other nuclei and copper data are taken at > 1-5 GeV^ 
scaling is observed even versus Bjorken x. These results suggest that Fermi motion 
of nucleons inside nuclei are doing the local averaging of the structure function F2 
over for the nuclear resonances and produces a structure function consistent with 
the high scaling limit of the structure function. 

The structure functions F2, Fi, F^, and R obtained by Rosenbluth separation have 
allowed us to study quark-hadron duality on nuclear targets in both the transverse 
and longitudinal channels. Comparison made to target mass corrected proton MRST 
(NNLO) parametrization |113j with the EMC correction and isospin weighting ap- 
plied, showed a similarity between the structure functions. Qualitatively, the nuclear 
effects in the resonance region appear to be similar to those in the deep inelastic 
region. This may be viewed as another consequence of quark-hadron duality. 
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A APPENDIX 



1.1 Tables of Rosenbluth Separated R, Fi, F2 and 







X 


R 


AR 


Fi 




F2 


AF2 


Fl 


AFl 


1 


110 


2.096 


0.901 


0.353 


0.081 


0.032 


0.001 


0.033 


0.001 


0.020 


0.004 


1 


110 


3.773 


0.943 


0.155 


0.097 


0.009 


0.000 


0.010 


0.001 


0.003 


0.002 


1 


150 


2.079 


0.885 


0.271 


0.071 


0.035 


0.001 


0.034 


0.002 


0.017 


0.004 


1 


150 


3.753 


0.933 


0.170 


0.097 


0.009 


0.001 


0.011 


0.001 


0.003 


0.002 


1 


190 


2.062 


0.869 


0.175 


0.065 


0.040 


0.002 


0.036 


0.002 


0.012 


0.005 


1 


190 


3.734 


0.923 


0.212 


0.100 


0.009 


0.001 


0.012 


0.001 


0.004 


0.002 


1 


230 


2.045 


0.854 


0.329 


0.076 


0.038 


0.002 


0.039 


0.002 


0.022 


0.005 


1 


230 


3.714 


0.914 


0.326 


0.109 


0.010 


0.001 


0.013 


0.001 


0.006 


0.002 


1 


270 


2.028 


0.839 


0.427 


0.082 


0.039 


0.002 


0.042 


0.002 


0.028 


0.005 


1 


270 


3.694 


0.905 


0.261 


0.099 


0.011 


0.001 


0.014 


0.001 


0.005 


0.002 


1 


310 


1.999 


0.823 


0.295 


0.063 


0.045 


0.002 


0.043 


0.002 


0.022 


0.005 


1 


310 


3.704 


0.896 


0.309 


0.088 


0.011 


0.000 


0.015 


0.001 


0.006 


0.002 


1 


350 


1.984 


0.809 


0.192 


0.055 


0.050 


0.002 


0.045 


0.002 


0.016 


0.005 


1 


390 


2.009 


0.798 


0.275 


0.044 


0.049 


0.001 


0.047 


0.001 


0.021 


0.003 


1 


430 


2.016 


0.786 


0.281 


0.049 


0.050 


0.001 


0.048 


0.001 


0.022 


0.003 


1 


430 


3.640 


0.869 


0.326 


0.081 


0.013 


0.001 


0.018 


0.001 


0.007 


0.002 


1 


470 


2.000 


0.772 


0.292 


0.047 


0.053 


0.001 


0.051 


0.001 


0.024 


0.003 


1 


510 


2.003 


0.761 


0.297 


0.044 


0.054 


0.001 


0.052 


0.001 


0.024 


0.003 


1 


510 


3.598 


0.851 


0.285 


0.076 


0.015 


0.001 


0.020 


0.001 


0.007 


0.002 


1 


550 


1.984 


0.748 


0.306 


0.045 


0.056 


0.001 


0.055 


0.001 


0.026 


0.003 


1 


590 


1.952 


0.733 


0.361 


0.055 


0.058 


0.002 


0.059 


0.002 


0.031 


0.004 


1 


590 


3.556 


0.834 


0.217 


0.068 


0.018 


0.001 


0.021 


0.001 


0.006 


0.002 


1 


630 


1.936 


0.721 


0.363 


0.052 


0.062 


0.002 


0.062 


0.002 


0.032 


0.004 


1 


670 


1.921 


0.709 


0.319 


0.052 


0.067 


0.002 


0.065 


0.002 


0.030 


0.004 


1 


670 


3.513 


0.817 


0.351 


0.075 


0.018 


0.001 


0.024 


0.001 


0.011 


0.002 



Table 10: Structure functions of carbon. 



200 







X 


R 


Ai? 


Fi 




F2 


AF2 


Fl 


AFl 


1.710 


1.905 


0.697 


0.325 


0.050 


0.070 


0.002 


0.068 


0.002 


0.032 


0.004 


1.710 


3.466 


0.807 


0.341 


0.074 


0.020 


0.001 


0.026 


0.001 


0.011 


0.002 


1.750 


1.889 


0.685 


0.367 


0.050 


0.072 


0.002 


0.072 


0.002 


0.036 


0.004 


1.750 


3.445 


0.799 


0.278 


0.068 


0.022 


0.001 


0.027 


0.001 


0.010 


0.002 


1.790 


1.804 


0.665 


0.318 


0.046 


0.084 


0.002 


0.079 


0.002 


0.036 


0.004 


1.790 


3.542 


0.796 


0.224 


0.077 


0.022 


0.001 


0.026 


0.001 


0.008 


0.003 


1.830 


1.786 


0.653 


0.299 


0.043 


0.090 


0.002 


0.083 


0.002 


0.035 


0.004 


1.830 


3.522 


0.788 


0.252 


0.077 


0.023 


0.001 


0.028 


0.001 


0.009 


0.003 


1.870 


1.769 


0.641 


0.284 


0.045 


0.095 


0.002 


0.086 


0.002 


0.035 


0.005 


1.870 


3.553 


0.782 


0.264 


0.059 


0.024 


0.001 


0.029 


0.001 


0.010 


0.002 


1.910 


1.825 


0.639 


0.323 


0.045 


0.090 


0.002 


0.085 


0.002 


0.037 


0.004 


1.910 


3.387 


0.767 


0.243 


0.081 


0.028 


0.001 


0.033 


0.001 


0.010 


0.003 


1.910 


3.920 


0.792 


0.292 


0.098 


0.020 


0.001 


0.026 


0.001 


0.009 


0.002 


1.950 


1.803 


0.628 


0.346 


0.038 


0.095 


0.001 


0.091 


0.001 


0.041 


0.003 


1.950 


3.178 


0.748 


0.186 


0.084 


0.034 


0.001 


0.037 


0.002 


0.009 


0.003 


1.950 


3.795 


0.780 


0.340 


0.076 


0.022 


0.001 


0.029 


0.001 


0.011 


0.002 


1.990 


1.785 


0.617 


0.357 


0.037 


0.099 


0.001 


0.095 


0.001 


0.044 


0.003 


1.990 


3.154 


0.740 


0.322 


0.092 


0.034 


0.001 


0.041 


0.002 


0.016 


0.003 


1.990 


3.774 


0.773 


0.101 


0.060 


0.025 


0.001 


0.028 


0.001 


0.004 


0.002 


2.030 


1.700 


0.597 


0.342 


0.042 


0.113 


0.002 


0.104 


0.002 


0.046 


0.005 


2.030 


3.752 


0.766 


0.238 


0.067 


0.025 


0.001 


0.030 


0.001 


0.009 


0.002 


2.070 


1.683 


0.586 


0.300 


0.040 


0.121 


0.002 


0.107 


0.002 


0.043 


0.005 


2.070 


3.107 


0.723 


0.345 


0.089 


0.038 


0.001 


0.046 


0.002 


0.019 


0.004 


2.070 


3.731 


0.758 


0.238 


0.065 


0.026 


0.001 


0.032 


0.001 


0.010 


0.002 


2.110 


1.665 


0.575 


0.380 


0.043 


0.123 


0.002 


0.115 


0.002 


0.054 


0.005 


2.110 


3.084 


0.715 


0.336 


0.088 


0.040 


0.001 


0.048 


0.002 


0.019 


0.004 


2.110 


3.756 


0.753 


0.171 


0.060 


0.028 


0.001 


0.032 


0.001 


0.007 


0.002 


2.150 


1.648 


0.565 


0.319 


0.039 


0.133 


0.002 


0.118 


0.002 


0.048 


0.005 


2.150 


3.060 


0.707 


0.271 


0.083 


0.043 


0.002 


0.049 


0.002 


0.016 


0.004 


2.150 


3.706 


0.745 


0.150 


0.051 


0.030 


0.001 


0.034 


0.001 


0.007 


0.002 



Table 11: Structure functions of carbon. 



201 







X 


R 


Ai? 


Fi 




F2 


AF2 


Fl 


AFl 


2.190 


3.036 


0.699 


0.277 


0.082 


0.045 


0.002 


0.051 


0.002 


0.017 


0.004 


2.190 


3.683 


0.738 


0.228 


0.053 


0.031 


0.001 


0.036 


0.001 


0.010 


0.002 


2.230 


3.347 


0.713 


0.248 


0.049 


0.039 


0.001 


0.046 


0.001 


0.014 


0.002 


2.270 


1.597 


0.535 


0.374 


0.040 


0.149 


0.003 


0.134 


0.002 


0.060 


0.005 


2.270 


2.989 


0.683 


0.250 


0.065 


0.051 


0.001 


0.056 


0.002 


0.017 


0.003 


2.270 


3.648 


0.724 


0.180 


0.059 


0.035 


0.001 


0.039 


0.001 


0.009 


0.003 


2.310 


1.576 


0.524 


0.333 


0.033 


0.159 


0.002 


0.138 


0.002 


0.056 


0.004 


2.310 


2.965 


0.675 


0.233 


0.063 


0.054 


0.001 


0.058 


0.002 


0.017 


0.003 


2.310 


3.627 


0.717 


0.226 


0.061 


0.036 


0.001 


0.042 


0.001 


0.012 


0.003 


2.350 


1.571 


0.517 


0.370 


0.030 


0.162 


0.002 


0.144 


0.002 


0.062 


0.004 


2.350 


2.942 


0.667 


0.315 


0.067 


0.055 


0.001 


0.063 


0.002 


0.023 


0.003 


2.350 


3.650 


0.713 


0.347 


0.079 


0.034 


0.001 


0.044 


0.002 


0.017 


0.003 


2.390 


1.632 


0.520 


0.422 


0.038 


0.152 


0.003 


0.142 


0.002 


0.067 


0.005 


2.390 


2.918 


0.659 


0.241 


0.063 


0.059 


0.001 


0.063 


0.002 


0.019 


0.003 


2.390 


3.627 


0.706 


0.232 


0.071 


0.038 


0.001 


0.045 


0.002 


0.013 


0.003 


2.430 


1.617 


0.511 


0.413 


0.038 


0.160 


0.003 


0.148 


0.003 


0.068 


0.005 


2.430 


2.894 


0.651 


0.277 


0.061 


0.062 


0.001 


0.068 


0.002 


0.022 


0.003 


2.470 


0.987 


0.383 


0.318 


0.043 


0.323 


0.007 


0.214 


0.005 


0.079 


0.010 


2.470 


1.601 


0.502 


0.345 


0.034 


0.174 


0.003 


0.151 


0.003 


0.060 


0.005 


2.470 


2.832 


0.641 


0.266 


0.069 


0.067 


0.002 


0.072 


0.002 


0.023 


0.004 


2.470 


3.453 


0.685 


0.228 


0.080 


0.045 


0.002 


0.051 


0.002 


0.014 


0.004 


2.510 


0.974 


0.374 


0.443 


0.047 


0.314 


0.007 


0.225 


0.005 


0.104 


0.009 


2.510 


1.609 


0.497 


0.364 


0.033 


0.176 


0.003 


0.154 


0.003 


0.063 


0.005 


2.510 


3.430 


0.678 


0.284 


0.082 


0.046 


0.002 


0.055 


0.002 


0.018 


0.004 


2.550 


1.008 


0.377 


0.481 


0.047 


0.302 


0.007 


0.225 


0.005 


0.109 


0.009 


2.550 


1.594 


0.488 


0.359 


0.033 


0.182 


0.003 


0.158 


0.003 


0.064 


0.005 


2.550 


2.824 


0.628 


0.208 


0.056 


0.074 


0.002 


0.075 


0.002 


0.019 


0.003 


2.590 


0.956 


0.359 


0.486 


0.042 


0.324 


0.007 


0.234 


0.004 


0.113 


0.008 


2.590 


2.764 


0.618 


0.188 


0.059 


0.080 


0.002 


0.079 


0.002 


0.018 


0.004 


2.630 


0.944 


0.351 


0.466 


0.039 


0.338 


0.006 


0.238 


0.004 


0.110 


0.008 



Table 12: Structure functions of carbon. 



202 







X 


R 


AT? 


Fi 




F2 


AF2 


Fl 


AFi 


2.630 


1.538 


0.468 


0.395 


0.037 


0.198 


0.003 


0.172 


0.003 





073 


0.005 


2.630 


3.098 


0.639 


0.286 


0.055 


0.062 


0.001 


0.070 


0.002 





023 


0.003 


2.670 


1.000 


0.359 


0.430 


0.038 


0.332 


0.006 


0.234 


0.004 





102 


0.008 


2.670 


3.076 


0.632 


0.183 


0.048 


0.069 


0.001 


0.070 


0.002 





016 


0.003 


2.710 


0.921 


0.335 


0.494 


0.038 


0.356 


0.006 


0.249 


0.004 





118 


0.007 


2.710 


2.337 


0.561 


0.234 


0.048 


0.116 


0.003 


0.109 


0.003 





030 


0.006 


2.710 


3.054 


0.625 


0.271 


0.051 


0.069 


0.001 


0.076 


0.002 





023 


0.003 


2.750 


0.909 


0.327 


0.431 


0.034 


0.378 


0.006 


0.250 


0.004 





107 


0.007 


2.750 


2.500 


0.572 


0.285 


0.056 


0.101 


0.002 


0.102 


0.003 





033 


0.005 


2.750 


3.032 


0.618 


0.163 


0.045 


0.076 


0.001 


0.076 


0.002 





015 


0.003 


2.790 


0.881 


0.316 


0.498 


0.034 


0.383 


0.006 


0.259 


0.004 





120 


0.006 


2.790 


2.479 


0.565 


0.257 


0.053 


0.108 


0.002 


0.106 


0.003 





031 


0.005 


2.790 


3.009 


0.612 


0.155 


0.045 


0.078 


0.001 


0.077 


0.002 





015 


0.003 


2.830 


0.925 


0.322 


0.492 


0.051 


0.374 


0.009 


0.257 


0.006 





118 


0.010 


2.830 


2.459 


0.558 


0.290 


0.053 


0.113 


0.003 


0.112 


0.003 





036 


0.005 


2.870 


0.882 


0.307 


0.467 


0.049 


0.400 


0.010 


0.262 


0.006 





115 


0.010 


2.870 


2.438 


0.551 


0.317 


0.054 


0.114 


0.003 


0.115 


0.003 





040 


0.005 


2.910 


0.870 


0.300 


0.521 


0.050 


0.399 


0.009 


0.267 


0.006 





125 


0.010 


2.910 


2.418 


0.544 


0.298 


0.051 


0.120 


0.003 


0.119 


0.003 





039 


0.005 


2.950 


0.850 


0.291 


0.483 


0.044 


0.418 


0.009 


0.268 


0.006 





118 


0.009 


2.950 


2.397 


0.537 


0.297 


0.051 


0.124 


0.003 


0.122 


0.003 





040 


0.005 


2.990 


0.839 


0.284 


0.501 


0.044 


0.427 


0.009 


0.272 


0.006 





121 


0.009 


2.990 


2.376 


0.530 


0.299 


0.051 


0.128 


0.003 


0.125 


0.003 





041 


0.005 


3.030 


0.807 


0.273 


0.516 


0.035 


0.447 


0.007 


0.279 


0.004 





126 


0.007 


3.030 


1.781 


0.453 


0.299 


0.068 


0.209 


0.007 


0.175 


0.006 





057 


0.010 


3.030 


2.424 


0.530 


0.253 


0.047 


0.129 


0.003 


0.122 


0.003 





035 


0.005 


3.030 


2.957 


0.579 


0.143 


0.052 


0.098 


0.003 


0.093 


0.003 





016 


0.006 


3.030 


3.584 


0.625 


0.164 


0.075 


0.065 


0.002 


0.069 


0.002 





013 


0.004 


3.070 


0.815 


0.271 


0.613 


0.043 


0.424 


0.008 


0.281 


0.005 





141 


0.008 


3.070 


1.765 


0.446 


0.221 


0.062 


0.223 


0.007 


0.174 


0.006 





044 


0.010 



Table 13: Structure functions of carbon. 



203 







X 


R 


Ai? 


Fi 




F2 


AF2 


Fl 


AFl 


3.070 


2.378 


0.521 


0.315 


0.046 


0.133 


0.002 





130 


0.002 


0.044 


0.004 


3.070 


2.937 


0.573 


0.235 


0.057 


0.096 


0.003 





097 


0.003 


0.026 


0.006 


3.110 


0.813 


0.267 


0.461 


0.025 


0.466 


0.005 





278 


0.003 


0.115 


0.005 


3.110 


1.749 


0.440 


0.220 


0.061 


0.229 


0.007 





177 


0.006 


0.044 


0.010 


3.110 


2.355 


0.514 


0.332 


0.046 


0.138 


0.002 





135 


0.002 


0.047 


0.004 


3.110 


3.642 


0.620 


0.300 


0.088 


0.062 


0.002 





073 


0.003 


0.023 


0.004 


3.150 


0.786 


0.257 


0.607 


0.038 


0.453 


0.008 





289 


0.004 


0.141 


0.007 


3.150 


1.741 


0.434 


0.228 


0.058 


0.233 


0.006 





180 


0.006 


0.046 


0.009 


3.150 


2.356 


0.509 


0.257 


0.046 


0.144 


0.003 





133 


0.003 


0.038 


0.005 


3.150 


3.618 


0.614 


0.268 


0.084 


0.065 


0.002 





074 


0.003 


0.022 


0.004 


3.190 


0.761 


0.248 


0.508 


0.036 


0.494 


0.008 





287 


0.004 


0.124 


0.007 


3.190 


1.593 


0.408 


0.319 


0.063 


0.260 


0.007 





204 


0.006 


0.068 


0.010 


3.190 


2.370 


0.507 


0.256 


0.038 


0.147 


0.002 





135 


0.002 


0.038 


0.003 


3.190 


3.021 


0.567 


0.175 


0.047 


0.100 


0.002 





097 


0.003 


0.020 


0.004 


3.190 


3.614 


0.610 


0.157 


0.073 


0.071 


0.002 





073 


0.002 


0.014 


0.004 


3.230 


0.751 


0.242 


0.591 


0.038 


0.483 


0.008 





292 


0.004 


0.138 


0.007 


3.230 


2.346 


0.500 


0.206 


0.035 


0.155 


0.002 





136 


0.002 


0.032 


0.003 


3.230 


2.999 


0.561 


0.156 


0.045 


0.106 


0.002 





101 


0.003 


0.019 


0.005 


3.230 


3.658 


0.609 


0.211 


0.067 


0.068 


0.002 





074 
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Table 14: Structure functions of carbon. 
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3.394 
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0.263 


0.059 
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0.005 
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0.005 
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0.050 


0.006 
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0.004 
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0.047 


0.007 



Table 15: Structure functions of carbon. 
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Table 16: Structure functions of carbon. 



206 



T/I/2 
vV 


r)2 


X 


It 


A R 
LAJri 






TP 


A P 


jp 

-fL 


A F 


i.iiU 


noR 
z.uyo 


n oni 
u.yui 


n Qoci 
U.oyo 


n 1 ni 
U.iUi 


n nQn 

U.UOU 


n nno 
U.UUz 


n nQi 

U.UOi 


n nno 
U.UUz 


n noi 

U.Uzi 


n nnc; 
U.UUO 


i.iiU 


o. / 66 


n o/io 
u.y4z 


n 07/I 
U.z / 4 


n 1 Q7 
U.io/ 


n nnc 
U.UUo 


n nni 
U.UUi 


n ni 1 
U.Uii 


n nni 
U.UUi 


n nn/i 
U.UU4 


n nno 
U.UUz 


i.ioU 


o n70 

z.u/y 


U.ooo 


U.o04 


n noi 

u.uyi 


n nQO 
U.Uoz 


n nno 
U.UUz 


n nQQ 
U.Uoo 


n nno 
U.UUz 


n non 
U.UzU 


n nni^ 
U.UUO 


i.ioU 


Q 71 /I 
6. 1 i4 


n OQO 
U.yoz 


n 070 
U.z / z 


n 1 QQ 
U.iOO 


n nno 
U.uuy 


n nni 
U.UUi 


n ni o 
U.Uiz 


n nni 
U.UUi 


n nnc; 
U.UUO 


n nno 
U.UUz 


1 1 on 
i.iyu 


Z.UDz 


n SRO 
u.ooy 


n Qi Q 
U.OiO 


n ns7 

U.Uo/ 


n nQci 
U.UOO 


n nno 
U.UUz 


n nQci 
U.UOO 


n nno 
U.UUz 


n ni o 
u.uiy 


n nnc; 
U.UUO 


1 1 on 
i.iyu 


o.oyo 


n ooQ 
u.yzo 


n oon 
u.zyu 


n 1 Qi 

U.iOi 


n nno 
U.uuy 


n nni 
U.UUi 


n ni o 
U.Uiz 


n nni 
U.UUi 


n nnc; 

U.UUO 


n nno 
U.UUz 


i.zoU 


O n/l c: 
Z.U40 


U.O04 


n Qnc 
U.oUo 


n ncQ 
U.Uoo 


n nQQ 
U.Uoo 


n nno 
U.UUz 


n nQQ 
U.Uoo 


n nno 
U.UUz 


n non 
U.UzU 


n nnc; 
U.UUO 


i.z / u 


noQ 

z.Uzo 


n QQO 

U.ooy 


n OKR 

U.ZDD 


n n77 
U.U / / 


n n/i 1 
U.U4i 


n nno 
U.UUz 


n nQO 
U.uoy 


n nno 
U.UUz 


n ni Q 
U.Uio 


n nnc; 
U.UUO 


i.oiU 


1 ooo 
i.yyy 


n soQ 

U.oZO 


n /inn 

U.4UU 


n nsn 

U.UoU 


U.U4U 


n nno 

U.UUZ 


n n/i Q 
U.U40 


n nno 

U.UUZ 


n no7 
u.uz / 


n nnc; 

U.UUO 


1 Qi n 

i.OiU 


o.oyi 


n QOR 

u.oyo 


U.OoD 


n 1 n7 

U.iU/ 


n ni 1 

U.Uii 


n nni 

U.UUi 


n ni 

U.UiO 


n nni 

U.UUi 


n nn7 

0.00 1 


n nno 
U.UUz 


i.ooU 


1 OQ/I 

i.yc54 


n cno 

u.c5uy 


n QnQ 
U.oUo 


n n^jo 
U.Uby 


n n/i 
U.U4D 


n nno 
U.UUz 


n n/i /I 
U.U44 


n nno 
U.UUz 


n noo 
U.Uzz 


n nnc; 
U.UUO 


i.ooU 


O.D/ i 


n QQ7 


n Q7Q 

U.o/ o 


n 1 n^j 
U.iUD 


n ni 1 
U.Uii 


n nni 
U.UUi 


n ni c: 
U.UiO 


n nni 
U.UUi 


n nn7 
U.UU/ 


n nno 
U.UUz 


1 Qon 


nno 
z.uuy 


n 70fi 
U. / yo 


n 07S 
U.z/o 


n n/iR 
U.U4D 


n n/i7 
U.U4 / 


n nni 
U.UUi 


n n/i R 
U.U4D 


n nni 
U.UUi 


n noi 
U.Uzi 


n nnQ 
U.UUo 


i.oyu 


o.u4y 


n 877 
U.o ( ( 


n Q/i/i 

U.044 


n nos 
u.uyo 


n ni 

U.UiZ 


n nni 

U.UUi 


n ni 7 

U.Ui ( 


n nni 

U.UUi 


n nn7 

U.UU ( 


n nno 

U.UUZ 


1 /I Qn 


o ni R 

Z.UiD 


n 7SR 
U. ( oO 


n ooo 
u.zyy 


n nc;o 
U.UOz 


n n/1 fi 
U.U40 


n nni 
U.UUi 


n n/1 7 
U.U4 / 


n nni 
U.UUi 


n noQ 
U.UzO 


n nnQ 
U.UUO 


1 /I Qn 


O.OZO 


U.ooo 


n oc;7 
U.zO / 


n noi 

U.uyi 


n ni Q 
U.Uio 


n nni 
U.UUi 


n ni 7 
U.Ui 1 


n nni 
U.UUi 


n nnfi 
U.UUo 


n nno 
U.UUz 


1 /1 7n 
i.4 / U 


o nnn 
Z.UUU 


n 770 
U. M Z 


n oi /I 
U.zi4 


n n/i fi 
U.U40 


n nc;Q 
U.UOO 


n nni 
U.UUi 


n n/i o 
u.U4y 


n nni 
U.UUi 


n ni fi 
U.Uio 


n nn/1 
U.UU4 


1 A 7n 

i.4 ( U 


Q Rn7 
o.DU ( 


u.ooy 


n Q1 s 
U.Oio 


n no/i 
u.uy4 


n ni /I 
U.Ui4 


n nni 

U.UUi 


n ni s 
U.Uio 


n nni 

U.UUi 


n nnfi 
U.UUO 


n nno 

U.UUZ 


1 c;i n 

i.OiU 


nnQ 
z.UUo 


n 7fii 

U. ( Oi 


n ofifi 
U.zDO 


n n/i/1 

U.U44 


U.U04 


n nni 

U.UUi 


n nc;i 

U.UOi 


n nni 

U.UUi 


n noo 

U.UZZ 


n nnQ 

U.UUO 


1 1 n 

i.iDiU 


o.OoO 


n c 1 
U.oOi 


n 07f; 
U.z /D 


n nfifi 
U.Uoo 


n ni 
U.UiO 


n nni 
U.UUi 


n non 
U.UzU 


n nni 
U.UUi 


n nn7 
U.UU 1 


n nno 
U.UUz 


i.OOU 


i.yo4 


n 7/1 s 
U. ( 40 


n 0/1 Q 
U.Z40 


n n/i 1 
U.U4i 


n nc;7 
U.UO ( 


n nni 
U.UUi 


n nc;Q 
U.UOO 


n nni 
U.UUi 


n noi 
U.Uzi 


n nnQ 
U.UUO 


i.OOU 


O.004 


n S/10 
U.04z 


n QOfi 
U.OzO 


n non 
U.uyu 


n ni 
U.UiO 


n nni 
U.UUi 


n noi 
U.Uzi 


n nni 
U.UUi 


n nno 
U.uuy 


n nno 
U.UUz 


i.oyu 


i.yoz 


n 7QQ 
U. ( 00 


n ooR 
U.zzo 


n nc;i 
U.UOi 


n nfiQ 
U.UOO 


n nno 
U.UUz 


n nc;7 
U.UO ( 


n nno 
U.UUz 


n noi 
U.Uzi 


n nnc; 
U.UUO 


i.oyu 


o.04z 


U.ooo 


n ooc; 

u.zyo 


n ncc; 
U.UoO 


n ni 7 
U.Ui ( 


n nni 
U.UUi 


n noo 
U.Uzz 


n nni 
U.UUi 


n nnc 
U.UUo 


n nno 
U.UUz 












n nfifi 

W . WVJVJ 




n 060 


002 


022 


004 


1.630 


3.521 


0.825 


0.313 


0.084 


0.018 


0.001 


0.023 


0.001 


0.009 


0.002 


1.670 


1.921 


0.709 


0.271 


0.051 


0.068 


0.002 


0.064 


0.002 


0.026 


0.004 


1.670 


3.500 


0.816 


0.465 


0.097 


0.017 


0.001 


0.025 


0.001 


0.013 


0.002 


1.710 


1.905 


0.697 


0.193 


0.047 


0.074 


0.002 


0.065 


0.002 


0.020 


0.005 


1.710 


3.444 


0.806 


0.335 


0.080 


0.020 


0.001 


0.026 


0.001 


0.011 


0.002 


1.750 


1.889 


0.685 


0.310 


0.050 


0.073 


0.002 


0.070 


0.002 


0.031 


0.004 



Table 17: Structure functions of aluminum. 
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0.003 


0.044 


0.006 


2.150 


3.741 


0.747 


0.298 


0.069 


0.028 


0.001 


0.035 


0.001 


0.012 


0.002 


2.190 


3.046 


0.699 


0.323 


0.085 


0.044 


0.001 


0.052 


0.002 


0.020 


0.004 


2.190 


3.718 


0.739 


0.202 


0.063 


0.030 


0.001 


0.035 


0.001 


0.009 


0.002 


2.230 


3.377 


0.714 


0.182 


0.050 


0.039 


0.001 


0.043 


0.001 


0.010 


0.002 


2.270 


1.597 


0.535 


0.284 


0.041 


0.157 


0.003 


0.132 


0.003 


0.048 


0.006 


2.270 


2.999 


0.683 


0.284 


0.079 


0.049 


0.002 


0.055 


0.002 


0.019 


0.004 



Table 18: Structure functions of aluminum. 



208 







X 


R 


AT? 


Fi 




F2 


AF2 


Fl 


AFl 


2.270 


3.684 


0.726 


0.157 


0.065 


0.034 


0.001 


0.038 


0.002 


0.008 


0.003 


2.310 


1.576 


0.524 


0.334 


0.035 


0.159 


0.002 


0.138 


0.002 


0.056 


0.004 


2.310 


2.976 


0.675 


0.253 


0.076 


0.052 


0.002 


0.057 


0.002 


0.018 


0.004 


2.310 


3.663 


0.719 


0.222 


0.068 


0.034 


0.001 


0.040 


0.002 


0.011 


0.003 


2.350 


1.571 


0.517 


0.324 


0.032 


0.165 


0.002 


0.141 


0.002 


0.055 


0.004 


2.350 


2.953 


0.668 


0.235 


0.073 


0.056 


0.002 


0.060 


0.002 


0.017 


0.004 


2.350 


3.670 


0.714 


0.254 


0.078 


0.035 


0.001 


0.042 


0.002 


0.013 


0.003 


2.390 


1.632 


0.520 


0.401 


0.041 


0.154 


0.003 


0.141 


0.003 


0.064 


0.006 


2.390 


2.930 


0.660 


0.359 


0.079 


0.056 


0.002 


0.066 


0.002 


0.027 


0.004 


2.390 


3.648 


0.707 


0.225 


0.075 


0.037 


0.001 


0.043 


0.002 


0.012 


0.003 


2.430 


1.617 


0.511 


0.342 


0.038 


0.164 


0.003 


0.144 


0.003 


0.057 


0.006 


2.430 


2.906 


0.652 


0.271 


0.072 


0.061 


0.002 


0.067 


0.002 


0.021 


0.004 


2.470 


0.962 


0.377 


0.375 


0.046 


0.322 


0.008 


0.220 


0.006 


0.091 


0.010 


2.470 


1.601 


0.502 


0.306 


0.037 


0.174 


0.003 


0.147 


0.003 


0.054 


0.006 


2.470 


2.835 


0.641 


0.271 


0.079 


0.067 


0.002 


0.072 


0.003 


0.023 


0.005 


2.470 


3.453 


0.685 


0.270 


0.090 


0.043 


0.002 


0.050 


0.002 


0.016 


0.004 


2.470 


3.915 


0.711 


0.107 


0.101 


0.036 


0.002 


0.039 


0.002 


0.005 


0.004 


2.510 


0.951 


0.369 


0.413 


0.045 


0.328 


0.008 


0.228 


0.005 


0.100 


0.010 


2.510 


1.609 


0.497 


0.352 


0.036 


0.174 


0.003 


0.151 


0.003 


0.061 


0.005 


2.510 


2.624 


0.617 


0.214 


0.056 


0.083 


0.002 


0.082 


0.002 


0.022 


0.005 


2.510 


3.165 


0.660 


0.208 


0.056 


0.055 


0.002 


0.059 


0.002 


0.015 


0.003 


2.550 


0.983 


0.370 


0.459 


0.047 


0.315 


0.007 


0.229 


0.005 


0.107 


0.010 


2.550 


1.594 


0.488 


0.413 


0.037 


0.174 


0.003 


0.158 


0.003 


0.070 


0.005 


2.550 


2.836 


0.629 


0.343 


0.073 


0.068 


0.002 


0.077 


0.002 


0.029 


0.004 


2.550 


3.984 


0.705 


0.205 


0.115 


0.035 


0.002 


0.041 


0.002 


0.010 


0.004 


2.590 


0.939 


0.354 


0.450 


0.042 


0.341 


0.008 


0.238 


0.005 


0.109 


0.009 


2.590 


2.767 


0.618 


0.195 


0.062 


0.078 


0.002 


0.077 


0.002 


0.019 


0.004 


2.590 


3.960 


0.698 


0.189 


0.111 


0.036 


0.002 


0.042 


0.002 


0.010 


0.004 


2.630 


0.928 


0.347 


0.463 


0.042 


0.346 


0.007 


0.241 


0.005 


0.111 


0.009 


2.630 


1.538 


0.468 


0.389 


0.039 


0.195 


0.004 


0.169 


0.003 


0.071 


0.006 



Table 19: Structure functions of aluminum. 



209 







X 


R 


AT? 


Fi 




F2 


AF2 


Fl 


AFl 


2.630 


3.098 


0.639 


0.235 


0.060 


0.064 


0.002 


0.069 


0.002 


0.019 


0.004 


2.670 


0.977 


0.353 


0.418 


0.039 


0.343 


0.007 


0.237 


0.005 


0.101 


0.009 


2.670 


3.076 


0.632 


0.238 


0.057 


0.065 


0.002 


0.070 


0.002 


0.020 


0.003 


2.710 


0.907 


0.331 


0.491 


0.041 


0.362 


0.007 


0.251 


0.005 


0.118 


0.008 


2.710 


2.337 


0.561 


0.319 


0.056 


0.110 


0.003 


0.110 


0.003 


0.039 


0.006 


2.710 


3.054 


0.625 


0.301 


0.058 


0.066 


0.002 


0.074 


0.002 


0.025 


0.003 


2.750 


0.897 


0.324 


0.513 


0.041 


0.368 


0.007 


0.256 


0.005 


0.122 


0.008 


2.750 


2.500 


0.572 


0.267 


0.060 


0.102 


0.003 


0.101 


0.003 


0.031 


0.005 


2.750 


3.032 


0.618 


0.351 


0.059 


0.069 


0.002 


0.080 


0.002 


0.030 


0.003 


2.790 


0.874 


0.314 


0.523 


0.034 


0.382 


0.006 


0.261 


0.004 


0.125 


0.006 


2.790 


2.479 


0.565 


0.296 


0.061 


0.105 


0.003 


0.106 


0.003 


0.035 


0.006 


2.790 


3.009 


0.612 


0.238 


0.053 


0.076 


0.002 


0.080 


0.002 


0.022 


0.003 


2.790 


3.694 


0.659 


0.183 


0.111 


0.050 


0.002 


0.055 


0.003 


0.012 


0.005 


2.830 


0.909 


0.318 


0.570 


0.045 


0.362 


0.007 


0.259 


0.005 


0.131 


0.008 


2.830 


2.459 


0.558 


0.308 


0.060 


0.109 


0.003 


0.110 


0.003 


0.038 


0.006 


2.830 


3.669 


0.653 


0.154 


0.108 


0.053 


0.002 


0.056 


0.003 


0.011 


0.005 


2.870 


0.870 


0.304 


0.586 


0.045 


0.381 


0.008 


0.268 


0.005 


0.136 


0.008 


2.870 


2.438 


0.551 


0.363 


0.061 


0.111 


0.003 


0.116 


0.003 


0.044 


0.005 


2.910 


0.859 


0.298 


0.578 


0.043 


0.395 


0.008 


0.272 


0.005 


0.136 


0.008 


2.910 


2.418 


0.544 


0.267 


0.055 


0.121 


0.003 


0.116 


0.003 


0.035 


0.006 


2.950 


0.844 


0.290 


0.550 


0.038 


0.409 


0.007 


0.272 


0.004 


0.130 


0.007 


2.950 


2.397 


0.537 


0.228 


0.054 


0.126 


0.003 


0.117 


0.003 


0.031 


0.006 


2.990 


0.833 


0.283 


0.591 


0.039 


0.412 


0.007 


0.278 


0.004 


0.138 


0.007 


2.990 


2.376 


0.530 


0.320 


0.055 


0.127 


0.003 


0.125 


0.003 


0.043 


0.005 


3.030 


0.806 


0.273 


0.593 


0.033 


0.430 


0.006 


0.282 


0.003 


0.139 


0.006 


3.030 


1.781 


0.453 


0.243 


0.063 


0.213 


0.007 


0.171 


0.006 


0.047 


0.011 


3.030 


2.424 


0.530 


0.349 


0.057 


0.124 


0.003 


0.126 


0.003 


0.046 


0.005 


3.030 


2.957 


0.579 


0.254 


0.070 


0.090 


0.003 


0.093 


0.004 


0.026 


0.007 


3.030 


3.584 


0.625 


0.220 


0.088 


0.062 


0.002 


0.069 


0.003 


0.017 


0.004 


3.070 


0.812 


0.271 


0.632 


0.040 


0.426 


0.008 


0.286 


0.005 


0.146 


0.007 



Table 20: Structure functions of aluminum. 



210 







X 


R 


Ai? 


Fi 




F2 


AF2 


Fl 


AFl 


3.070 


1.765 


0.446 


0.284 


0.064 


0.216 


0.007 





177 


0.006 


0.055 


0.011 


3.070 


2.378 


0.521 


0.327 


0.049 


0.133 


0.003 





131 


0.003 


0.045 


0.005 


3.070 


3.559 


0.619 


0.144 


0.083 


0.066 


0.002 





068 


0.003 


0.012 


0.004 


3.110 


0.810 


0.266 


0.520 


0.027 


0.458 


0.006 





284 


0.003 


0.127 


0.005 


3.110 


1.749 


0.440 


0.264 


0.062 


0.224 


0.007 





179 


0.006 


0.052 


0.010 


3.110 


2.355 


0.514 


0.319 


0.051 


0.137 


0.003 





133 


0.003 


0.045 


0.005 


3.110 


3.642 


0.620 


0.154 


0.089 


0.064 


0.002 





066 


0.003 


0.012 


0.005 


3.150 


0.787 


0.257 


0.604 


0.037 


0.454 


0.008 





289 


0.004 


0.141 


0.007 


3.150 


1.741 


0.434 


0.316 


0.062 


0.225 


0.006 





186 


0.006 


0.062 


0.009 


3.150 


2.356 


0.509 


0.230 


0.050 


0.145 


0.003 





131 


0.003 


0.034 


0.006 


3.150 


3.618 


0.614 


0.138 


0.086 


0.067 


0.002 





069 


0.003 


0.011 


0.005 


3.190 


0.766 


0.249 


0.630 


0.037 


0.462 


0.008 





292 


0.004 


0.145 


0.007 


3.190 


1.593 


0.408 


0.313 


0.062 


0.260 


0.008 





203 


0.006 


0.066 


0.011 


3.190 


2.370 


0.507 


0.315 


0.042 


0.144 


0.002 





139 


0.002 


0.046 


0.004 


3.190 


3.021 


0.567 


0.250 


0.058 


0.095 


0.003 





098 


0.003 


0.027 


0.005 


3.190 


3.626 


0.611 


0.225 


0.086 


0.067 


0.002 





073 


0.003 


0.018 


0.004 


3.230 


0.757 


0.244 


0.673 


0.038 


0.463 


0.008 





296 


0.004 


0.152 


0.007 


3.230 


2.346 


0.500 


0.187 


0.039 


0.154 


0.002 





133 


0.002 


0.029 


0.004 


3.230 


2.999 


0.561 


0.254 


0.056 


0.099 


0.003 





102 


0.003 


0.028 


0.005 


3.230 


3.656 


0.609 


0.149 


0.080 


0.068 


0.002 





070 


0.002 


0.012 


0.004 


3.270 


0.726 


0.233 


0.604 


0.029 


0.500 


0.006 





296 


0.003 


0.141 


0.005 


3.270 


1.686 


0.414 


0.342 


0.063 


0.243 


0.007 





199 


0.006 


0.069 


0.010 


3.270 


2.288 


0.489 


0.359 


0.052 


0.155 


0.003 





151 


0.003 


0.055 


0.006 


3.270 


2.977 


0.555 


0.118 


0.067 


0.109 


0.004 





099 


0.004 


0.014 


0.007 


3.270 


3.631 


0.603 


0.237 


0.084 


0.069 


0.002 





076 


0.003 


0.020 


0.004 


3.310 


1.670 


0.407 


0.274 


0.058 


0.257 


0.007 





198 


0.006 


0.057 


0.010 


3.310 


2.266 


0.482 


0.285 


0.049 


0.163 


0.003 





148 


0.003 


0.045 


0.006 


3.310 


2.974 


0.550 


0.225 


0.072 


0.105 


0.004 





104 


0.004 


0.026 


0.007 


3.310 


3.606 


0.598 


0.156 


0.079 


0.073 


0.002 





075 


0.003 


0.014 


0.004 


3.350 


1.654 


0.401 


0.376 


0.063 


0.251 


0.007 





207 


0.006 


0.076 


0.010 



Table 21: 



Structure functions of aluminum. 



211 







X 


R 


AT? 


Fi 




F2 


AF2 


Fl 


AFl 


3.350 


3.015 


0.550 


0.353 


0.078 


0.099 


0.004 


0.109 


0.004 


0.038 


0.007 


3.350 


3.585 


0.592 


0.333 


0.093 


0.072 


0.002 


0.084 


0.003 


0.028 


0.005 


3.390 


2.220 


0.469 


0.372 


0.046 


0.172 


0.003 


0.164 


0.003 


0.060 


0.005 


3.390 


2.792 


0.527 


0.218 


0.062 


0.122 


0.004 


0.116 


0.004 


0.028 


0.007 


3.390 


3.365 


0.573 


0.254 


0.087 


0.085 


0.003 


0.091 


0.003 


0.025 


0.005 


3.430 


2.198 


0.463 


0.322 


0.045 


0.180 


0.003 


0.164 


0.003 


0.053 


0.005 


3.430 


2.712 


0.515 


0.305 


0.070 


0.127 


0.004 


0.127 


0.005 


0.040 


0.008 


3.430 


3.306 


0.565 


0.243 


0.085 


0.090 


0.003 


0.095 


0.003 


0.025 


0.005 


3.470 


2.175 


0.456 


0.290 


0.043 


0.187 


0.003 


0.165 


0.003 


0.050 


0.005 


3.470 


2.692 


0.510 


0.215 


0.063 


0.137 


0.004 


0.127 


0.005 


0.030 


0.008 


3.510 


1.591 


0.377 


0.333 


0.057 


0.287 


0.007 


0.220 


0.006 


0.072 


0.009 


3.510 


2.153 


0.450 


0.318 


0.044 


0.191 


0.003 


0.170 


0.003 


0.055 


0.005 


3.510 


2.927 


0.527 


0.173 


0.070 


0.123 


0.004 


0.114 


0.005 


0.022 


0.008 


3.510 


3.488 


0.570 


0.191 


0.079 


0.087 


0.003 


0.089 


0.003 


0.019 


0.005 


3.550 


1.575 


0.371 


0.308 


0.055 


0.300 


0.007 


0.223 


0.006 


0.069 


0.009 


3.550 


2.130 


0.444 


0.289 


0.043 


0.199 


0.003 


0.172 


0.003 


0.051 


0.005 


3.550 


2.657 


0.499 


0.189 


0.069 


0.145 


0.005 


0.129 


0.005 


0.027 


0.009 


3.550 


3.267 


0.550 


0.282 


0.085 


0.097 


0.003 


0.103 


0.003 


0.030 


0.005 


3.590 


1.719 


0.388 


0.304 


0.054 


0.272 


0.006 


0.210 


0.005 


0.064 


0.009 


3.590 


2.342 


0.464 


0.298 


0.045 


0.174 


0.003 


0.158 


0.003 


0.048 


0.005 


3.630 


1.618 


0.370 


0.332 


0.052 


0.291 


0.006 


0.221 


0.005 


0.072 


0.008 


3.630 


2.320 


0.458 


0.352 


0.046 


0.175 


0.003 


0.164 


0.003 


0.056 


0.005 


3.670 


2.238 


0.445 


0.362 


0.036 


0.191 


0.003 


0.177 


0.003 


0.062 


0.005 


3.710 


2.175 


0.435 


0.317 


0.035 


0.204 


0.004 


0.179 


0.003 


0.056 


0.005 


3.710 


2.804 


0.498 


0.180 


0.059 


0.144 


0.004 


0.129 


0.004 


0.026 


0.007 


3.710 


3.383 


0.544 


0.225 


0.060 


0.100 


0.002 


0.102 


0.003 


0.024 


0.005 


3.750 


2.175 


0.431 


0.287 


0.043 


0.210 


0.005 


0.179 


0.004 


0.052 


0.007 


3.750 


2.760 


0.490 


0.312 


0.072 


0.142 


0.005 


0.140 


0.005 


0.043 


0.008 


3.750 


3.352 


0.539 


0.222 


0.056 


0.103 


0.002 


0.104 


0.003 


0.025 


0.004 


3.790 


2.154 


0.425 


0.339 


0.045 


0.212 


0.005 


0.186 


0.004 


0.061 


0.007 



Table 22: Structure functions of aluminum. 
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vV 
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Table 23: Structure functions of aluminum. 
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Table 24: Structure functions of iron. 
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Table 26: Structure functions of copper. 
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Table 27: Structure functions of copper. 
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Table 28: Structure functions of copper. 



HMS Kinematic Settings 



HMS Central Momentum P (GeV) 


HMS central angle 9 
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0.707 
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0.813 


75.00 



Table 29: Beam energy is 2.097 GeV. 



219 



HMS npnfrfll MnmpTifiim P ITIpVI 

±±iVJ.kJ y^Kjllbl CXil iVJ-UlllCilL Llili L I VJC V J 


T-TA/rS ppnfrfll flncrlp ft 


1 184 


1 nn 




1 nn 


1 n9Q 


95 00 


1 184 




1 '^Rl 

1.0L)± 


O'^ 00 






n 77R 


'^'^ 00 


yJ.OiJO 


'^'^ 00 


1 090 


00 




3^ 00 
oo.yjyj 


1 '^fil 


35 00 


n fi78 


45 00 

'-tO.yJyJ 


n 77R 


oo.yjyj 




oo.yjyj 


X . UZjC/ 


d^ nn 

'TO.yJyJ 


1 1R4 

X . XOt: 


d^ nn 


X . OU X 


d^ nn 




nn 

TO.yJyJ 




Rn nn 


U. OO X 


Rn nn 


n R7R 


Rn nn 

uu. uu 


0.778 


59.99 


0.895 


60.00 


0.439 


74.99 


0.439 


74.99 


0.505 


74.99 


0.581 


75.00 


0.678 


75.00 


0.778 


75.00 



Table 30: Beam energy is 3.116 GeV. 
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HMS Central Momentum P (GeV) 


HMS central angle 9 
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1.148 
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1.318 


48.00 


1.318 


48.00 


1.318 


39.00 



Table 31: Beam energy is 4.134 GeV. 
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HMS Central Momentum P (GeV) 


HMS central angle 6 
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Table 32: Beam energy is 5.151 GeV. 
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1.3 Cross Sections 



E = 2.10GeV,e = 45 
= 0.5-0.9 GeV^ 




0.5 1 1.5 2 2.5 3 3.5 4 



(GeV^) 

Figure 89: Extracted differential cross section for carbon compared to the model 
cross section. 
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Figure 90: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 91: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 92: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 93: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 94: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 95: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 96: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 97: Extracted differential cross section for carbon compared to tlie model 
cross section. 
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Figure 98: Extracted differential cross section for iron compared to tlie model cross 
section. 
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Figure 99: Extracted differential cross section for iron compared to tlie model cross 
section. 
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Figure 100: Extracted differential cross section for iron compared to tlie model cross 
section. 
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Figure 101: Extracted differential cross section for iron compared to tlie model cross 
section. 
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Figure 102: Extracted differential cross section for iron compared to tlie model cross 
section. 
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Figure 103: Extracted differential cross section for copper compared to the model 
cross section. 
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Figure 104: Extracted differential cross section for copper compared to the model 
cross section. 
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Figure 105: Extracted differential cross section for copper compared to the model 
cross section. 



239 



600 
400 





400 
^ 300 
CD 200 

O 

^ 100 

(/) 

~o 



200- 



100 — 


100 



- E = 5.15GeV,e = 27 

- Q2 = 2.2-3.8 GeV^ 















































1 1 1 1 1 1 1 1 1 

: Old model 

: New model 

: • Data 
































: E = 5.15GeV,e = 30 
: = 2.3-4.1 GeV^ 































; E = 5.15GeV,e = 33 
- Q2 = 2.8-4.2 GeV^ 




























, ,...ri 


, 










: E = 5.15GeV,e = 38 
- Q2 = 3.2-4.2 GeV^ 




























— r— r i i 


1 r--i 1 








r-v 



2 3 4 

(GeV^) 



Figure 106: Extracted differential cross section for copper compared to the model 
cross section. 



